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I. INTRODUCTION 



At a first glance it is surprising that a change in the shock-normal angle 0Bn by just a few degrees from, say, 

©Bn = 50° to Qbu ~ 40° should completely change the character of the supercritical shock. We have seen in the 
Q last chapter, when having discussed the conditions for reflection of ions from the shock ramp that such a change in 

the shock properties is theoretically predicted. The critical shock-normal angle Qbu — 45° does indeed separate two 
c/3 completely different phases of a supercritical shock. At this angle the shock experiences a 'phase transition' from 
I ^ I the quasi-perpendicular to the quasi-parallel shock state, with the shock-normal angle Qbu having the property of a 

'critical control parameter'. 

T— I We have already learned that the reason for the different behaviour of the two shock phases is that in a quasi- 

^ perpendicular shock all reflected particles from the foot region of the shock return to the shock after not more than a 
few gyrations when they have picked up sufficient energy in the upstream convection electric field Ei = — Vi x Bi to 
ultimately overcome the shock ramp potential, pass the shock ramp and to merge into the downstream flow. We have 
not discussed what happens to the accelerated ions in the downstream region as this is not of primary importance 
in the shock formation mechanism which to good approximation depends only on the upstream conditions. This 
question will be treated in a separate section on shock-particle acceleration. 

In contrast, in a quasi-parallel shock the combined geometries of the upstream magnetic field and generally curved 
shock surface prevent the shock-reflected particles from immediate return to the shock. The reason is that their 
gyro-orbits, after having suffered reflection from the shock ramp, lie completely upstream, outside the shock ramp, 
such that they do not touch the shock ramp again after reflection. Since, in addition, their upstream velocities have a 
large component parallel to the upstream magnetic field, which increases the more the shock-normal turns parallel to 
^ the upstream magnetic field, the refiected particles are enabled to escape upstream from the shock along the magnetic 
^ field thereby forming fast upstream particle beams. 

A quasi-parallel supercritical collisionless shock thus populates the upstream space with a refiected particle compo- 
nent. This population moves a long distance away from the shock along the magnetic field. Was the upstream fiow, 
in the case of the quasi-perpendicular shock, completely uninformed about the presence of the shock up to a distance 
of the mere width of the shock foot so, at the quasi-parallel shock, it receives a first signal of the presence of a shock 
already at quite a large upstream distance when the first and fastest refiected particles arrive on the magnetic field 
lines that connect the fiow to the shock. It is, in fact, only these particles that can inform the fiow about the presence 
of a supercritical shock, because any low-frequency plasma wave cannot propagate far upstream for supercritical Mach 
numbers M > Mc, while any electromagnetic radiation that is generated at the shock has frequency LUrad > ^pe- For 
it the plasma fiow presents a vacuum. The upstream fiow recognises the refiected particles in its own frame of reference 
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FIG. 1: Schematic one-dimensional profile taken along the nominal instantaneous shock normal of a supercritical quasi-parallel shock 
as seen in the magnetic field component Bz. This is the analogue to the quasi-perpendicular shock profile. It shows the main features in 
the vicinity of the quasi-parallel shock transition: the large amplitude upstream waves with the turbulent fluctuations on top of the 
waves, the formation of shocklets, i.e. steep flank formation on the waves exhibiting small-scale fluctuations on top of the wave, which 
act already like small shocks, very-large amplitude pulsations (magnetic pulsations or SLAMS) which turn out to be the building blocks 
of the shock, multiple shock-ramps at the leading edges of the pulsations belonging to diverse ramp-like steep transitions from upstream 
to downstream lacking a clear localisation of the shock transition (Note that the entire figure is, in fact, the shock transition, as on this 
scale no clear decision can be made where the shock ramp is located.), and their attached phase-locked whistlers. Not shown here are the 
out of plane oscillations of the magnetic field that accompany the waves. Also not shown is the particle phase space. 



as a high-speed magnetic-field ahgned beam. Thereby a beam-beam configuration is created which leads to a number 
of beam-driven instabilities. These excite various plasma waves that fill the space in front of the shock and modify 
its properties. 

Figurejl] shows a sketch of the magnetic profile of a supercritical quasi-parallel shock which contrasts the profile of 
a quasi-perpendicular shock that has been given in the previous chapter. The quasi-parallel magnetic shock profile 
is much stronger distorted than that of a quasi-perpendicular shock, such that it becomes difficult to identify the 
location of the genuine shock ramp on the profile. 

The main difference between quasi-perpendicular and quasi-parallel shocks is that quasi-perpendicular shocks pos- 
sess a narrow ~ Ird wide foot region that is tangential to the shock surface, while quasi-parallel shocks possess an 
extended foreshock region. Interestingly, in curved shocks which arise, for instance, in front of spatially confined 
obstacles, both phases of a supercritical shock can co-exists at the same time, being spatially adjacent to each other. 
An example is shown in Figure ?? in the sketch of the curved Earth's bow shock. Dealing with quasi-parallel shocks 
means to a large extent dealing with the processes that are going on in the foreshock. It will thus be quite natural 
to start with a discussion of the properties of the foreshock. This discussion will occupy a substantial part of this 
chapter. 

However, before continuing we point out that in spite of the strict distinction between the quasi-perpendicular 
and quasi-parallel shocks there is also a close relation between the two. Both, being supercritical, can exist only 
because they reflect ions; and both possess an upstream region in front of the shock transition that is populated by 
the reflected ions. That this region is narrow in the case of quasi-perpendicular shocks is a question of the ions being 
tied to the magnetic field, which also holds in the case of the quasi-parallel shock. At quasi-perpendicular shocks the 
reflected ions do readily return to the shock. At the quasi-parallel shock they ultimately do also return to the shock, 
but only after having been processed far away from the shock in the foreshock, having coupled to the flow, having 
passed several stages in this processing, and having become the energetic component of the main flow. As such they 
finally arrive at the shock together with the stream. In between their main duty was to dissipate the excess energy, 
which they possessed when arriving for the first time at the shock and which could not be dissipated in the narrow 
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FIG. 2: Comparison of average bulk plasma parameters in quasi-perpendicular and quasi-parallel shocks [after [T6] . The figure shows 
(on the left) the mean magnetic field B, bulk flow velocity V, average magnetic fluctuation amplitude at ultra-low frequencies, and (on 
the right) plasma density N, density of high energy ions Nui, and electron temperature Te. The shaded regions are the downstream 
parts of the AMPTE IRM crossings of the bow shock. The data have been obtained by normalising the time with respect to crossing the 
nominal shock ramp by using the measured normal component (not shown) of the bulk flow velocity. This for the many observations 
included in this figure implies a stretching (or squeezing) of each individual shock crossing, causing some uncertainty, in particular for 
quasi-parallel shock crossings as there the shock ramp is not well defined. However, this figure serves for an immediate overview of the 

differences in both types of shocks. 



shock-ramp transition region. This could been achieved only in the broad extended foreshock which, seen from this 
point of view, is already the shock. It belongs inextricably to the quasi-parallel shock transition. Here, a substantial 
fraction of the energy of the incident flow is dissipated in a way which is completely different from the flow being 
shocked. 

It is these dissipation processes that cause the main difference between the quasi-perpendicular and quasi-parallel 
states of a coUisionless supercritical shock. To stress the analogy with phase transitions a little further, we may say 
that quasi-perpendicular shocks are in the solid - or ordered - shock state, while quasi-parallel shocks are in the fluid 
- or partially disordered - shock state. 



II. THE (QUASI-PARALLEL SHOCK) FORESHOCK 

Quasi-parallel shocks are abundant in space because in most cases when shocks develop the flow direction is inde- 
pendent of the direction of the magnetic field. Moreover, as pointed out earlier, when a bow shock forms around an 
obstacle (planet, magnetosphere, moon ...) this bow shock is curved around the obstacle, and the shock is quasi- 
perpendicular only in a certain region on the shock surface that is centred around the point where the upstream 
magnetic field touches the shock tangentially. Farther away the shock turns gradually to become quasi-parallel. On 
the other hand, when a shock survives over a long distance in space as for instance in supernova remnants then it 
sweeps the upstream magnetic field and pushes it to become more tangential to the shock surface. In this case the 
shock is about quasi-perpendicular. We will later provide arguments that any quasi-parallel supercritical shock on 
the small scale close to the shock surface, i.e. on the electron scale, behaves quasi-perpendicularly while on the larger 
ion scale it remains to be quasi-parallel. This has consequences for the differences in the dynamics of electrons and 
ions during their interaction with quasi-parallel shocks. 

Figure|2] shows at one glance the main differences in the (average) bulk plasma parameters between quasi- 
perpendicular and quasi-parallel shocks as measured with the AMPTE IRM spacecraft at many crossings of the bow 
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shock. The data used in this figure have been stapled, averaged and plotted with respect to the time normalised to 
the shock ramp crossing. For such a normalisation one uses the shock-normal upstream velocity to recalculate the 
time. This procedure is not very certain for quasi-parallel shocks since - as we will see later - the shock ramp is ill 
defined in a quasi-parallel shock. However, for a simple comparison of the main differences this uncertainty is less 
severe. 

The shaded area in the figure corresponds to the downstream region. Shown are - in pairs of quasi- 
perpendicular/quasi-parallel values - the magnetic field B, bulk velocity V, average fiuctuation amplitude in the 
ultra-low frequency waves |b|, plasma density N, high-energy ion density Nih of energy > 15keV, and electron tem- 
perature Te. The general conclusion from this figure is that all quantities in the quasi-perpendicular case exhibit a 
much sharper transition than in the quasi-parallel case. Moreover, the quasi-perpendicular averages are quieter than 
those of the quasi-parallel case. Also, in general, the quasi-parallel levels are higher than the quasi-perpendicular. In 
almost all cases the pre-shock levels are enhanced in the quasi-parallel shock case with over the pre-shock levels of 
quasi-perpendicular shocks. This is seen most impressively in the energetic ion density, which is nearly constant over 
this distance/time scale at quasi-parallel shocks and much higher than that in quasi-perpendicular shocks, signifying 
on the one hand the importance of energetic particles in quasi-parallel shock dynamics, on the other hand the capa- 
bility of quasi-parallel shocks to accelerate particles to substantial energies. The presence of energetic ions (particles) 
far in front of the quasi-parallel shocks and the enhanced pre-shock levels indicate the importance of foreshocks in 
quasi-parallel shock dynamics. In the following we will therefore first concentrate on the foreshock. 

The physics of quasi-parallel shocks cannot be understood without reference to the foreshock. The foreshock is that 
part of the upstream shock region that is occupied with refiected particles. At a curved shock, like the Earth's bow 
shock, the foreshock starts on the shock surface at the location where the upstream magnetic field shock-normal angle 
exceeds Os„ ^ 45°. From that point on electrons and ions escape along the magnetic field in upstream direction. 
Since electrons generally move at a larger parallel velocity than ions they are less vulnerable to the convective motion 
of the upstream magnetic field line to which they are tied, and so there is generally a region closer to the foreshock- 
boundary magnetic field line where only upstream electrons are found. This region is confined approximately between 
the line that marks the electron foreshock boundary and the more inclined line line that marks the ion foreshock 
boundary. An example of this geometry was depicted in Figure?? for the bow shock of the Earth. 

More schematically this is shown in a simplified version in Figurc[3] for the particular case that the upstream 
magnetic field forms an angle of 45° with the symmetry axis of the shock. In this case half of the shock is quasi- 
perpendicular and the other half is quasi-parallel. The figure also shows the directions of three shock normals, the 
narrow foot region in front of the quasi-perpendicular shock, and the two (electron and ion) foreshocks. Particles 
escape from the quasi-parallel shock along the upstream magnetic field. The magnetic field is convected toward the 
shock by the perpendicular upstream velocity component Vj^ as shown in the figure. This component adds to the 
velocity of the upstream particles leading to an inclined foreshock boundary. Since the ions have much smaller speed 
than the electrons, the ion foreshock boundary is more inclined than the electron foreshock boundary. 

In discussing the properties of the foreshock one thus has to distinguish of which foreshock is the talk. However, the 
properties of the electron foreshock are not as decisive for the formation of a quasi-parallel shock as are the properties 
of the ion foreshock. Because of this reason we will, in the following, refer to the ion foreshock as the foreshock. The 
electron foreshock properties we will mention only later. 

A. Ion foreshock 

The ion foreshock is not a homogeneous and uniform region. The refiected ion component evolves across the ion 
foreshock from the ion foreshock boundary to the centre of the ion foreshock and from there towards the shock. 
Speaking of a refiected ion component that can unambiguously identified as being refiected, i.e. streaming into the 
upward direction, makes sense only in the immediate vicinity of the ion foreshock boundary. Here the refiected ions 
appear as a fast ion beam the source of which can be traced back to the shock. Deeper in the foreshock the beam 
component cannot be identified anymore. 

The ion foreshock boundary region 

First identifications of refiected beam protons in space in the magnetic fiux tube connected to the Earth's bow shock 
wave were reported by Gosling et al. ,40] and Paschmann et al. |96j who distinguished those beams by their poorly 
resolved distribution functions from more diffuse protons deeper in the foreshock. Interestingly, observations in the 
foreshock of interplanetary travelling shocks did not show any indication of such beams but only the diffuse ion 
component. Figure|4] gives an observational example of such a reflected ion beam that propagates very close to the 
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FIG. 3: Schematic of the the relation between a curved shock and its foreshock in dependence on the direction of the upstream 
magnetic field B, shock-normal n, and shock-normal angle ©s„ for the special case when the magnetic field is inclined at 45° with 
resepct to the symmetry axis of the shock. In this case the upper half of the shock becomes quasi- parallel {0Bn < 45°), the lower half is 
quasi-perpendicular (0s„ > 45°). The velocity of reflected particles is along the magnetic field. However, seeing the flow the fleld-line to 
which they are attached displaces with perpendicular velocity. This velocity shifts the foreshock boundary toward the shock as shown for 
electrons (light shading) and ions (darker shading). The ion foreshock is closer to the shock because of the lower velocity of the ions than 
the electrons. For the electrons the displacement of the electron foreshock boundary is felt only at large distances from the shock. 



1S57<»ia59!00 UT ISEE 1 Inbound 19 Nov 1977 




FIG. 4: ISEE 1 observation of a reflected ion beam on November 19, 1977 propagating along a magnetic fleld line that was connected 

to the quasi-parallel Earth's bow shock. Left: A two-minute average pseudo-threedimensional ion velocity-space profile in the 
{vx,Vy)-plane showing the undisturbed and cold (narrow) plasma infiow in negative i^^i-direction, and the fast and warm (broad) beam of 
refiected ions propagating in positive u^^-direction and spreading in Vy. This beam is quite anisotropic in temperature. Velocities are in 
kms~^. The scale on the right is count rates, and background count rates were suppressed by choosing only values above 50 s~^. Right: 
Contour plot of a similar beam a little earlier showing that the beam is centred on the magnetic field that connects to the shock, is quite 
narrow along the field and about 2-3 times as broad perpendicular to the field [after 96.. The cross indicates the origin (zero velocity), 
the dot the bulk flow centre. The 10~^^ s"^ cm~® level flux contour has been marked. 



foreshock boundary upstream away from the shock. The bulk flow is the narrow cold beam in the left part of the 
figure which is displaced in negative i;a;-direction (note that in this figure the positive direction points away from the 



6 




74.5' 



CLUSTER FOM 
jl-3-20(t] 



1S:45 



dowrifitreatTi 



TS:4S 




18:54 



Clusler SC I CO^IF H ' 31-3 -2001 

2WKIJ ; loot 





1«D -lIHg -1D« »« 
V [hnilEn] 



FIG. 5: Ion phase space at three locations along the shock-connected magnetic field line in a high-Mach number supercritical 
quasi-perpendicular shock the magnetic field of which is shown in the top panel. The shock-normal angle is ©s„ = 74.5°. The three 
lower plots show downstream, shock ramp/shock foot, and distant upstream phase-space plots. The ramp/foot plot shows the presence 
of the inconing flow (SW) and gyrating ions in the foot. The upstream plot shows the usptream field-aligned relatively hot (large 

velocity spread) beam well separated from the inflow (SW) [from [64] . 



shock). The reflected beam is less dense (lower count rates) but much more energetic. It is displaced in -l-Ujj-direction, 
i.e. streaming away from the shock, and has also a —w,; -component, i.e. it constitutes a gyrating bunch of ions moving 
away from the shock. In the right part of the figure it is seen that the beam is moving away along the magnetic field 
line that is connected to the shock, while the bulk of the plasma fiows in positive direction. 

These beams along the foreshock boundary play some role in the foreshock dynamics as they seem to represent a 
source population for the entire ion foreshock. Whether and why this is really so is not yet been fully understood 
as the shock should refiect ions at almost every place in its quasi-parallel state. However, it seems as that only the 
group of ions that escape from the shock along the foreshock boundary can form such beams. This points on a further 
interesting relation between quasi-perpendicular and quasi-parallel shocks at a curved shock surface with a smooth 
transition from quasi-perpendicular to quasi-parallel as sketched in FigurejS] and realised in space, for instance, at 
planetary bow-shocks. It seems as so these beams escape from the quasi-perpendicular region of the shock along the 
nearly tangential field lines. This would also be in agreement with the observation \41} that the foreshocks of extended 
interplanetary shocks do not show any signs of reflected ion beams. They are only very weakly curved being nearly 
planar, and not possessing a recognisable quasi-perpendicular area on the surface. 

[64i analysed ion distributions along magnetic field lines that were connected to the quasi-perpendicular area of 
Earth's bow shock. Investigating the origin of those beams these authors found that, indeed, the observed ion beams 
at the foreshock boundary result from reflection at the quasi-perpendicular shock and that, without the presence of 
a quasi-perpendicular region, there would presumably be no distinct foreshock boundary and no ion beams escaping 
into the foreshock. If this is really the case, then the population in the quasi-parallel ion foreshock is indeed provided 
by two different ion sources, the beams from the quasi-perpendicular shock region and the genuine foreshock ion 
population. The latter has no beam character but is rather a diffuse ion component [see also i82) . 

estimate that roughly 2% of the ion inflow leaves the shock ramp upstream in the form of a beam along 
the magnetic field. They argue that the ions which escape along the magnetic field, are reflected from the very 
ramp/overshoot region where they have been in resonance with low-frequency plasma waves, which they assume 
to be large amplitude Alfven-whistler waves. These ions experience pitch-angle scattering and pitch-angle diffusion 
towards small pitch angles, and subsequently can escape along the magnetic field in the upstream direction. Since 
the conditions for escape depend in the first place on the pitch-angle scattering process, the beams should be highly 
variable in time and location. The mechanism might still sound a bit speculative as long as no simulation proves its 
reality, but any mechanism which is able to pitch-angle scatter ions along the magnetic field in a quasi-perpendicular 
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FIG. 6: Left: The reduced field-aligned ion distribution functions in the foreshock region showing the field-aligned refiected ion beam 
(FAB) at 1220-1225 UT on February 18, 2003 which arrives at the spacecraft location along the magnetic field line that is connected to 

the quasi-perpendicular region of the Mach number Ma ^ 8 supercritical shock (smooth solid line). The upstream bulk velocity was 
Vi ~ 660kms~^. Also shown is the reduced parallel distribution function (dashed) for the diffuse ion distribution observed at 1130-1135 

UT the same day (for the geometry see FigurcjTji. It is clearly seen that the latter is about symmetric to the magnetic field direction 
indicating the symmetry of the ring of diffuse ions, while the foreshock-boundary field-aligned beam is flowing in the direction away from 
the shock into the upstream medium (negative velocities). Note also that the diffuse distribution appears as a smooth tail on the full ion 
distribution. The small gap on the left is uncertain as it dips into the 1-count level. Right: The two-dimensional phase space plot for the 
time interval 1220-1225 UT when the field-alighned beam was observed. Indicated are the bulk upstream flow, fleld-ligned ion beam, and 

the gyro-phase bunched residue of the diffuse upstream ions [after l57l . 

shock will naturally cause ion beams to escape from the ramp both in the upstream and in the downstream directions. 
Such simulations require a three-dimensional treatment which is not in reach yet. 

A measured example of the ion foreshock-boundary field-aligned ion beam distribution [57 is shown in Figure[6]for 
an upstream flow velocity Vi ~ 660 km s^^, Mach number A4a = 8, and average shock-normal angle Ob„ ~ 15° [2], 
at an upstream distance from the shock (in this case again the Earth's bow shock). This distribution is a so-called 
reduced distribution; it is the integrated over pitch-angle (p Etnd perpendicular velocity v± magnetic field-aligned 
phase space distribution function /(wy) ~ 27r J v'j_dv±f{v±,v^\) that has been appropriately binned and smoothed. 
The information to be taken out of this figure is that the reduced ion-beam distribution (solid line) is narrow in 
velocity, maximising at a speed (in absolute terms \vt\\\ — 800km/s) that is only slightly larger than the parallel 
flow velocity (Vni ~ 640km/s, when taking into account the shock-normal angle), i.e. \vb\\ \ — 1.25 Viy. It is directed 
opposite to the flow. For comparison a reduced parallel diffuse ion distribution is shown in the same plot taken deeper 
in the foreshock. This distribution is about symmetric to the magnetic fleld direction, indicating the about circular 
phase-space distribution of the diffuse ion component which appears as an energetic tail on the main ion distribution 
(note that the small gap on the left of the dashed diffuse-ion distribution curve is questionable as it dips below the 
1-count level). [82 reported CLUSTER observations of foreshock-boundary ion beams simultaneously with diffuse ions. 
They found that the nominal ion-beam velocity Vb — 1.7 Vi had no relation to any known shock-reflection mechanism 
like specular reflection, a conclusion which supports pitch-angle scattering as the beam-injection process as this is 
independent of the speed of the inflow. 

|57j have carefully analysed the relation between the observed ion distribution at the upstream spacecraft position 
and the calculated shock-normal angle &Bn, determined from the local upstream magnetic fleld direction and the 
predicted shape of the bow shock. The result is shown in Figure[7]for three successive times on February 18, 2003 when 
the CLUSTER spacecraft was outside the bow shock. At 1130 UT and 1355 UT the shock normal angles Qbu ^ 15° 
and shock-spacecraft distances ~ 7.5 resp. ~ 6.7 Re were similar. Despite this similarity the observed ion phase- 
space distributions were completely different. At 1130UT no low-energy gyrating ions were observed, while they 
were present at the later time 1355UT. Hence, at 1130 UT the spacecraft must have been closer to the ion foreshock 
boundary, i.e. the foreshock boundary was more inclined than at the later time such that beam particles scattered 
from the foreshock boundary have not arrived at the location of the spacecraft. Due to the velocity fllter effect they 
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FIG. 7: Reconstruction of the shape and location of the ion-foreshock boundary from the measured upstream plasma properties 
(direction of magnetic field, speed and density) for three successive times in order to explain the ion-phase space observations in Figure[6] 

[from [57] . The cold magnetic field-aligned ion beam is observed when the spacecraft is located in the vicinity of the ion-foreshock 
boundary at 1220 UT, while at 1130 UT and 1355 UT no beam was detected. It is assumed that the beam is generated along the 
foreshock boundary at the position where at the quasi-perpendicular shock surface the shock-normal angle is roughly about Qbu ~ 60°. 

have been separated out at spacecraft distance. The reconstruction of the position of the ion-foreshock boundary using 
the measured upstream conditions and shape of the shock for this period is shown in Figure[7] Indeed, at 1220 UT 
the spacecraft was close to the foreshock boundary and, as expected, detected the ion beam (as seen in Figure|6|. 

This observation supports the above advocated view that the ion-foreshock beam is generated in the transition 
region from quasi-perpendicular to quasi-parallel shock. The remaining questions to be answered are: what pitch- 
angle scattering mechanism is responsible for the generation of such a beam, and what is the fate of the foreshock- 
boundary beam-ions during their propagation along the foreshock boundary? Do they contribute to the foreshock ion 
population and if, in what way? Currently we are not able to answer either of these questions definitely. In particular, 
the pitch-angle scattering mechanism is unknown or at least uncertain. When discussing wave generation, we will 
touch on the problem of the fate of the beam. Below we present evidence for the scattering of the ion-beam ions and 
merging into the upstream foreshock diffuse-ion population. 

Diffuse ions 

Reflected upstream magnetic field-aligned ion beams are observed at the foreshock boundary only. The second (and 
main) ion component encountered in the foreshock is the diffuse ion population which is detected there as the energetic 
extension of the inflow plasma. It is widely assumed that the origin of this component is also at the shock as there is 
no other energetic particle source available. However, close investigation of the diffuse ion component in the foreshock 
of the Earth's bow shock wave has demonstrated that these ions are not produced in a specular reflection process at 
the shock. Rather their origin is of diffusive nature. 

So far it has not been possible to identify the source of these ions though some models have been proposed which 
we will briefly discuss below in relation to the appearance of foreshock waves. We will also return to these particles in 
the next chapter on shock particle acceleration. Here, we merely discuss some of their properties and provide evidence 
that they are indeed coming from the shock, proving that the shock is an energetic ion source. 

The foreshock-boundary ion beams also merge into the foreshock particle distribution by scattering from the fore- 
shock boundary off their self-generated wave spectrum and subsequent being convected downstream by the bulk flow. 
[96J have mapped this merging process by following the evolution of the foreshock-boundary beam distribution shown 
in Figure|4] during the convection. In that flgure the beam was detected close to the foreshock boundary. Figure[8] 
show its form deeper in the foreshock when it has spread substantially in angle evolving into half of a ring distribution 
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FIG. 8: The evolution of the ion foreshock-boundary ion beam during its convection downstream into the foreshock as has been seen 
by ISEE 2 on 04 November 1977. Left: A one-minute average pseudo-threedimensional ion velocity-space profile in the (v^;, ?;j,)-plane 
showing the spreading of the beam in angle around the bulk flow without merging into the bulk flow. Velocities are in kms~^. The scale 
on the right is count rates, and background count rates were suppressed by choosing only values above 50 s~^. Right: Contour plot of 
the partial ring distribution. The direction of the magnetic field is also shown [after l96l . The 10~^® and lO"^'^ s^ cm~^ level flux 

contours have been marked. 
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FIG. 9: ISEE 1 observation of diffuse ion beam on November 19, 1977 propagating deep inside the foreshock away from the ion 
foreshock boundary. Left: A two-minute average pseudo-threedimensional ion velocity-space profile in the (vx, Vy )-plane showing the 
undisturbed and cold (narrow) plasma inflow in negative ti^^-direction about centred and surrounded by a ring distribution of fast and 
warm (broad) diffuse ions that have been reflected from the shock but have been processed in the foreshock region when propagating 
from the foreshock boundary into the foreshock. These ions are hollow in the sense that they separate from the bulk distribution but 
have a nearly isotropic distribution function. Velocities are in kms~^. The scale on the right is count rates, and background count rates 
were suppressed by choosing only values above 50 s~^. Right: Contour plot of a similar ring a little earlier showing that the ring centre 
(star) is slightly displaced from the bulk flow (dot) on the magnetic field that connects to the shock and from the centre (cross) of the 
phase space frame [after [961 . The 10~^^ and 10~^^ s"^ cm~^ level flux contours have been marked. Note the near isotropy of the ring 

distribution. 



already. Even deeper inside the foreshock the reflected ion distribution assumes the shape of a full ring around the 
bulk distribution as is shown in Figure[9] 

The three observations depicted in Figures |4] |8] and [9] are from different times; it has, however, been checked that 
they are at distances corresponding to increasing distance from the foreshock boundary such that the assumption 
of the convectively processed beam evolution is well founded (or at least reasonable) even though it has not been 
directly proven. It is interesting to note that a gap remains between the original ion-phase space beam distribution 
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FIG. 10: Partial-density gradient of the diffuse foreshock-ion component along the shock-connected magnetic field, and density 
e-folding lengths, determined from CLUSTER observations in the Earth's bow shock foreshock [after l56l . Left: The parallel partial 
diffuse-ion density as function of distance from the shock along the magnetic field flux tube connecting the spacecraft to the shock. The 
diffuse ion density falls off exponentially with distance suggesting an ion-energy dependent diffusive process being responsible for 
transport of the ions upstream of the shock. Right: The e-folding distance of diffuse ions along the magnetic field as determined as 
function of energy from the exponential decay of the ion density. The e-folding distance increases linearly with ion energy. 



and the bulk-flow distribution, which is another indication that the evolving distribution is part of the evolution of 
the ion-foreshock beam. Of a distribution that does not evolve out of a beam one expects a less regular behaviour 
and, generally, no such well expressed gap between bulk flow and beam in velocity space. In fact, the main energetic 
ion component in the foreshock is irregular and lacks a well expressed gap. 

The discrepancy between the smooth no-gap foreshock distributions |131j and the gap-observations of [SB] has, in 
fact, been noted much earlier |112l 1115] without giving an explanation but suggesting a continuous ion source at the 
parallel shock. Below we provide further arguments for the two-source, foreshock-boundary beam and continuous 
extended shock-surface source hypothesis. Nonetheless, this conclusion must be taken with care because Figure[9] 
shows the gap progressively closing. Being sufficiently far, i.e. even farther away from the ion foreshock boundary 
then in this figure, it will not anymore be possible to distinguish between beam-evolved and genuine diffusive-ion 
distributions. Ultimately, both distributions will have merged indistinguishably. Still there is no agreement whether 
and where, i.e. at what distance from the ion foreshock boundary this merging of the two populations occurs. 

[151' , using AMPTE IRM measurements of diffuse ion densities upstream of the quasi-parallel (bow) shock found that 
the diffuse ion density decreases exponentially with shock distance. This investigation was substantially improved 
by [56] in an attempt to infer about the source of the upstream diffuse energetic ions. These authors determined 
the partial-density gradient of diffuse ions in the energy range from 10 to 32 keV as a function of distance from 
the (bow) shock. This investigation was made possible due to the availability of the CLUSTER spacecraft, a four 
identical-spacecraft mission which during this measuring period had an inter-spacecraft separation distance between 
~ 1 and ~ 1.5 Re. They used a nominal bow shock model ^98j, based on the measured upstream ffow parameters 
(basically the dynamic pressure of the flow), providing the shock-spacecraft distances along the magnetic field flux 
tube for the individual CLUSTER spacecraft and the local shock-normal angles Qbu- The average Qbu over the whole 
10 hours of observation time was 20° ± 8° proving that CLUSTER was in front of the quasi-parallel shock far away 
from the ion-foreshock boundary, and the Mach number was A4a 

Diffuse ion partial densities were determined as function of distance from the shock and in four consecutive energy 
bands every 32 s at two spacecraft. This allowed to determine the partial diffuse ion-density gradients along the 
magnetic fleld as function of energy from the density differences between two spacecraft and the differences of the 
spacecraft distances along the magnetic fleld to the shock intersection point. Perpendicular density gradients were 
neglected. The obtained parallel gradients were attributed to the average CLUSTER location. These gradients were 
then used to flnd the e-folding distance of the density variation. 

The results of this investigation are shown in Figure fTO] It is learned from this flgure that in the deep foreshock, 
i.e. at distances far away from the ion-foreshock boundary, the diffuse ion component is densest close to the shock 
with density decaying exponentially with increasing upstream distance from the shock along the magnetic fleld. This 
spatial decay of the diffuse partial ion density Ni{£, z) ~ exp[—z/L{£)] is different for particles of different energy £. 
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The e-folding distance L{£) ~ £ turns out to increase linearly with energy, i.e. low energy particles are confined to the 
shock. The higher the ion energy the deeper can the ions penetrate into the upstream plasma. The proportionality 
constant determined from these data under the special conditions of the Earth's bow shock is ~ 0.14RE/keV. This 
behaviour of the energetic foreshock ions provides indisputable evidence for the extended parallel shock-surface origin 
of the diffuse ion component. The source of the diffuse ions lies at the quasi-parallel shock. In order to be found at 
a distance upstream of the shock the ions undergo a diffusion process along the magnetic field. These ions arc thus 
completely different from the beam ions found at the ion-foreshock boundary. 

The e-folding distance for the diffuse ions is given by L{£) = k^^{£)/Vi, with spatial diffusion coefficient K||(£) = 
where is the diffusion length (parallel ion mean free path) and v the particle velocity (note that the 
diffusion coefficient has the correct dimension — m^s~^; justification of the diffusion assumption will be given in 
the chapter on particle acceleration). From balance between convective inflow and diffusion into upstream direction, 
one can write e\\{£) = 3L{£)^/&[/£ - ^/T^. The diffusion length increases as the root of the product of particle 
energy £ and upstream flow energy £i. In the solar wind the flow energy is a few keV, and a 20-keV diffuse ion will 
have a typical parallel diffusion length (or mean free path) of £|| ^ (1 — 2) Rg. This is a rather short distance, orders 
of magnitude shorter than the coUisional mean free path of an ion. Hence, the diffusion estimate suggests that strong 
wave-particle interactions can be held responsible for the scattering and acceleration of the diffuse particle component, 
which enables it to diffuse and escape upstream from the shock and populate the foreshock. The diffusion process is 
energy dependent with the most energetic ions diffusing fastest. 

These interactions should take place in the quasi-parallel shock transition because, as we have shown above, the 
diffuse upstream-ion density maximises closest to the shock. It is interesting to estimate the corresponding upstream- 
ion collision frequency Vcui — For the 20 keV-upstream ions this yields i^cui ~ 0.2 Hz. This value is comparable 

to the ion cyclotron frequency LUi.i/2TT — (0.1 — 0.3) Hz in the B ~8nT upstream to B ~30nT shock ramp magnetic 
field |56j during the time of observation. 

It seems that waves, electromagnetic and/or electrostatic, related to the ion-cyclotron frequency arc involved into 
the process of upstream ion diffusion. Since this diffusion is energy dependent, this process is not a simple pitch-angle 
diffusion as in the case of the generation of the ion beam that propagates along the foreshock boundary. The diffuse 
ion component must have experienced a substantial acceleration in this process, and this acceleration is located at or 
around the shock transition and contrasts with the ion-beam acceleration which is a scattering process followed up 
by pick-up acceleration when the upstream-propagating beam ions are subject to the effect of the main-bulk-stream 
convection-electric field in which they become accelerated in the direction perpendicular to the magnetic field to 
roughly four times the energy of the bulk flow, thereby evolving into the ring distribution that characterises their 
phase space distribution. 

B. Low-frequency upstream waves 

In the frame of the upstream bulk flow the two ion components that populate the ion foreshock carry a substantial 
amount of free energy which is subject to dissipation. Since this dissipation is coUisionless it can proceed only through 
the excitation of waves and wave turbulence through instability upstream of the quasi-parallel supercritical shock. 
On the other hand it is obvious that the presence of neither of the components can be understood without complete 
knowledge of the waves in the foreshock and their interaction with the particles. 

Since their first detection by [HI], |110j and [31], observation of shock-upstream waves has been a long-standing 
issue. Their existence was predicted by |149j . followed by hydromagnetic JB^ and kinetic |46j theories of electromagnetic 
wave excitation and propagation upstream of a coUisionless quasi-parallel supercritical shock. |163j suggested that 
they might develop into discrete wave-packets as had been inferred from observation by |110j . The ISEE 1-3 spacecraft 
allowed for a more elaborate investigation of the properties of upstream waves [cf., e.g., |35J HHJ ISHl EZl 1101 11071 11081 
11091 1131L 11321 11331 11441 among others] . More recently, CLUSTER measurements have been used to investigate the 
temporal and spatial structure of upstream waves and wave turbulence [cf., e.g., [23l |24l El EHl ESI EH UHl EHl among 
others] . We will briefly review the properties of the upstream waves in the ion-beam and diffuse ion region in view 
of the observations and mechanisms of their generation. Wave generation is coupled to particle-energy loss and to 
particle scattering both playing a substantial role in particle acceleration. It will therefore be quite natural that in 
the next chapter on particle acceleration at shocks we will return to the upstream-wave problem. 

|llj gave a comprehensive review of the various types of waves encountered upstream of quasi-parallel shocks in 
the ion foreshock. In his words, "upstream particles cause upstream waves .... Once a wave is created .., it then 
propagates, and so its continued existence relies on it remaining in a region where it is undamped. Its properties 
might even change as it propagates. Observationally, the wave propagation is superimposed on the convection of the 
plasma frame, which introduces Doppler shifts in frequency, and possible reversal of polarisation sense.... And that 
is not the end of the story, since one must take account of the feedback of the waves on the particle distribution 
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function..., and even the possibility that the ... shock injection of particles into the foreshock os modulated either by 
intrinsic processes or even by the foreshock waves themselves." 

Burgess' review was organised by the observed wave frequencies. He distinguishes between Low Frequency Waves 
(5mHz- few 100 Hz) and High Frequency Waves (> IkHz), the latter covering the electrostatic waves from ion- 
sound to electron plasma waves, as well as radiation. Radiation generation will be discussed in detail later. Here, 
we only note the almost continuous presence of waves in the ion-acoustic band which have been known since |105j 
to populate the complete foreshock region. These spectra might be composed of several different modes, ion-sound, 
electron-acoustic, Buneman modes, electron-cyclotron harmonics, and others. Their generation mechanism is not clear 
yet. They might, via a number of different instabilities, be the result of the presence of the hot foreshock-electron 
component, which also invades the ion foreshock, or they are excited by unresolved narrow electron beamlets that 
emanate from the quasi-parallel shock. They might also be excited by plasma inhomogeneities, spatial inhomogeneities 
in the electron distribution, or they are the result of nonlinear wave-wave interaction which is expected to take place in 
the foreshock. Currently these questions are difficult to answer and await further observation, simulation and theory. 
Little has changed so far since Burgess' remarks concerning high frequency waves. 

However, there has been substantial progress in the understanding of the low frequency waves and their role in 
quasi-parallel shock dynamics. The waves that are most important in shock formation propagate in the ultra-low 
frequency range < 0.1 Hz. Usually they have large (magnetic) wave amplitudes, around |b|/i? ~ 0.2 — 1.0, which 
identifies them as highly nonlinear. These large amplitude waves had already been observed by |110l to have wave 
forms from monochromatic to solitary waves, frequently with steep edges resembling shocklets and suggesting that 
the waves have experienced nonlinear steepening during their evolution and propagation. They sometimes show the 
typical signs of fluctuations that are connected to these edges and obviously propagate in the whistler mode. Thus 
these forms are indeed little shock-like structures. 

In the same frequency window, large- amplitude pulsations have been identified. These are very typical for quasi- 
parallel shocks. |126j and |125j . identifying them in the AMPTE magnetic field measurements, coined the (somewhat 
ugly as in German it means mud) name SLAMS for them, which stands for 'Short, Large Amplitude Magnetic Struc- 
tures'. We prefer to call them upstream pulsations, here. Their duration is 10-20 s; they have very large amplitudes 
|b|/i? ~ 5, indeed, but appear as a more coherent structure that is embedded in the ultra-low frequency wave tur- 
bulence. Like the ultra-low frequency waves, they propagate in upstream direction in the plasma rest frame while 
being swept toward the shock by the convective flow. Their polarisation is mixed with - possibly - left-hand polarisa- 
tion (in the plasma rest frame) slightly dominating, suggesting their ultra-low-frequency wave origin. Sometimes the 
polarisation is different on both sides of the upstream pulsation, indicating that they have evolved by some process 
which produces both kinds of polarisation, which is similar to a solitary wave. 

The propagation velocity of the upstream pulsations exhibits an interesting amplitude dependence. The upstream 
directed pulsation speed increases with amplitude, which also is a solitary wave-like property. Moreover, they grow 
when approaching the shock and entering the increasing density gradient of diffuse ions, and they play an important 
role in shock reformation. |125j suspect that these pulsations are the 'building blocks' [an expression used bv ll26| of 
quasi-parallel shocks. 

Another interesting property is that upstream pulsations contain thermal plasma with properties of the upstream 
flow, while being surrounded by the hot foreshock plasma. One would therefore believe that their source region 
is located at the ion-foreshock boundary. Being created there by an ion-ion beam plasma interaction they might 
grow nonlinearly until reaching quasi-equilibrium like solitary structures, having captured the upstream plasma, and 
afterwards being convected toward the shock into the heart of the ion foreshock. We will return to these interesting 
structures when discussing simulations below. 

Ion-beam waves 

Each of the two different upstream-ion populations is responsible for the excitation of its own instabilities. In this 
section we deal only with those waves which are excited by the foreshock-boundary ion beam. From Figure[6]we obtain 
that the situation is that of an ion-ion beam (when for the moment neglecting the gyrophase-bunched diffuse-ion 
component). In the frame of the upstream flow the reflected foreshock-boundary ion beam propagates upstream along 
the magnetic field at parallel speed ~ —(2 — 3)Vijj, where Vi\\ ~ 600kms~^. The upstream ion temperature (in 
energy units) for this case was ~ (1 — 2) keV, yielding roughly a thermal ion velocity of Wj ~ 200kms~^. Moreover, 
the beam can be taken as warm with thermal speed vt^th ^ Hence, Wf,j| > with magnetised background and 
beam ions. On the other hand, the electrons are hot with ~ 100 eV. 

a. The expected wave modes. Since the large parallel speed of the beam corresponds to a large parallel temper- 
ature anisotropy, it is clear that the beam can excite long-wavelength negative-helicity Alfven waves via the firehose 
instability. In addition, because the beam is mildly warm, it can excite the resonant left-hand ion-ion beam instability, 
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FIG. 11: Reduced parallel distribution functions along the ion-foreshock boundary showing the cold main flow ion distribution _Fi(i;||), 

the hot main flow electron distribution _Fe(f||) shifted to the left into the direction of the foreshock ion beam for keeping zero-current 
conditions, and the ion-foreshock ion beam distribution Fii,{v^^). This configuration is unstable with respect to ion-ion beam instabilities 

and the ion-beam driven ion-acoustic instability. 

which is possible for vi,\\ > Va ^ 100 km s^^. And under conditions, when the beam thermal speed can be considered 
to be small, it excites the right-hand resonant ion-ion beam mode. Both waves propagate with the beam upstream 
along the magnetic field on the background of the upstream flow. They are not as fast as the beam, however, and are 
thus subject to downstream convection with the flow towards the shock. The firehose mode, at the contrary, moves 
against the beam and thus by itself moves downstream in the direction of the shock position, when excited. 

In all three cases the foreshock-boundary beam will lead to the excitation of low-frequency Alfven and ion cyclotron 
waves which in the shock frame approach the shock while having their source on the foreshock-boundary field line. 
During this shock-directed convection and/or propagation they populate the foreshock region with low frequency 
electromagnetic fluctuations, which might further interact with the diffuse foreshock-ion component. On the other 
hand, since the phase and group velocities w/fcj| ~ {dLo/dk\\ ) ~ Va <C Vi,\\ of these waves are of the order of the Alfven 
velocity and are, thus, much less than the beam and flow velocities, downstream convection will quickly remove them 
from the foreshock-boundary source region. Hence, their further evolution in the foreshock is determined by the 
competition between nonlinear wave steepening and interaction with the diffuse foreshock-ion component. 

Including the electrons (while so far neglecting the electron foreshock component) leaves us with an ion-acoustic 
unstable phase-space configuration. In the upstream ion-plasma frame a relatively dense foreshock-boundary ion- 
beam is propagating upstream on a cold ion-hot electron plasma. In order to keep the plasma current-free the 
electron component is slightly retarded creating conditions under that ion-acoustic waves can be excited. On the 
other hand, the configuration is not able to excite neither the Buneman-two stream instability nor - because the 
beam propagates solely parallel to the magnetic field - the modified-two stream instability. 

This is shown in Figure [TT| The instability is excited by the velocity difference between the ion foreshock-boundary 
ion beam Fn,{v\\) and the slightly shifted to the left hot (Tg ^ 100 eV) electron distribution because the velocity 

difference AV between ion beam and electron component > AV ~ lO'^kms"^ > Cia ''^ lOOkms"^ exceeds the ion 
acoustic speed while being less than the electron thermal velocity v^- These ion acoustic waves occupy a relatively 
broad spectrum with downstream parallel phase velocities < Vib and become swept towards the shock. Since they 
quickly leave the foreshock boundary and since the beam is not hot, ion-Landau damping in the source region plays 
no role. However, when entering the diffuse-ion foreshock region these waves encounter the hot diffuse ion component, 
will interact with it, and will thereby experience Landau damping. 

b. Observations. Observationally, it is difficult to distinguish between foreshock-boundary waves and diffuse- 
ion generated waves. Two types of waves that can be related to the presence of the upstream beams have been 
reported in the vicinity of the bow shock. [31 found monochromatic large-amplitude discrete wave packets in the 
frequency range lojl-K ~ 0.4 Hz which [47 ^ could show to be right-hand polarised whistlers at w ~ lOw^i propagating 
on the plasma rest frame and being unable to escape far upstream. The other class of waves is of smaller amplitude 
\b\IB ~ 0.1 and frequency of the order of w/27r ^ IHz. These waves form trains which are directly tied to the 
upstream foreshock-boundary ion beams |50| . In fact, as |50j demonstrated they occur only in the presence of the 
observed reflected 2 — 5 keV foreshock-boundary ion beams [H] which |57] have shown to evolve from a beam into a 
gyrating particle component that in the deep foreshock is superimposed on the diffuse ion component before it merges 
into it. However, as [? ] have shown from measuring the electron distribution function, it is not the ion beam who 
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FIG. 12: Left: The combined power spectral density of magnetic fluctuations excited by the ion foreshock-boundary ion beam with 
central frequency around ^ 1 Hz [lower part, after 50 and at < 0.1 Hz [left upper part, after 167] as measured by the ISEE spacecraft. 
This spectrum is obtained in the ISEE spacecraft frame. The low frequency spectrum is shown at three different times corresponding to 
(from below upward) increasing distance from inside the foreshock to the foreshock boundary. It is about five orders of magnitude more 

intense than the high-frequency waves. It is worth noting that the high frequency spectrum of |50| toward lower frequencies in the 
overlap with the 1671 spectrum does not show any indication of the lower frequency peak. As indicated by the thin straight line, it would 
smoothly continue into the 1930 UT branch at 0.01 Hz. This lack of the low frequency peak is probably accidental and due to the 
particular conditions at the time of measurement. Right: The wave dispersion relation for the high frequency ~ IHz whistler- mode 
waves determined from a tentative estimate of the relavant wave numbers and transformed into the plasma rest frame. It is seen that 
these waves in the observation range have about linear dispersion and frequencies uj ~ (20 — 100) a; ci, far above the ion-cyclotron 

frequency u> ci [after |50) . 



excites these whistlers. The probabihty that they are driven by the particular electron distribution in the foreshock 
is rather higher. 

This is a very favourable case as these waves carry information about the generation mechanism. The waves 
propagated obliquely (< 60°) with respect to the magnetic field, at an average angle of ~ 45°. The observed spectrum 



and the dispersion relation determined from the measurements are shown in the high-frequency part of Figure 12 



This figure has been combined from the high-frequency observations of [50| and the low-frequency observations of [67^ , 
both obtained from the ISEE spacecraft. Tentative wave numbers have been determined for the high-frequenc waves 
from measuring the time delay of the wave front arrivals at the two spacecraft ISEE 1 and ISEE 2 spacecraft yielding 
(surprisingly) short wavelengths < 100 km. Knowing the wave number, the frequency has been back-Doppler shifted 
to = uj ISEE ~ k • Vi into the plasma frame yielding unusually high frequencies lu ~ (20 — 100) w^i. This procedure 
determines the dispersion relation (naturally with large errors as indicated by the bars), which is found to be about 
linear in the narrow high frequency range of the waves in the rest frame of the plasma. The waves that have been 
left-hand polarised in the spacecraft frame turn out to become right-hand polarised in the plasma rest frame. 

From these properties it was initially concluded that the waves propagate in the whistler branch and have been 
excited by the cool ion-ion beam instability, which would be consistent with our initial discussion. This conclusion is, 
however, questionable. The high rest-frame frequencies are not in agreement with model calculation [311 1131] using 
the realistic observed beam properties. These yield wave frequencies of the order of lo 0.1 Ud for both the ion-ion 
beam and firehose modes. In fact, [27 report CLUSTER observations of similar waves but with much lower frequency 
Lu/uJci ~ 0.1 and wavelength the order of ~ 1 Re ~ 6000 km, which is in excellent agreement with the theoretical 
predictions. These waves are cold ion beam-excited fast-kinetic (magnetosonic) whistlers, in the terminology of [38] 
and |63] . The high-frequency waves ~ 1 Hz whistlers are instead most probably driven by the particular electron 
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FIG. 13: Left: The excellent correlation between the density and magnetic field variations in the ion-beam excited foreshock-boundary 
low-frequency waves as measured by CLUSTER [after [24]. The variations are practically in phase thus identifying the fluctuations as fast 

magnetosonic. 



distribution in the foreshock [Sj rather than by the cold-ion beam instabihty. This claim is also supported by the fact 
that deeper in the ion foreshock, where the electron distribution becomes more isotropic, these waves do not occur 
separate from smaller magnetic field structures (large-amplitude magnetic pulsations of the SLAMS type or shocklets) . 

In quite good agreement with the CLUSTER measurements are the ISEE foreshock waves analysed by |H7j who found 
the spectral peak at frequencies uj ~ 10~^ Hz with the spectrum broadening with increasing distance from the ion- 
foreshock boundary. This is shown in the upper left part of the spectrum in Figure 12 where the measurements of [SD] 
and [67j have been combined. Even though the conditions on the two observation times were not identical, one sees 
that the high-frequency whistlers identified by [50j occur on the approximate high-frequency extension of the [57] 1930- 
2030 UT spectrum which is closest to the ion foreshock boundary. Note, however, that the low-frequency part of the 
|50j-spectrum (not shown here other but indicated by a thin straight line in the figure) was flatter and in the (relatively 
short) overlap region with the 67J-spectrum did not show any indication of the magnetosonic foreshock-boundary 
waves, which must be due to the particular conditions prevailing during this observation period. Thus the [50 j waves 
are a (high frequency) wave species that is different from these more common low- frequency/long- wavelength fast- 
magnetosonic waves. These latter waves are also of much higher (~5 orders of magnitude) spectral density as is 
seen from Figure 12 and has been confirmed by the later CLUSTER measurements. Moreover, the magnetosonic waves 
are left-hand polarised in the plasma frame, and their compressive fast-magnetosonic character is proved from the 
in-phase variation of the density and magnetic field fiuctuations shown in Figure [T3| 

The presence of these low-frequency ion-beam generated waves implies that the ion beams interact with the fiuc- 
tuating electromagnetic field. In this process they become scattered and diffuse in phase-space. This is the reason 
for the spreading of the ion beam in phase space and the final merging into the diffuse background distribution. [57] 
have followed this evolution of the beam as we have described above. It is responsible also for the gradual spreading 
of the spectrum in the upper left part of Figure[T2] A spectrum similar to that given there can be found in the paper 
of [ST]. [2], from CLUSTER observations, inferred the nature of these ultra- low frequency (ULF) waves and showed 
that their correlation lengths along the wave vector direction k is of the order of 1 — 3 Re, while it can be a factor of 
three larger in the direction perpendicular to k, rendering these waves oblate though nearly planar. 

During times of hot beams |23j observed upstream propagating left-hand low-frequency waves which have been 
excited on the Alfven-ion-cyclotron branch of the kinetic dispersion relation. These waves are kinetic Alfven waves 
which have been excited by the hot upstream propagating beam similar to the one shown on the right in Figure|6] 
In the spacecraft frame these waves because they are swept downstream by the fiow, appear as right-handed waves. 
[23J report that with onset of the waves the beam gets more diffuse. This can either be interpreted as the reaction 
of the waves on the beam or that the spacecraft enters the region where the initially cold beam enters the foreshock, 
spreads in velocity and after becoming hot enough generates the observed waves. 
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FIG. 14: Left: The dispersion relation of low- frequency waves as measured by CLUSTER in the foreshock away from the foreshock 
boundary [after l85l . determined from the "wave-telescope" analysis method. The upper panel shows the dispersion relation tij(fc) 
consisting of several parts and distinguished by the sense of polarisation as right-hand, left-hand and linear. The different branches 
cluster. There are forward (positive frequency parallel to the magnetic field) and backward (negative frequency antiparallel) branches. 
The lower panel shows the wave propagation angle for the shorter wavelengths being mostly oblique ear 9 (30 — 25)°, for longer 
wavelengths being more perpendicular to the magnetic field. Right: The five low-frequency modes in a cold beam plasma system of 
similar conditions as the foreshock plasma and for the propagation angle 8 = 24° showing similarity to the observed dispersion relation 

m- 



The question arises whether there are simulation studies available of the evolution of the reflected ion beam and the 
generation of upstream waves. This question cannot be definitely answered at present. It is clear that a simulation of 
this kind should reconstruct the region of the foreshock boundary which requires that a curved shock surface must be 
assumed from the beginning. In other word, this problem can be investigated only in a two-dimensional simulation. 
Two-dimensional simulations are available in hybrid form but have have not been applied to curved shocks. Hence, 
this problem remains an open simulation problem. On the other hand, in the following we will extensively discuss 
the role diffuse particles play in shock formation and upstream wave generation for planar shocks in one and two 
dimensions and hybrid as well as PIC simulations. 



Diffuse ion waves 



Long-period ~ 30 s waves are the rule in the foreshock. They occur together with the diffuse foreshock-ion component 
[cf., e.g.. [SUl FllSj which is located deeper inside the foreshock. Because of this reason, any waves that are excited by 
the diffuse ion component are restricted to the interior of the foreshock. One even has defined some fuzzy boundary 
of these waves called ULF-wave boundary jlSJ 1109] , which is even more inclined against the upstream magnetic field 
than the ion-foreshock boundary and outside of which the ULF-wave activity should be weak. The reason is that ULF 
waves, if propagating upstream, can move at most with fast-magnetosonic velocity which for a supercritical shock is 
smaller than the stream and also less than the reflected ion beam velocities. The advection by the flow will blow them 
downstream towards the shock and confine them to a region relatively close to the shock [5T] bounded to upstream 
by the ULF-boundary. 

|85] determined the dispersion relation of the low-frequency waves in the foreshock from the three-dimensional 
observations of CLUSTER. Their result is shown in Figure[l4] The left part of the figure shows the dispersion relation 
in the plasma rest frame and the angle of wave propagation with respect to the local average upstream magnetic field. 
The scatter is quite large. Nevertheless it is surprising to find a well expressed nearly linear part on the dispersion 
relation even though the dispersion relation is composed from the contributions of several distinct modes with different 
polarisation. The representation is linear, which, on the fc-axis, emphasises the short wavelengths (larger k values). 
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FIG. 15: Left: Short wavelength spectra for the magnetic fluctuations measured by CLUSTER in the foreshock away from the 
foreshock boundary. These spectra are obtained by apphcation of the "wave-telescope" analysis method. The perpendicular fluctuations 
are more than one of magnitude more intense than the parallel fluctuations while the spectra decay about according to the Kolmogorov 
— |-law of stationary inertial turbulence. Indication of a cut-off is seen fat wave numbers larger than the inertial wave number ki„ where 
strong dissipation sets on. From the maximum at fc ~ 10~^ one concludes that in this wavelength range the spectral energy is injected 
by instability of fast and Alfvenic magnetosonic waves which cascade nonlinearly forward towards shorter wavelengths. Right: The 
probability distributions of the parallel and perpendicular magnetic fluctuations. The dashed distributions are log-normal of same 
maximum. These distributions exhibit extended tails and thus indicate non-stationary and probably intermittent not fully developed 

magnetic turbulence I90| . 



From linear theory one expects that the long wavelengths (small k) should be more pronounced. Unfortunately, this 
region could not be resolve sufficiently at the available CLUSTER separation distances at that time |89]. These low 
frequency waves, in the plasma frame, seem to propagate close to perpendicular to the magnetic field and, therefore, 
are probably in the fast magnetosonic mode. Moreover, it seems that a straight beam dispersion relation contributes 
to the perpendicular propagating waves. It also seems that in the low-frequency (magnetosonic) waves have negative 
frequency. Linear dispersion theory for the parameters during the observation time yields the curves on the right in 
the figure showing the coupling between the presumable beam mode and the four plasma modes at wave numbers 
kVA/uJci < 0.3. The experimentally determined dispersion relation resembles this clustering of couplings at small 
k and very low frequencies uj. However, application of linear theory is dangerous if not questionable because the 
foreshock plasma is highly disturbed, the low frequency waves have rather large amplitudes and can barely be related 
linearly, and the plasma is very inhomogeneous exhibiting steep gradients in density and field. For an Alfven speed 
Va — 30kms~^ and an ion-cyclotron frequency a;/27r ~ IHz, wavelengths around the linear wave coupling are of the 
order of fc/27r ~ 100 km and should thus be affected by the plasma inhomogeneities (note that decreasing the reference 
cyclotron frequency decreases the wavelength even further). Nevertheless, at the very low frequencies and very long 
waves the waves are probably in the fast magnetosonic wave band. 

More mysterious are the higher frequency-short wavelength waves. From comparison with the linear dispersion 
relation they seem to fit on the right-handed whistler branch (high frequency magnetosonic whistler R"*"). However, 
inspection of the dispersion relation indicated that all kinds of polarisation are scattered along the dispersion curve. 
The curve itself is very irregular even though the propagation angle of the waves seems to be about constant at weakly 
oblique angles. This wave composition suggests that we are not dealing here with one single wave mode but rather 
with the short wavelength part of a turbulent spectrum which has generated all kinds of short wavelength fluctuations 
with different almost randomly distributed polarisation in a forward cascading process from long to short wavelengths. 
This idea has been elaborated in more detail in [90, who found that the shorter wavelength spectrum is indeed about 
featureless and power-law, close to a — |-Kolmogorov-spectrum of stationary turbulence with most of the power in 
the perpendicular (non-compressive) magnetic component providing another argument for Alfvenic and magnetosonic 
turbulence. This is shown on the left in Figure [T5l The figure also indicates the respective gyro- and inertial wave 
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numbers, kd = Vi/ujci and kin = c/(^pi- There is indication that the spectra change slope at around these number 
due to onset of ion viscosity and inertia. The large maximum on the perpendicular power curve near k ~ 10~^km~^ 
corresponds to a wavelength of A ~ 6000km reported earlier [e.g., US] and is at the right position for energy injection 
by instability. 

However, the conclusion that the deep foreshock is subject to fully developed fast magnetosonic/Alfvenic turbulence 
is not fully justified. This becomes obvious from considering the probability distribution of the field fluctuations as 



given on the right in Figure 15 Would the fluctuations be normally distributed then their distribution functions 
would have the shape of the dashed curves. Instead the distributions are skewed into tails. In the case of the 
parallel component the distribution lacks symmetry. Such tails may indicate that the turbulence is non-stationary, 
intermittent or inhomogencous. All three cases might hold in the foreshock, in particular as the foreshock is a 
rather limited spatial region which is bounded from two sides and is subject to plasma injections and plasma losses. 
Nevertheless, the spectra determined contain signs of strongly nonlinear and turbulent interactions, and the waves, 
particularly the low-frequency waves, are of large amplitude and interact with the plasma ion component as also with 
other waves. These waves affect the upstream backstreaming ion component scattering, heating, and accelerating it. 
On the other hand the upstream ion component is responsible for the existence of the waves as it is the ultimate 
energy source of the waves. And, to close the cycle, the waves cannot escape upstream very far from the foreshock 
since in a supercritical shock the Mach number of the flow is higher than any Mach number based on the wave 
speed and is thus high enough for the flow to advect the wave spectrum towards the shock. As long as the waves 
do not completely dissipate their energy during this convection, the electromagnetic wave energy accumulates at the 
location where they arrive at the shock. The consequence of this accumulation is that the waves affect the shock, 
cause instability of the shock surface, and reorganise the shock front. This makes quasi-parallel shocks non-stationary 
and subject to some kind of irregular reformation. One may thus expect that a supercritical quasi-parallel shock does 
not represent a solid shock surface over a large area. It consists of a more or less dense patchwork of areas which 
together form a shock but which also appear and disappear in an irregular manner, come and go, and organise the 
shock in a certain volume where the entropy increases but where a multitude of very difficult to handle processes 
takes place that can be investigated only experimentally or with the help of properly designed numerical simulations. 
These are, however, more difficult to design than in the quasi-perpendicular case, because of the greater variability 
of the conditions at a quasi-parallel shock. Since waves and particles are tied to each other, simulations cannot be 
discussed separately for waves and for particles. We therefore delay the discussion of the simulation results to the 
later section on quasi-parallel shock reformation. There it will become obvious what the upstream waves are good for 
and what their role is in the quasi-parallel shock process. 

An important question that has not been addressed anywhere in the investigation of the ion-foreshock dynamics 
concerns the role of the foreshock electron component which, in the large ion-foreshock domain, consists of two 
populations, the ~ 100 eV upstream electron population which belongs to the quasi-neutral upstream ffow, and the 
hot ~ several kev-electron population which is the product of the shock reffected electrons we are going to discuss 
in the next paragraphs. This component is isotropic but irregular and might contribute to ion-foreshock instabilities 
thereby affecting the foreshock turbulence and shock formation. Since these questions are difficult to be treated 
theoretically one has to wait until numerical simulations will become capable of including them into full particle codes 
in a similar way as has been done [cf. [75l|76] for the electron-ion populations in the feet of quasi-perpendicular shocks. 



C. Electron foreshock 



Even though the ion-foreshock occupies the larger part of the foreshock and, which is more important, plays the 
decisive role in the dynamics of the quasi-parallel supercritical shock, it would be an ignorant attitude not to mention 
the part that aside of the main flow is populated solely by the reflected higher energy electron component, i.e. the 
electron-foreshock region. Schematically its extent and relation to the ion-foreshock is shown in Figures[3| and \TE\ 

The electron foreshock has been identified already by |114| in the OGO 5 satellite observations. fll4| noted the 
occurrence of enhanced electron fluxes near the bow shock, which were related to observations of high-frequency 
electric-field spikes at 30 kHz. They concluded that these spikes resulted from Langmuir waves that had been excited 
by electron beams arriving from the bow shock along shock-connected magnetic field lines. (A detailed overview of 
the early observations can be read in |SH].) It could be confirmed that these electron fluxes were magnetic-field aligned 
and of higher energy than the bulk electrons in the main flow. This was interpreted as electron beams emitted into 
upstream direction from the shock, even though no mechanism was known that could provide the required shock- 
electron acceleration - in fact, a de Hoffman- Teller-frame shock-reffection mechanism had already been proposed by 
|138j and was reinvented and elaborated on much later by |164j . 
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FIG. 16: Synoptical view of the two foreshocks: the electron and ion foreshocks, respectively. 



Electron beams 



In analogy to the ion-foreshock boundary the electron-foreshock boundary (which is the ultimate upstream boundary 
of the foreshock) carries narrow bursts of electron beams, which escape into upstream direction from the shock along 
the magnetic field and which are slightly displaced by the convective fiow into downstream direction. 

An example of a narrow energy spectrum of such a beam measured by the IMP-8 satellite is shown on the left in 
Figure fT?! The spacecraft was in the electron foreshock boundary for a relatively short time only. Being near apogee, 
the spacecraft was about standing. Thus the short contact time with the beam most probably indicates that the 
electron-foreshock boundary is fairly narrow. The electron beam also occupies only a small volume in velocity space 
corresponding to a narrow bump on the distribution [1]. Only one energetic electron-gyroradius deeper in the electron 
foreshock the electron beam becomes depleted and just forms an energetic tail of hot halo electrons on the electron 
distribution function of the flow. 

Figure 18 shows electron phase-space observations from the electron instruments on ISEE 1 & 2 during crossings 



of the electron-foreshock boundary. The left part of the figure [taken from |36] is a much higher time and velocity 
(respectively energy) resolution plot than that in the former figure. It is nicely seen how the electron distribution 
evolved from a field-aligned (nearly) Maxwellian distribution prior to contact with the foreshock boundary, through a 
field-aligned-beam-like distribution at contact, into a distribution with an energetic tail along the magnetic field, where 
the electron beam has been completely washed out. This transition takes place within a time interval of 15s, just 
allowing to obtain the three measurements. The pseudo-threedimensional plot of the electron distribution function on 
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FIG. 17: Left: The electron energy distribution as seen by IMP-8 at three times when crossing the foreshock boundary. Right: The 
high-frequency wave spectrum of the electromagnetic radiation as measured some distance away from the shock near the foreshock 
boundary. Two types of emissions are visible, the so-calle 'harmonic' emission / = 2/pe at twice the electron plasma frequency 
/pe = a;pe/27r, and the weaker so-called 'fundamental' emission at /pe. [after [44]. 

the right [taken from I was obtained on another day when the spacecraft passed the electron-foreshock boundary and 
a series of distributions were recorded from upstream, across the foreshock boundary, and into the electron foreshock. 
Here only the foreshock-boundary distribution is shown when the narrow electron beam has evolved flowing into 
upstream anti-shockward direction along the magnetic field. Note the narrow angular extension of the beam, its 
weakness, and its comparably large velocity which allows to identify it just outside the bulk of the hot ^ 100 eV 
background electrons which are much hotter than the beam (indicated by the small beam width in v\\). 

The conclusion that can be drawn from these observations, which in the follow-up time had also been confirmed by 
measurements of other spacecraft like AMPTE IRM and AMPTE UKS, is that similar to the ion foreshock, the electron 
foreshock consists of a very narrow region along about the shock-tangential upstream-magnetic field line (or fiux tube) 
which is populated by a cold, weak, and fast electron beam, and a broad electron region where the electron distribution 
exhibits tails and in addition contains a hot, fairly isotropic electron component. The electron beam (like the ion 
beam at the ion-foreshock boundary) has its origin near the shock-tangential magnetic field line which is connected to 
the quasi-perpendicular part of the shock. The mechanism of its generation has not yet been unambiguously clarified. 
Presumably if consists of a combination of the de Hoffman- Teller frame mechanism of |138j and the bending of the 
shock surface plus some kind of stochastic acceleration in the electric fields which evolve in the shock transition, ramp 
and overshoot. Such fields have indeed been reported recently on scales below the electron skin depth Ae = c/uipe to 
be very large, of the order of < lOOmVm"^ parallel to the magnetic field and < lOOmVm^^ perpendicular to the 
magnetic field [5 . The presence of the hot electron component in the foreshock, on the other hand, requires another 
mechanism that heats and accelerates electrons at a quasi-parallel shock over a large area of its surface. 

Langmuir waves 

So-called 'gentle electron beams', which are just those beams one observes in the electron-foreshock boundary - fast, 
but not too fast (V^few times Vi), parallel (Vf, = Vf,||,V(,j_ = 0), weak (A^f, ^ N^), cool (wefc < f e ^ H), are know to 
be the drivers of Langmuir waves with dispersion relation 

Lj\k) - UjI^ + ivie, V,,s > VSVe (1) 

via the resonant kinetic gentle-electron beam instability. The Langmuir resonance conditions is w — k • Vres — 0, where 
Vrcs ~ is the resonant electron velocity. Combining it with the dispersion relation, it is easy to see that in order 
for the wave number k to be real the resonant velocity must satisfy the condition on the right in the above equation. 
Below we will discuss the properties of these waves in more detail in connection with observations. 
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REDUCED PARALLEL ELECTRON DISTRIBUTION 




FIG. 18: Left: Three successive reconstructions of the reduced parallel electron distribution function during crossing of the electron 
foreshock boundary. The distributions are only seconds apart, at 1650:37 UT before touching the foreshock boundary, at 1650:43 UT 
just crossing it, and at 1650:52 UT being behind it in the electron foreshock. The beam is visible only during the short crossing time. 
Afterwards the electron distribution shown a heated non-symmetric tail, indicating that the spacecraft sees the depleted beam in the hot 
electron- foreshock plasma [after [36]. Right: The full two-dimensional electron distribution function during another crossing of the 
electron-foreshock boundary at the moment when the beam is visible as a narrow enhancement of the distribution in the negative 
parallel velocities (upstream along the magnetic field). An instant later, when the spacecraft entered into the electron foreshock the 

beam had disappeared [after [l]. 



c. Stability of the electron beam. However, there are a number of severe limitations to this instabihty which 
theoretically would inhibit its application to the electron-foreshock boundary beam. The most stringent of these 
limitations is that it is believed that quasi-linear saturation of the gentle-beam instability should quickly, i.e. within a 
time as short as about Atujpe ^ 10, deplete the beam and transforming it into a plateau on the distribution function. 
For an upstream plasma density iVg ~ 5 x 10^ m~'^ the plasma frequency is cjpe/27r ~ 20 kHz. Hence, the beam should 
become depleted within At ~ 0.08s. A beam of velocity Vb ~ lO^'kms^^ (which satisfies the above condition on the 
resonant velocity in an inflow of Tg 100 eV) would thus propagate just 800 km upstream of the shock along the 
magnetic field before being depleted, with the implication being that the beam could not escape from its shock-source 
region. This distance is much less than the distances from the shock at which both the electron-foreshock boundary 
beams and the beam-excited plasma wave spectra have been observed. Either quasilinear saturation of the Langmuir 
instability does not work, or the observed plasma waves are no Langmuir waves. The latter is probably not true. So 
the question arises of how the beam can avoid becoming depleted due to quasilinear effects. 

There are several possibilities of inhibiting beam depletion. It has for long time been believed that collapse of 
Langmuir waves resulting from modulation instability (also known as oscillating two- stream instability) could shift the 
beam-excited Langmuir waves out of resonance and remove them from reaction onto the beam. This is a very elegant 
mechanism which is based on the ponderomotive force the Langmuir waves exert on the plasma. The ponderomotive 
force is the gradient of the wave pressure which is strongest right in the place where the Langmuir waves are most 
intense, i.e. in the beam region. For the inert ions of the plasma background the intensity gradient V|ep of the 
fast large-amplitude oscillating Langmuir electric field presents a pressure force which pushes them out of place thus 
decreasing the plasma density locally. Since the electrons react immediately and follow the expelled ions in order to 
maintain quasi-neutrality, the wave-ponderomotive-pressure force drills holes into the plasma which are filled with 
Langmuir waves. The density variation produced 

^'"'I'^ia 

is nothing else but an ion-acoustic wave. The Langmuir waves, when becoming intense enough, can drive ion-acoustic 
waves of speed Cia, which exist only due to the presence of the Langmuir wave pressure gradient. The broad Langmuir 
waves in this way organise into a large number of density cavities, and the waves are removed from resonance with the 
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beam electrons. Enhancement of the outside pressure then forces the holes to shrink. This shrinkage of the hole size 
shortens the Langmuir wavelength, increases the wave number, reduces the phase velocity and, in this way, shifts the 
waves into the main electron distribution where they can undergo Landau damping and dissipation. Unfortunately, 
this kind of modulation instability and collapse of Langmuir waves could never be experimentally approved at shocks. 
It is therefore quite unlikely (though not impossible) that it occurs and stabilises the beam. 

Another possibility of stabilising the beam by pushing the Langmuir waves out of resonance is scattering them on 
thermal ions. This process reads i -\- % ^ I' -\- i' . For purely elastic scattering only the momenta /ik^ of the Langmuir 
wave i and Pi of the ion i are changed according to the momentum conservation equation 

m^g + p» ^ Sk^ + p ■ (3) 

where the primed quantities are after collision of wave and ion. Loss of momentum decreases the wave number kg 
and shifts the wave phase velocity to higher values out of resonance with the beam; gain of momentum increases 
ki and brings the wave phase velocity down into the bulk of the electron distribution where it is again dissipated. 
The efficiency of this process is, however, low and presumably insufficient for complete stabilisation of the beam. In 
addition it strongly depends on the available particle number; it is thus most probable to work in the immediate 
vicinity of the shock only. It does, however, produce a new Langmuir wave I' that plays a role in the generation 
of radio-radiation from a shock. Because of this capacity it nevertheless remains to be of interest in the physics of 
coUisionless shocks. 

[81] investigated this nonlinear wave scattering off ions including the polarisation cloud which accompanies each ion. 
They found that of the ion distribution the ions with velocity v ~ Ivi are the most effective scatterers of Langmuir 
waves. The inclusion of the polarisation cloud increases the scattering efficiency by a factor of three which does not 
change the above conclusion. However, the wave number of the scattered wave, its direction of propagation and other 
wave properties change substantially. Nevertheless, the effects are not strong enough to alone explain the long survival 
of the electron beam over distances very far from the shock. 

A review of all the available theories has been given by [83j who performed simulations of a spatially bounded 
beam. It seems that many effects act together in order to allow the beam to survive. Weak scattering off ions removes 
some waves from the beam. The remaining waves are only partially reabsorbed because the fast beam front runs 
them out such that they stay in the trail of the beam only. In this way only the trail part of the beam is depleted 
and retarded, which affects only the low velocity component of the beam distribution, assuming that the beam has a 
small but sufficiently large natural beam spread Augb <C We, as is suggested by the observations shown, for instance, 
in FigurejlSj If injection is continuous, this implies that the beam starts pulsating, which may be a reason for the 
observed strong Langmuir wave variations along the electron foreshock-boundary field line. Finally, the Langmuir 
waves, riding on the beam, are generally slower than the beam; since v^^^ > Sv^, their group velocity, which carries 
the energy that can be absorbed, is at most V£g < v^y/i. So they are more vulnerable to the downstream convection 
than the beam electrons. The convective VL-niotion of the upstream flow thus sweeps them readily out of the beam 
region at the electron foreshock-boundary field line. All these arguments taken together explain why the beam can 
survive over long distances and, in addition, why the region across the electron foreshock boundary where the beam 
is detected is as narrow as found in the observations. 



d. Electron foreshock-boundary waves. Figure 19 taken from |18j . shows a beautiful double passage of the 
CLUSTER-spacecraft quartet through the Earth's bow shock on December 22, 2000 in the light of the CLUSTER plasma 
wave spectrum recorded by the WHISPER plasma wave instrument aboard the four CLUSTER spacecraft. Figure[20| 
shows an example of these high-resolution measurements 

CLUSTER came from the magnetosphere, entered the downstream region of the shock, passed the shock to upstream, 
re-entered the shock, and escaped downstream again. The relevant signatures of each region are seen in all four 
spacecraft. The uppermost light blue trace is the electron plasma frequency fpe = tUpe/^ir which maps the local 
plasma density. The shock appears as the intense broadband emission with maximum intensity in the low-frequency 
waves. In the density it is mapped as a steep drop in fpe to the low upstream density values. 

Close to the shock a number of intensifications in /pe can be recognised. These occur at the time of contact with 
the electron foreshock boundary when intense plasma waves are excited in the Langmuir mode. It is interesting to 
note that the intensification occurs in spots and is not necessarily narrow-band. It may exceed /pe, or it may also drop 
below it: the plasma frequency may have 'hair' or 'beards'. Moreover, sometimes intense lower-frequency emissions 
are detected half way between the plasma frequency and the low frequency emissions in the foreshock. The most 
intense of these emissions belong to the region closer to the shock where the density gradient has not yet settled to 
the upstream values, but weaker emissions of the same kind occur farther away in the upstream low density domain. 

A detailed analysis of all these waves detected by CLUSTER has not yet been undertaken and thus is not yet available. 
It is, however, highly suggestive that the observed intense spots in the plasma frequency, with the 'hair' parts exceeding 
and the 'beard' parts hanging down from the local /pe, are related to magnetic field-aligned electron beams emanating 
from the shock and propagating upstream along the electron foreshock-boundary magnetic field as these are the same 
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FIG. 19: A CLUSTER spacecraft passage across the Earth's bow shock region on December 22, 2000. The spacecraft are coming from 

the magnetosphere, pass across the downstream magnetosheath region with its enhanced density, cross the shock into the upstream 
region and back again and remain in the downstream magnetosheath. The hole region is seen in the light of the broadband plasma wave 

electric field spectrum from about kHz up to 80 kHz. The figure shows the four panels of all four CLUSTER spacecraft (with their 
somewhat childish names Rumba, Salsa, Samba, Tango). The uppermost light blue emissions in each panel belong to the electron plasma 

frequency fpg, = u>pg/2iT which on most of the path is thermal noise mapping the local plasma density. The passage of the shock is 
signalled by the appearance of intense broadband waves starting at low frequencies and being correlated with a fairly steep drop in the 

plasma frequency respectively the plasma density. Outside the shock in the low density upstream region the strong spots of 
intensification of the plasma frequency indicate contact with the electron foreshock boundary beam and excitation of Langmuir waves, 
the intensification is highly structured, sometimes stretching above, sometime hanging down below fpe- Also seen are intermediate 
frequency emissions near the shock. The lowest panel shown the estimated from fpe plasma density for all four spacecraft [after [Ts] . 



signatures as those observed with the ISEE and AMPTE spacecraft, though with much better frequency resolution and 
higher sensitivity here. For instance, the latter two spacecraft were unable to resolve the plasma frequency which 
in the CLUSTER data is nicely distinguishable even though it represents just thermal noise. It is most interesting 
that these contacts with the electron foreshock boundary occur quite irregularly throughout the hole period when 
CLUSTER was upstream of the shock and that they do not seem to exhibit a one-to-one correlation between the four 
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FIG. 20: An example of ISEE 1 high frequency wave observations during contact and passage of the electron foreshock boundary. The 
thin faint spotty line around 30-40 kHz is the local plasma frequency fpe = cjpe/27r ~ 30 kHz, seen here only as thermal noise with lesser 

instrumental sensitivity than in the CLUSTER observations. The occasional intense dark spots as the one marked by the white arrow 
are brief contacts with the foreshock boundary field line when the beam occurs, deeper in the foreshock the spectrum broadens 

non-symmetrically, evolving 'hair' and 'beards'. Some moving structures can be identified in the hair. The emissions deeper inside the 
foreshock are correlated with the occurrence of low frequency emissions below a few kHz which are probably ion-sound waves [after 1301 . 



spacecraft. Taking into account that the upstream plasma frequency away from the shock is fairly constant this 
implies that the shock-normal angle Qsn changes on a irregular and fast time scale, i.e. that the direction of the 
upstream magnetic field is highly variable even under conditions of apparent quiescence in the upstream medium. 
The conclusion drawn from this would not be that the upstream magnetic field is subject to violent variations as in 
the electron foreshock there is no reason for the magnetic field to be strongly affected. Rather this variability points 
on the temporal variation of the shock area near the tangential field line. Even a small change in the shock normal 
there will cause a large variation in the location of the tangential magnetic flux tube at a distance upstream of the 
shock. Reasons for such a variability have been given in the section on quasi-perpendicular supercritical shocks, but 
to repeat here, the main reason is that the super-critical shock is neither in thermal nor in thermodynamic equilibrium 
and is thus by its very nature subject to changes in all its physical parameters. 

In addition to the plasma frequency emissions, which hint on contact of the spacecraft with the shock-tangential 
field line that makes up the electron foreshock boundary, there are a number of most interesting features seen in the 
wave spectra of Figure[T9] One of them we merely note at this place because we will return it farther below. It is the 
broadband signals that are occasionally seen in many places, in the plasma frequency as signals covering a frequency 
band of up to ~ 20 kHz, in the lower frequencies near the shock, an in particular in the shock transition itself where 
the frequency range of the emission A/ > 80 kHz in all cases exceeds the whole range of the instrument. There are 
no known plasma waves of a comparable bandwidth reaching from ~ kHz deep into the range of free space radiation 
frequencies. Such spectra are, however, known to be produced from narrow spatially localised electric field structures. 
The broader the spectrum the narrower will the structures be. We therefore conclude that in all the places where 
such broadband emissions are encountered, and in particular in the shock transition, we are dealing with very narrow 
localised and intense electric fields which have been generated by nonlinear processes. Such structures are presumably 
solitons and/or phase space holes either of electronic or ionic nature. 

The so far most elaborate investigation of the width of the electron-foreshock boundary beam and its relation to the 
excitation of plasma waves has been performed by [30 and subsequently refined and supported by theory by [65 . An 
example of their ISEE lobservations is shown in Figure[20] Since the instrument was much less sensitive than WHISPER 
the thermal noise at the plasma frequency around fpg ~ 30 kHz appears only as a very faint spotty line. Occasional 
strong narrowband intensification of this line has been identified as contacts with the electron foreshock-boundary 
beam by inspecting the simultaneous electron measurements > 300 eV. These emissions were indeed found to be 
very narrow band, just 1-2 kHz wide at fpe- Deeper inside the foreshock the bandwidth of the emissions broadens 
substantially in a non-symmetrical way. First, it becomes very noisy, consisting of short emissions clumped together 
in groups. Second, long 'hair' emissions evolve of roughly up to 10 kHz bandwidth, while 'beard' emissions are also 
found being generally weaker, but sometimes they extend to low frequencies. Low frequency waves below a few kHz 
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are also detected in connection with these deeper-foreshock emission and were found to smoothly merge from below 
into the beards. 

|30| and subsequently [65j could show that the spectra changed from very narrowband in the foreshock boundary 
flux tube to broadband in the foreshock, just a few electron gyroradii away from the foreshock boundary. They also 
showed that the electron beam distribution flattened over this distance. The waveform of the waves at the plasma 
frequency in the broadband region away from the foreshock-boundary exhibited modulations, which group the waves 
into groups of length of a few 10ms, a frequency roughly corresponding to the frequency of the low frequency waves 
seen in the dynamic spectrum of Figure [20| 

e. The nature of electron-foreshock waves. The narrowband electron foreshock-boundary waves detected at the 
foreshock-boundary flux tube are clearly electron-beam driven Langmuir waves of the surprisingly large amplitude of 
a few mVm"-*^ at the edge of the foreshock and a few % bandwidth in frequency. 

It is considerably more difficult to infer about the nature of the waves deeper in inside the electron foreshock. These 
waves reach large amplitudes of a few lOmVm"-'^ to a few 100 mVm"-'^ and bandwidths of > 30%. (Note that the 
stationary convection electric field in a ^ 500 kms^^ flow in a ^ 5nT magnetic field is just E ^ 2.5mVm^^!) The 
wave electric field is practically parallel to the upstream magnetic field, and the wavelength is long below and short 
above fpe- Moreover, these waves are modulated and seem spectrally to connect to the intense low frequency waves. 
For plasma frequencies of ^ 30 kHz the ion plasma frequency is 0.7kHz. Thus, accounting for the Doppler-broadening 
of the low-frequency wave spectrum in the fast fiow, the low frequency waves can be tentatively identified with ion- 
acoustic waves, which accompany the high frequency waves at the plasma frequency. Obviously these wave modulate 
the latter, it is however not clear whether or not they are created by the high-frequency waves via the modulation 
instability or whether they are generated in a different interaction between the depleted beam and the upstream ion 
flow via an ion-acoustic instability. We have already argued that the modulation instability is unlikely under the 
conditions of the electron-foreshock boundary beam. Moreover, the weak modulation of the wave form of the high 
-frequency wave noted above does not argue in favour of the modulation instability and caviton formation as the 
waves are not really bundled in localised groups and the wavelength is not changed appreciably. On the other hand, 
the spiky broadband nature provides a weak argument for some localisation which can, however, be objected owing 
to the pronounced asymmetry of the spectrum with respect to fpe- 

|65j . ignoring the correlation with the low-frequency waves, interpret the high-frequency waves around fpe on the 
basis of the kinetic Langmuir wave instability including the beam plasma, i.e. referring to a non-gentle beam of finite 
temperature and appreciable density. Their simplified dispersion relation then becomes 
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with Nb, Vb, Vbe the respective beam parameters, density, velocity, and thermal speed, and A^je Ve/wpe is the Debye 
length of the upstream plasma. The first two terms in this expression are the ordinary Langmuir wave. The second 
term on the right is the beam contribution, for small gentle-beam densities the contribution to the real frequency can 
be neglected. However, for larger beam densities we see that Nb/N adds to the unity on the left. Hence the frequency 
of the wave should decrease for small k. Thus long wavelength waves will have frequencies to < tOpe- This agrees 
nicely with the observation. Short wavelength waves, where the term 3/c^A|,g comes into play, will have frequencies 
Lu > ojpe- This is also in accord with observation. Thus, the inclusion of the dense beam into the dispersion relation 
does reproduce the basic observed spectral properties of the beam excited waves at least qualitatively in a simple 
linear way without reference to nonlinear effects like the modulation instability and collapse. The occurrence of 
the low-frequency ion-acoustic waves is then probably due to the electron-ion velocity difference of the hot electron 
component and the cold upstream ion flow via the electron-ion acoustic instability. 

In a series of one-dimensional particle simulation papers, |21| has investigated the dynamics of beams at the 
foreshock boundary in an attempt to quantify the above qualitative theoretical conclusions. He considered a pure 
electron beam on an otherwise neutralising ion background. As expected, a gentle weak beam evolves basically 
according to quasilinear theory into a plateau which for long times remain unchanged. Even in extremely long 
simulation times there is practically no evolution in wave energy, and the distribution function remains stable once 
the quasilinear plateau is formed. This happens approximately after a few hundred plasma periods. For a plasma 
frequency fpe — 30 kHz this amounts to a relaxation time of the order of just ^^Ims. Hence, gentle beams at the 
foreshock should be stabilised within this time, and any nonlinear effects will evolve only at much later times if at all. 
Excluding the ion dynamics, the simulations have been run until ~ 3000 Wp^^ without any susceptible change in the 
distribution and wave level, in physical times until ^ 10 ms. Moreover, the saturation level of the waves is remarkably 
low, much lower indeed than the theoretical estimates suggested. One does not expect that weak or strong turbulence 
effects will evolve which could change or prevent plateau formation. For this to happen, of course, ion dynamics 
should be included into the runs and should produce small k waves, while in the plateau formation the wave number 
increases because the plateau widens and the beam front proceeds to smaller and smaller velocities until the plateau 
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FIG. 21: Results of a one-dimensional electron-electron beam particle simulation [after [211 intended to explain the high-frequency 
plasma waves in the electron foreshock. Left: The electron distribution used. Initially (thin line) it consists of a broad (warm) 
Maxwellian background distribution plus a narrow (cold) electron beam sitting on the main distribution. All quantities are normalised, 
velocity is normalised to nominal beam velocity. At end time tuipe = 1053 (solid line) the beam has become depleted and heated but has 
not disappeared. Right: Real dispersion relations uj(k) (solid lines) and growth rates 7(fc) (thin lines) at different simulation times tujpe 
as indicated by the numbers. Only the unstable part of the real dispersion relations (corresponding to positive 7 > 0) is plotted. The 
Langmuir wave is shown separately with its growth rate labelled /pe which is weekly positive only at very long wavelengths. The main 

unstable waves are shown to propagate in the beam mode (about straight lines starting at origin of uj and k. Maximum unstable 
frequencies are well below the plasma frequency ojpe.With progressing time the wavelength of the unstable modes increases (decreasing 
k). The unstable domain in frequency and wave number shrinks. The unstable frequency decreases as a consequence of the decreasing 

beam velocity which leads to a decrease in the slop of the beam mode. 



is completed and Landau damping stops the further evolution. Therefore, based on these simulations, the dynamics 
of a clod beam cannot explain the richness of the wave observations in the foreshock boundary. 

For a very cold beam the evolution is a bit more complicated as the cold beam initially is not subject to the Langmuir 
instability but rather to the reactive beam instability which generates a broad spectrum below fpe as inferred above 
|65j . The maximum growth is close to but below /pg, and the spectrum is sharply cut off just above /pe- Afterwards 
the beam evolves readily into the kinetic Langmuir stage described above, forming a plateau and stabilising. This 
transition proceeds due to intermediate electron trapping in the wave, which heats the beam to a temperature when 
the kinetic instability can take over. One might thus conclude, that the waves observed below /pe indicated the 
passage of the cold narrow beam front. This, however, is in contradiction to the observation that these waves are 
not observed in contact with the electron-foreshock boundary but deeper in the foreshock where no cold beams are 
present. Hence, another interpretation for the waves is needed. 



One solution is to take into account the bulk streaming distribution [21 . This is shown in Figure 21 Then the 



interaction becomes an electron-electron beam mode interaction with the possibility to destabilise the electron acoustic 
wave mode. This mode has frequency sufficiently far below the plasma frequency, uj < tOpe- A condition is that the 
beam velocity spread is small. Initially the beam also excites frequencies substantially above /pe with large growth 
rate, but these are quickly stabilised. The unstable modes are actually beam modes with resonance condition to < kVb 
which are destabilised by Landau damping from bulk electrons, i.e. the beam speed must enter the range of bulk 
velocities for Landau damping and should thus not be displaced far from the bulk distribution like in the gentle beam 
case. Plateau formation takes very long time in this case such that the wave can reach quite large intensity of the 
order of W/NT^ < 10"'^. For an upstream density of ~ 5 x 10^ m~'^ and electron temperature of Te ~ 100 eV this 
yields an rms electric field amplitude |e|rms ^ lVm~^ in very good agreement with observation. For larger beam 
temperatures the general trend as in Figure[2T] remains valid. It is interesting that the unstable frequency readily 
drops from higher than ujpe to lu < Upe while the wave number shrinks, i.e. the wavelength of the maximum unstable 
waves increases. One should thus observe falling tones in the frequency of the emission. It is also interesting that 
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FIG. 22: Parallel distribution functions in the deep electron foreshock consist of the cold bulk ion distribution and the hot bulk 
electron distribution on top of which the parallel part of the hot diffuse, i.e. warm upstream electron distribution. This distribution in 

the upstream frame moves away from the shock. The upstream electron distribution adjusts to the current-free condition by being 
retarded by the small amount \Ve — Vi\ = Vi,Ni,/N. Its maximum is shifted in the direction of the upstream diffuse beam distribution. 
Such a configuration should be unstable with respect to electron-acoustic waves (dense hot electrons, less dense cool beam electrons). 

during the evolution of the instabihty the slope of the beam mode decreases which simply reflects the retardation of 
the beam during depletion. The final state is low frequency-long wavelength. 

It seems that this calculation explains the observation of electron foreshock wave emissions. However, the problem 
about this model is that it should work only right at the electron-foreshock boundary where only emission at the 
plasma frequency is seen like in the classical case of a gentle warm low density beam. The broadband waves below fp^ 
axe observed inside the electron foreshock, where no electron beams have ever been detected. Moreover, it does neither 
explain the high intensities and the broadband nature of the waves above fpe, nor does it explain the connection of 
the waves below fpe to the low frequency ion-acoustic waves. Therefore we conclude, that probably a two-temperature 
counter streaming electron-component plasma with both components warm will be more appropriate to the interior 
of the foreshock, and probably the cold ion component of the bulk upstream flow must also be included. On the other 
hand, taking these simulations for serious, the deep electron foreshock region should be flUed with many cold and 
not overwhelmingly fast electron beams propagating along the magnetic fleld upstream. Possibly the resolution of 
the current instrumentation is still unable to resolve them. The observations on the right in FigurejlS] might indeed 
indicate the presence of such small beamlets that are distributed over the gyro-angle. So far, however, they can only 
been regarded as measurement fluctuations. 

The parallel configuration of the distribution in the foreshock is shown in Figure[22]for a warm beam that is the field- 
aligned part of the diffuse electron component in the electron foreshock. Because of the vanishing-current condition the 
bulk electron component will be slightly retarded at the small amount |AVe| = ViNb/N . This configuration consists 
of a hot dense bulk electron component and the warm dilute beam component. In addition to the beam instability it 
should be unstable with respect to electron-acoustic waves at frequency oj ~ kve{Ni,/N)2 . Electron-acoustic waves are 
long-wavelength waves below the electron plasma frequency. Thus the low frequency waves find several explanations, 
none of them completely satisfying though. For the high frequency waves which exceed the plasma frequency the 
only available interpretation is that they result from localised electric fields. In frequency space these localised waves 
have a broadband signature. The distributions in the perpendicular direction are non-symmetric halo distributions 
with the backstreaming electrons populating the halo tails. They still await an in-depth treatment for inferring their 
contribution to the wave spectra, particle scattering and plasma heating. 

We close this section by presenting in Figure|23] an average synoptic spectral view of the waves detected in the shock 
transition as was provided by [105 and [45_ from consideration of a large number of shock spectra. There is no new 
information in this figure except that it summarises at a glance the main features in the higher frequency electric and 
magnetic field wave spectra. The cut-off of the magnetic spectra at the electron-cyclotron frequency is no surprise. 
In the electric spectra there is a large temporal variability as can be seen from a comparison of the average and peak 
values (measured within 1 s measuring time) . The variations cover up to two orders in magnitude. Of interest is also 
that the electric spectra exhibit a maximum at the lower-hybrid frequency, an absorption at the electron cyclotron 
frequency, and show a broad maximum at the ion plasma/Buneman two-stream frequency before being steeply cut 
off towards higher frequencies. This cut-off occurs due to Landau damping at the Doppler-shifted frequency of ion- 
acoustic waves when the waves are shifted into Landau resonance with bulk electrons. In this frequency range all 
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FIG. 23: A representative distribution of electric shock spectra. Left: Average (1 s-averages) spectral values, showing the broad peak 
in the electric field spectra and the exponential decay of the spectrum toward high frequencies. For intermediate spectral intensities 
thepeak is well developed, while when the spectral intensity is very high it smoothes out. Centre: The same spectra but in the peak 
values (30 ms resolution) and not in the averages The spectra are similar but much more variable and up to two orders of magnitude 
more intense. This points to the high variability of the electric wave emissions at shocks. In addition indications of the plasma frequency 
are seen at the high frequency end, suggesting that the emissions in fpe are highly time variable [after [TOSl . Right: A schematic summary 
of shock spectra showing the magnetic spectra being cut off at the electron cyclotron frequency fee , and the electric spectra containing 
several maxima at the lower hybrid frequency fi^n at its low frequency end, around the ion plasma fpi and Buneman two-stream /s 
frequencies, with an absorption dip at the electron cyclotron frequency, a steep exponential decay caused by Doppler-shifted Landau 
damping of the waves in the Buneman and ion-acoustic modes at frequency fpeVi/ve = Vi/27rA£i, and the little bump at the electron 

plasma frequency [after [45] . 



kinds of waves overlap, reaching from ion-acoustic waves and electron-cyclotron harmonics to Buneman waves, localised 
structures like Bernstein-Green-Kruskal modes and solitons, while the modified two-stream instability provides the 
large maximum at the lower hybrid frequency. The high variability in the peak values results from the presence of these 
localised structures which recently have been observed [5] in situ but could have been concluded also just from the high 
variability of these spectra in the frequency range of the Buneman-modes respectively Doppler-shifted ion-acoustic 
waves. These are nothing but the signatures of many microscopic phase-space holes that obviously accumulate in the 
shock transition region. Support to such an interpretation has been given long ago by the ISEE measurements of the 
electron distribution function during shock transition which [45 made responsible already for the high variability 
of the peak spectral values. In these measurements the electron phase-space distribution function transforms from 
the upstream streaming Maxwellian to the shock-ramp and downstream flat-topped heated electron distribution, 
which just in the short time interval of crossing the shock -ramp magnetic-overshoot exhibited a clear signature of an 
electron beam that was sitting on the upstream edge of the flat top of the distribution (as depicted in Figure ??). The 
electron beam in that flgure is caused by acceleration in the shock potential and has an upstream directed velocity of 
a few 1000 km s^^. Compared to the bulk electron temperature this beam is cool. Thus it may excite electron-beam 
waves and electron-acoustic waves. However, in combination with the bulk ion flow it is also capable of exciting the 
Buneman two-stream instability which will readily evolve into electron holes and plays a major role in heating the 
plasma and making the main electron distribution flat-top. One should also remember that the scale of the overshoot 
is narrow, of the order of a skin depth, and is most probably due to an intense electron current flowing in the ramp 
or ramp-transition region which, presumably, is also related to the observed beam. 

At high frequencies the spectra also exhibit the small peak caused by the beam excited Langmuir electron plasma 
waves, while the lower frequency foreshock waves are buried in the fat bump of the Doppler-shifted ion-acoustic waves. 
Of course, more cannot be concluded from a picture like this. The more detailed discussion of the spectra requires 
much higher spectral resolution. To the extent as it was available at the time of writing, we have given it here. But 
much more work is to be done until the various modes which can be excited in the foreshock and shock transition will 
be understood better. 
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Radiation 



Shocks are frequently referred to as sources of radiation. Famous examples are supernova shocks, which are visible in 
almost all wavelengths, from radio through visible light up to x-rays [e.g., [19], and solar type II shocks (cf. Chapter 
8) with their main radiation signatures seen in the radio waves. 

Supernova shocks are relativistic shocks which are not treated here. Their Mach numbers range from weakly 
relativistic to highly relativistic, but the energy per particle in them remains to be less than the rest energy of an 
electron, < mgC^ ~ 0.511 MeV, which allows to treat them classically. This does not hold anymore for the central 
supernova engine which drives the flow and which in some cases results in the generation of ultra-relativistic jets. 
There the shocks become non-classical, and not only radiation losses but also particle generation must be taken into 
account in their description. 

/. Observations. The radiation that is occasionally emitted from nonrelativistic shocks is restricted to the radio 
wave range. They do not generate x-rays or visible emissions because, first, of their comparably low energy per particle, 
which is less than the rest energy of an electron <g; rrieC^, second, because of their low energy transmission rate (particles 
are not retarded from flow speed to rest) and, third, because of the low 'emission measure' EM = J dr^AN^{r), i.e. 
the number of radiating particles of density AA^e in the volume that contributes to the emission of x-rays is too low 
for providing any measurable intensity. Moreover, since magnetic flelds are weak and the ratio LUpe/cOec > 1 of electron 
plasma to electron cyclotron frequency is larger than one, gyro-synchrotron emission is unimportant. Thus, the only 
means of how free-space radiation of frequency lo > LUpe can be produced is via plasma wave emission. 

Radio emission from coUisionless shocks in the heliosphere is a widely studied fleld including type II solar radio 
burst, interplanetary type II burst, CME-driven radio bursts and radio emissions from planetary bow shocks. While 
the solar bursts, because of their high frequency can be observed from Earth, most of the other emissions have been 
discovered only from aboard spacecraft. Of the enormous wealth of such observations made by the ISEE, AMPTE, Polar, 
WIND, GEOTAIL and other spacecraft, many of them never published, we pick here just a more recent observation 
from CLUSTER |154j . Figure 24 shows two cases of such observations on December 22, 2000 and three months later on 
March 13, 2001, when CLUSTER was crossing the bow shock and moved into the electron foreshock. On December 22, 
2000 the upstream density was relatively high with plasma frequency near 40 kHz. The intense spots in /pe between 
1200-1300 UT indicate touching of the electron foreshock boundary field line as has been discussed above. Broadband 
wave modes of around ^ \ fpe frequency indicate beam modes though no frequency drift is detectable. However, the 
density during this time is low enough for a faint emission to occur at about 80kHz, just two times fpe- At this 
frequency emission can be only in the electromagnetic free space mode. We are thus witnessing local generation of 
radio emission from the foreshock boundary. At later times when the spacecraft moved deeper into the foreshock - 
as visible from the widening of the plasma wave spectrum to both sides, up and down from /pe - the radio emission 
becomes more diffuse and more broadband shifting to lower frequencies than the second harmonic of the plasma 
frequency. This suggests that possibly here the lower frequency modes close to fp^ do actively participate in the 
generation of emission. On the other hand, radiation from the shock ramp might also contribute to these emissions. 

The lower panel in Figure 24 recorded on March 13, 2001 during a day of much lower upstream density, shows the 
more typical case of Langmuir waves at fpe and radio emission at almost precisely 2 fpe- Note that near 80 kHz a faint 
indication of the presence of an emission at the third harmonic ^ 3/pe can be made out. The dashed vertical lines 
in this panel indicate the times when the magnetic field direction changed abruptly with the changes in magnetic 
field magnitude and density remaining comparably small. The emissions in the plasma frequency and harmonic are 
well correlated during the entire event. Obviously CLUSTER was constantly close to the foreshock boundary as only 
intensification in fp^ is seen but no 'hair' nor 'beards' neither evolve except during a period shortly after 1200 UT in 
the three yellow spots in /pe when the harmonic emission extends to lower frequencies together with the development 
of a little beard. 

Two other interesting features can be read from this panel. At early time in the panel, just before the short active 
phase of the topside sounder at 1020 UT, splitting of the harmonic emission into two narrow bands is seen, which 
reminds at the band-splitting in type II bursts. The other interesting feature is the large number of drifting emissions 
with increasing frequencies, which suddenly evolve right after 1200 UT (following the drop in plasma density at 1140 
UT). Both features, the sphtting of 2/pe and these drifting bursts are not understood yet. The latter might be related 
to the abrupt changes in density and plasma frequency seen in this panel. These density changes are accompanied also 
by changes in the magnetic field, which are not shown here. It is worth noting that these banded drifting emissions 
cannot come from remote simply because the low frequencies arrive first. Radiation from any remote source should be 
visible first at high frequency. The drifting emissions must be related to a nearby source, most probably the shock or 
foreshock. Understanding its generation mechanism should provide valuable information about the radiation source. 

g. Interpretation. Theory of shock-emitted radiation is based on plasma processes which under the prevailing 
coUisionless conditions in the shock and foreshock plasmas refer to wave-wave coupling as the main generation mecha- 
nism. Direct emission from particles is unimportant, because the energy losses a particle experiences when becoming 
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FIG. 24: CLUSTER observations of electromagnetic radiation from the electron foreshock, following the multiple shock crossings 
between 1120-1210 UT that is marked by the intense low frequency noise and its broadband extension through the entire frequency range 
from 0-80 kHz. Upper panel: Near electron-foreshock boundary typical fpe and lower frequency emissions from 1210-1250 UT. These are 

accompanied by a 2 fpe radiation band near 80 kHz and another weak emission band just above fpe the upper cut-off of which is 
decreasing in frequency. At later times the spacecraft is deeper in the shock as indicated by the broadband electrostatic emissions at fpe- 

But these are correlated with broadening of the harmonic radio radiation towards lower frequencies. Lower panel: A sequence of three 
hours of low upstream density when the spacecraft remained close to the foreshock boundary. With the exception of a few short periods 

there is no broadening of the plasma frequency, while all the time a harmonic emission of intensity not much less than at fpe 
accompanies the plasma frequency. During the time of broadening of the plasma frequency one also observes a downward broadening of 
the harmonic radiation. At high frequencies near 80 kHz a faint third harmonic band ~ 3 fpe can be identified. Of interest is the band 
splitting seen in the harmonic radiation between 1000-1010 UT. Finally, the many upward drifting equally spaced in frequency 
narrow-band radiation bands after 1140 UT are surprising. They do not have any counterpart in the plasma frequency and must thus be 
from a remote source that might be related to the small but sharp change in density (drop in fpe) at 1140 UT. The dashed vertical lines 
mark changes in the magnetic field (not shown here), when its direction turned abruptly. Some of these are correlated with small 
variations also in the magnetic field strength [after [T54l and Trotignon, private communication]. 



retarded or reflected at a shock, are not transformed into radiation. In all non-relativistic cases radiation losses can 
completely be neglected compared with all other energy losses. Nevertheless, the observed radiation is of interest 
because in many cases, where no measurements are possible to be performed in situ, radiation is the only direct and 
presumably identifiable signature a shock leaves, when seen from remote. The other signature is the generation of 
energetic particles, which will be treated in the next chapter, but energetic particles are a more diffuse indicator of 
a shock, because their propagation is vulnerable to scattering from other particles, obstacles and, in the first place, 
scattering by magnetic fields. They, thus, do not provide an image of the shock as clear as radiation does. Just 
because of this reason, investigation and understanding of the mechanisms of emission of radiation from shocks enjoys 
- and deserves - the large amount of attention it receives. 
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Expecting that - presumably - direct particle involvement into radiation is improbable (maybe with two exemptions, 
which we will note later) we are left with a small number of possible mechanisms, which all belong to the class of 
wave-wave interaction in weak plasma turbulence. The most probable of these are three-wave processes. These can be 
understood as 'collisions' between three 'quasi-particles'. Since only these three are involved, the interaction conserves 
both, energy and momentum, and can symbolically be written as 

L + L'^T, L = {c.L(kL),kL}, L' = VL(k'L),k'L},T={c.(k),k} (5) 

Here L stands for longitudinal, and T for transverse - as the emitted radiation of frequency uj and wave number 
k propagates in the free space mode and thus is a transverse electromagnetic wave while its two mother waves are 
assumed to be longitudinal (i.e. electrostatic) waves of sufficiently high frequency. Moreover, these frequencies tt'L(kL) 
etc. depend on wave number through the real parts of the electrostatic dispersion relation Di(LUi,\ii) ~ and may 
become quite complicated expressions. 

In the presence of a plasma there are two free space modes, the ordinary and the extraordinary mode. Naturally, 
in order to leave the plasma and escape in the form of radiation their frequencies must exceed some lower cut-off 
frequency, for the ordinary mode 

LU^ - CJ^, + k'c^ (6) 

this is the plasma frequency reached at very long wave lengths fc = 0. However, since the speed of light c 3> f e is 
so large, the radiated wave length is much longer than any of the wavelengths of the longitudinal waves involved. 
This is immediately recognised when comparing the above ordinary wave dispersion relation with the Langmuir wave 
relation, with L = £, 

= Lul^, + ikjvl (7) 

which is of exactly the same structure. Hence, as long as <C c/ ^/i we will have 3> fc, and the wave number of the 
radiated mode is practically zero. Momentum conservation of the three interacting quasi-particles becomes simply 
kL ~ k'|_, implying that the interaction selects counter-streaming electrostatic waves. As for an example, any process 
that is capable of generating Langmuir waves of comparable wavelengths, propagating in both directions along and 
opposite to the magnetic field, can in principle contribute to generating escaping radiation. From energy conservation 
hhJi + huj'^ — hw of the three 'quasi-particles' involved one immediately finds that uj « 2u;pe. This is the origin of and 
the simplest mechanism for the generation of the 2/j,e-second plasma harmonic radiation and has been proposed more 
than half a century ago by |39j as an explanation for the observation of solar type II and type III radio bursts. 

In this simple reasoning we have completely neglected not only the contribution of the Langmuir wave number (which 
turns out not to be important in magnitude, it just shifts the emitted frequency a tiny amount up in frequency) but 
also the fact that the electric field of Langmuir waves is polarised along k^ — zLkg'B/B and thus along the ambient 
magnetic field B, while the electric field of the transverse emitted electromagnetic radiation must necessarily be 
polarised perpendicular to k (because of the absence of space charges at frequencies sufficiently higher than fpe)- 
Since the electric field will, after collision and annihilation of the two Langmuir waves involved, remain to oscillate 
along the magnetic field, the emission is preferably directed perpendicular to the magnetic field. It turns out, then, 
that it is easier to radiate in the extraordinary than in the ordinary free space mode. The extraordinary mode has a 
slightly more complicated dispersion relation; in a dense plasma with Upe S> ujce (as is encountered in near Earth space 
where supercritical coUisionless shocks evolve) it has a slightly higher cut-off frequency. But the argument about the 
smallness of fc ^ fc^ holds also in this case. Radiation at the second harmonic uj « 2ujpe should therefore be polarised 
perpendicular to the magnetic field in the extraordinary mode. 

Unfortunately, foreshock emission has not been found to show any preference in polarisation [103 . Moreover, 
emission is not only in the second harmonic but has also been detected close to > and at the third harmonic 

3wpe . These emissions require different waves to be involved for which a number of mechanisms have been proposed 
[cf., e.g., [14]. So far none of them could be ultimately verified or even agreed upon, each having its merits and 
pitfalls. Since radiation is energetically negligible, as we have mentioned above, the whole problem could be put 
aside. However, since the assumptions made in every radiation mechanism contain important information about the 
source region, the problem of radiation production in coUisionless shocks remains to be tantalisingly urgent and awaits 
resolution. 

Since there is no agreement yet about the radiation mechanism, we merely note here some of the different proposals. 
The first is the above mentioned merging of two counter streaming Langmuir waves. The problems about this simple 
though suggestive mechanism are numerous. First, Langmuir waves are assumed to be generated by the gentle-beam 
instability. Ignoring the problem of beam survival during its propagation along the shock-tangential field line, which 
we have discussed already in detail, gentle beams excite only forward Langmuir waves, which requires some mechanism 
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that backscatters a substantial percentage of waves and inverts the direction of their wave numbers. There are three 
known elegant processes that are capable of doing this: modulation instability respectively collapse, scattering of 
Langmuir waves off thermal ions, and scattering off ion-sound waves, all three proposed long ago [for an early review 
of the latter two mechanisms cf., e.g., |156j. 

Modulation instability generates ion-sound waves via the ponderomotive pressure force Fpmf — —{e^ /mf.ui^^)V\eg\'^ 
of the high-frequency Langmuir wave, eg{v,t). These waves, when becoming locally large amplitude, structure the 
plasma into a chain of cavities in which the Langmuir waves become trapped. This process generates long wave 
lengths. It is described by the Zakharov equations for the combined evolution of the Langmuir wave field and the 
density variation 5N , respectively. 

The first of these equations is a driven wave equation for the density variation which for slow time variations, when 
the derivative with respect to time is neglected, just gives pressure balance between the ponderomotive pressure on 
the right and plasma pressure on the left, i.e. proportionality 5N ^ In other words, the density variation 

anti-correlates with the field pressure, which corresponds to caviton formation. The second equation is a nonlinar 
Schrodinger equation for the evolution of the wave amplitude. 

The Langmuir waves trapped in the cavitons must bounce back and forth, which naturally creates counter-streaming 
waves of equal intensity with opposite wave numbers. During collapse the cavities shrink in size, the wave numbers 
and momenta of the waves increase, and the wave energy density increases as well because of the shrinking volume. 
This yields both, the counter streaming Langmuir waves being localised in the same region and, in addition, a large 
radio emissivity. Unfortunately, we have already noted it, this process - as beautiful as it might be - has not been 
confirmed experimentally, neither in the observations nor in the full particle simulations. Observed wave intensities 
are too low in the electron-foreshock boundary and electron foreshock, and the density variations did not indicate the 
presence of the expected cavities. Simulations, on the other hand support quasilinear evolution and wave scattering 
off thermal ions. We note, however, that the most recent detection of the very strong electric fields in the shock ramp 

might indicate that it is not the electron foreshock where one should expect caviton and collapse to work and cause 
the most intense radiation, rather it might be the very shock transition where shock radiation is generated by such 
processes. It is, in this respect, most interesting to remind of the strange radiative behaviour reproduced in Figure [24| 
that was detected by CLUSTER. We also note that similar observations had been made much earlier with the wave 
experiment on AMPTE IRM in the Earth foreshock (R. A. Treumann & J. LaBelle, unpublished 1986). The observed 
band splitting and high intensities might have been caused by Langmuir caviton collapse |153] . 

Other possibilities to produce counter streaming Langmuir waves are scattering of Langmuir waves off thermal ions 
|84l investigated this process in full detail numerically including the ion polarisation cloud] , a mechanism known since 
the early sixties. The process reads symbolically L + i ^ L' -I- i*, where the primed quantities are after the collision, 
and the star on the ion indicates excitation of the ion as it is too heavy for changing momentum during the collision 
with the Langmuir wave. It is merely excited while the scattered Langmuir wave has changed direction and lost some 
of its momentum, i.e. attains a longer wave number and lower frequency. The same process does also work with 
ion-sound waves as L+ IS — > L'. The scattered Langmuir waves then also change direction by absorbing the ion sound. 
Both processes have been used for radiation generation |165j . 

Radiation at higher frequency, e.g. radiation at the third plasma harmonic can be generated by a four-wave process. 
This is also favoured by caviton formation and collapse since the waves in this case are all confined to one and the 
same volume. However, other mechanisms have also been proposed. All these processes are of the kind of wave-wave 
interactions and thus their efficiencies are proportional to the product of the involved relative wave intensities. Since 
the latter are usually low, the efficiencies are very small as a rule. An attempt to increase the growth of Langmuir waves 
has in the recent past been the idea to consider a statistical theory of growth called 'stochastic growth' [5nill04j . This 
attempt takes advantage of the statistically distributed density fiuctuations in the foreshock region like in a random 
medium. Since the Langmuir-wave growth rate is proportional to SN/N, an average growth rate over the volume 
of occupation by the Langmuir wave can be calculated. This might be more realistic than using the linear growth 
rate. Regions of decreased density contribute strongly to the average Langmuir-wave amplitude. In this way an in 
the average larger emission efficiency is obtained. Moreover, the averaging procedure introduces a statistical element 
which supports the incoherence of the relation between the detected Langmuir waves and radiation, which is in partial 
agreement with the observation. 

Generation of radiation at the fundamental co > in a three wave process requires the presence of a low frequency 
wave. Ion-acoustic waves are one possibility, other possibilities are lower-hybrid waves, Buneman waves, the modified- 
two stream instability, various kinds of drift waves, and also electron acoustic waves or electron beam waves. In 
particular the latter are present in the foreshock region and thus can combine with Langmuir waves to generate 
fundamental radiation slightly above the plasma frequency. 
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The two exemptions when particles become involved are the above mentioned scattering of Langmuir waves off 
ions, and the so-called electron-cyclotron maser instability [for a contemporary review see, e.g.. ,152j . Its advantage is 
that it operates directly on the free space mode avoiding any intermediate step like three-wave processes or particle 
scattering. However, it requires a particular form of the electron distribution with a velocity space gradient in the 
perpendicular direction dFe{v\\,v_\_) / dv^_ > 0, a hot electron distribution and low cold electron density. It is barely 
known whether such distributions are realised at the shock. However, if they are in some place, then the cyclotron 
maser instability will outrun all other mechanisms and directly feed the free-space electromagnetic radiation modes. 
Radiation will then be at a harmonic of the electron cyclotron frequency which is a severe restriction if the magnetic 
field is low and the density high. Therefore, regions of low density and stronger converging magnetic fields are the 
best candidates for this radiation source. 



III. QUASI-PARALLEL SHOCK REFORMATION 



In quasi-parallel supercritical shocks there is not such a stringent distinction between the region upstream of the shock 
and the shock itself like in quasi-perpendicular shocks. The foreshock, which we have discussed in some detail in the 
previous section, and the shock itself cannot be considered separately. This is due to the presence of the refiected and 
diffuse particle components in the foreshock. These, as we have seen are the source of a large number of waves. 

The interaction of these waves with the shock is one of the main issues in quasi-parallel shock physics. In the present 
section this will become clear when we will be dealing with the formation, behaviour and structure of quasi-parallel 
shocks as it has been inferred less from observation than from numerical simulations. The reason is that the real 
observations in space do not allow to separate the particles and waves from the shock. They all occur simultaneously 
and are interrelated and can never be observed in their initial state. The observations to which we will nevertheless 
occasionally refer will leave the impression of large-amplitude noisy fiuctuations. In simulations, on the other hand, 
it is at least to some degree possible to prepare the system in such a way that a single effect can be studied. For 
instance, in one-dimensional simulations the direction of wave propagation can be prescribed which allows studying 
just waves in one direction and their effect on the shock and particles. Moreover, treating the electrons as a neutralising 
Boltzmannian fiuid suppresses their effect on the ion motion and wave genereation. Treating them as an active fiuid 
allows taking account of electron-ion fiuid instabilities. Finally, full particle PIC simulations can be performed with 
low or realistic mass ratios in order to investigate different time and spatial scale dependence and the excitation or 
coupling to higher frequency waves. Most simulations that have been performed in the past have taken advantage of 
these possibilities. 



A. Low-Mach number quasi-parallel shocks 



It is usually assumed that low-Mach number shocks are stable, i.e. show no substantial time variation or reformation. 
It is not completely transparent why this should generally be so. Firstly, the critical Mach number has been shown 
by |54| to become small at narrow shock normal angles Qbu ~^ 0; in which case even low-Mach number quasi-parallel 
shocks should become supercritical and reflect ions. Secondly, any fast ions of parallel velocity v\\ > V\ that have 
been heated in the shock can in principle escape from the quasi-parallel shock upstream along the magnetic field and 
should appear in the foreshock where they contribute to wave generation. Therefore, it makes sense to investigate the 
state of quasi-parallel shocks in view of their stability and wave generation even for low Mach numbers. In addition, 
any waves that are generated in the shock ramp or transition with upstream directed k and fast enough parallel phase 
or group velocities could also escape from the shock in upstream direction. This could, in particular, be possible just 
for low Mach number quasi-parallel shocks. 

To check this possibility [92j have performed one-dimensional hybrid simulations finding that initially the quasi- 
parallel shock consisted of phase-standing dispersive (magnetosonic) whistler waves with the last whistler wave cycle 
constituting the shock ramp. As expected, the wave vectors of these phase-locked magnetosonic whistlers are aligned 
with the shock normal. At later times, backstreaming ions along the upstream magnetic field excite a long-wavelength 
whistler wave packet upstream of the shock. In the one-dimensional simulation the wave vector is restricted to the 
shock normal while it is known from theory that the growth rate is largest along the magnetic field. These oblique 
whistler waves should thus show up in two-dimensional simulation and may be visible at larger distance for sufficiently 
large upstream phase velocites. 

Such two-dimensional (hybrid) simulations with a non-inertial electron fluid have been performed by |120] for a 
Mach number of Ma = 2.2 and angles = 20°, 30°, 45° and by [20] for an angle Qsn = 30° and high Mach 

number M.a = 5 (Figure 25 1. The lower Mach number simulations show the presence of a substantial number of 
backstreaming ions which cause an ion-ion instability in the upstream region. However, the excitation and properties 
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FIG. 25: Two-dimensional hybrid simulations of the evolution of upstream waves. Left column: Upstream wave in a low-Mach number 
Ma = 2.2 quasi-parallel shock [after 220] ■ The nominal shock is in the (y, z)-plane. Ions escaping to upstream generate the oblique 
upstream magnetosonic waves. The contour plot of the two normalised to the upstream magnetic field components of the magnetic 
fluctuations shown is taken at relatively early times uidt = 68. It shows the nearly plane magnetic wave fronts inclined against the shock 
in direction x and having wavelengths of ~ lO/ojpi in z while being much shorter in x. In the vicinity of the shock the wave fronts turn 

more parallel to the shock and produce a non-coplanar magnetic component \by\ which is of same order as the \bz\. Moreover, even 
though the shock has relatively low Mach number, it is not completely stable but shows structure in z direction which is produced by the 
presence of the reflected upstream particles and the upstream waves. At places it is impossible to identify one single shock ramp. Right 

column: Two-dimensional hybrid simulations of the evolution of giant magnetic pulsations (SLAMS) in front of a quasi-parallel 
supercritical shock [after [20l . Top: The simulation plane showing the structure of the (normalised) magnetic fluctuation field |b| at time 

tujci = 45 and shock normal angle @Bn = 30°. The accumulation of the growing wave fronts at the shock transition, their increasing 
amplitudes, and their turning towards becoming parallel to the shock is clearly visible from the rotation of the two wave fronts and their 
k vectors shown in white. Away from the shock transition the angle between k and the shock normal n is large. Close to the shock the 
two vectors are about parallel. The magnetic field is in the wave front, so 0b„ is close to 90° here. Bottom: Pulsation amplitude and ion 

phase space. The fluctuations evolve into large amplitude pulsations when approaching (and making up) the shock. The strong 
retardation of the upstream flow by the pulsations is visible in the shock-normal velocity component (Mach number Ma)- In the hybrid 

simulations this slowing done is accompanied by some ion heating. 

of the waves depend strongly on the shock normal angle. Initially, as in the case of [HI] phase-locked short-wavelength 
whistlers appear which are replaced at later times by upstream long-wavelength whistlers with phase velocity directed 
and amplitude growing towards the shock ramp but upstream directed group velocity, i.e. the shock radiates energy 
away towards upstream, as one would naively expect, because the shock being supercritical must reject the excess 
inflow of energy which it can do by both, reflecting particles and emitting waves into upstream direction. These waves 
are excited by the backstreaming ion component in a strongly nonlinear interaction process because of the evolving 
steep ion-density gradient, which is of the same scale as the whistler wavelength. The k-vector turns away from the 
magnetic field having comparable components parallel to B and parallel to the shock normal n. For small Qbu a 
remnant of the initial phase-locked whistlers survives but disappears &t Qbu = 30°. Close to the shock, where the 
backstreaming ion density is high, the waves have short wavelengths, and k is almost parallel to n. In the high Mach 
number simulations no shock is produced but instead reflected ions were artificially injected with same Mach number 
as the incoming flow but with much higher temperature Vi = 14.1 Va, forming a spatially uniform ion beam. The 
intention was to investigate the effect of the hot reflected ions. This is shown on the right in the above figure. The 
result resembles the former one where a shock was generated by reflection at a wall, but the effect in the injected beam 
case is stronger because of the higher Mach number. Hence it is the hot reflected ion component that is responsible 
for the wave dynamics and the shock dynamics. 

All this can be seen from the two-dimensional intensity contours of these waves in the foreshock, which are plotted 
at a relatively early time in the shock evolution in the simulations in Figure [25l On the right in this figure the geometry 
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FIG. 26: The evolution of the shock normal angle ©b„ on distance from the shock in two-dimensional hybrid simulations for two 
initial quasi-parallel shock-normal angles 0Bno — 2° and &Bno — 20°, respectively [simulation results taken from 12Q;. The horizontal 
line at 45° is the division between quasi-perpendicular and quasi-parallel shock normal angles. In both cases thetabn evolves from 
quasi-parallel direction into quasi-perpendicular direction. The shaded areas identify the quasi-perpendicular domains. 



is given, with the magnetic field fluctuation vector b in the (y, z)-plane. The bottom panel on the left shows contours 
of the bz fluctuations in the (a;, z)-plane. The nominal shock ramp is at a; « 145 Xi ion inertial lengths Xi = c/cupi at 
this time. The upper panel shows the non-coplanar component 6j^-contours in the same representation. Behind the 
shock the fluctuations are irregular and disorganised. However, in front of the shock a clear wave structure is visible 
with strongly inclined wave fronts and of roughly ~ 10 A.; wavelengths in z parallel to the nominal shock surface. The 
wavelength in x is about three times as short. 

These waves are seen in both components, bz,by, are low amplitude at large distance from the shock but reach 
very large amplitudes simultaneously in both components during shock approach while, at the same time, bending 
and assuming structure in z-direction that is different from the regular elongated shape at large distance. This 
deformation of wave front may be due to the residual whistlers near to the shock, but it implies that the shock has 
structure on the surface in both directions x and z and is not anymore as planar as was initially assumed. The shock 
becomes locally curved on the scale of the shock-tangential wavelength. The waves deform the shock and, in addition, 
being themselves of same amplitude as the shock ramp, become increasingly indistinguishable from the shock itself. 
The shock is, so to say, the last of the large- amplitude magnetic wave pulsations in downstream direction, and the 
shock-magnetic field is not anymore coplanar, because the waves have contributed a substantial component by that 
points out of the coplanarity plane. 

We have emphasised this phrase, because it expresses the importance of the low-frequency upstream magnetic 
waves in quasi-parallel shock physics. Contrary to quasi-perpendicular shocks where the reflected gyrating ions in 
combination with the reflected-ion excited modifled-two stream instability were responsible for the shock dynamics 
and different kinds and phases of shock reformation, quasi-parallel shock reformation and much of its physics is 
predominantly due to the presence of large-amplitude and spatially distinct upstream waves. These are the generators 
of the shock and, due to their presence, the shock changes its character. It is highly variable in time and position along 
the shock surface and is - close to the shock transition on a smaller scale ~ 'less quasi-parallel' (or more perpendicular, 
i.e. the shock-normal angle Qsn has increased on the scale of the upstream waves). The latter is due to the out- 
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of coplanarity-plane component of the upstream waves. In spite of concluding this from a hybrid simulation, this 
conclusion remains basically valid also in full particle simulations. It had been suggested already earlier on the basis 
of ISEE l& 2 observation of the evolution of the upstream ultra-low frequency wave component . 

The gradual evolution of the shock normal angle 0Bn has been demonstrated in other hybrid simulations by [120 and 
|20j who investigated the evolution of the shock normal angle in dependence on distance from the shock. This is shown 



in Figure 26 for two-dimensional hybrid simulations with initial shock-normal angles Qbhq = 2° and Qbho = 20°, 
respectively, which we anticipate here. In both cases Qbu evolves from quasi-parallel to quasi-perpendicular angles. 
Qualitatively there is little difference between the two cases. At the shock ramp Qsn is deep in the domain of quasi- 
perpendicular shocks. The only difference is that for the nearly parallel case the angle jumps to quasi-perpendicular 
quite suddenly, just before approaching the shock ramp, while the evolution is more gradual for the larger initial &Bn- 
In both cases the evolution is not smooth, however, which is due to the presence of large-amplitude foreshock waves. 
Transition to quasi-perpendicular occurs for the initially nearly parallel case at the nominal shock ramp while for the 
initially quasi-parallel case it occurs at an upstream distance of about < lOOAi from the shock. 

One notices that this transition is on the ion scale, implying that in the region close to the shock the ions experience 
the shock occasionally (because of the large fluctuations in &Bn) ~ f^nd when ultimately arriving at the shock - as 
quasi-perpendicular. It is thus not clear, whether the electrons do also see a quasi-perpendicular shock, here. However, 
the ISEE measurements of the electron distribution function by |34j at the shock do not show a difference between 
quasi-perpendicular and quasi-parallel shocks. This fact suggests, in addition that, close to the shock transition, 
quasi-parallel shocks behave like quasi-perpendicular shocks as well on the electron scale, which is just what wc have 
claimed. 

We finally note that the behaviour of the shock normal angle gives a rather clear identification of the location of the 
shock transition in the quasi-parallel case, as indicated in Figure|26]by shading. Three distinctions can be noticed: 

• first, at larger initial shock-normal angles the transition to quasi-perpendicular angles occurs earlier, i.e. farther 
upstream than for nearly parallel shocks. This is due to the stronger effect of the large amplitude upstream 
waves in this case; 

• second, at larger initial shock normal angles the quasi-perpendicular shock transition is considerably broader 
than for nearly parallel shocks, i.e. it extends farther downstream before the main quasi-parallel direction of 
the magnetic field in the downstream region takes over again and dominates the direction of the magnetic field: 

• third, at an initial shock-normal angle of 20°, this region is roughly ~ 150Ai wide, implying that the magnetic 
field direction behind a quasi-parallel shock remains to be quasi-perpendicular over quite a long downstream 
distance measured from the shock ramp. For the nearly parallel shock this volume is only about ~ SOA^ wide. 

This observation must have interesting implications for the physics downstream of quasi-parallel shocks. For instance, 
applied to the Earth's bow shock, where Ai ~ 10"^ km. both distances correspond to regions wider than the order of 
> 5 Re which is larger than the nominal width of the magnetosheath! Thus, behind the bow shock a substantial part 
of the magnetosheath plasma should behave as if the bow shock would have been a completely quasi-perpendicular 
shock. 



B. Turbulent reformation 



When speaking about turbulent reformation we have in mind that a quasi-parallel supercritical shock is basically 
a transition from one lower entropy plasma state to another higher entropy plasma state that is madiated ba a 
substantially broad wave spectrum. Such a transition has been proposed by |126j based on the detection of the large 
amplitude magnetic pulsations (SLAMS) in the foreshock. 

Figure 27 on the left shows their model assuming that somewhere upstream in the foreshock magnetic pulsations 
have been excited which become convected downstream toward the shock by the convective flow, grow in amplitude 
and number and accumulate at the shock transition to give rise to a spatially and temporarily highly variable transition 
from upstream of the shock to downstream of the shock. An important clue in this argument was the observation that, 
first, the pulsations grow in amplitude when approaching the shock and that, second, they slow down, this slowing 
down is effectively an increase in their upstream directed velocity on the plasma frame with growing amplitude such 
that their speed nearly compensates for the downstream convection of the flow. 

On the right of the figure, which is suggested by the observations of [S7], a larger volume is seen. Here the pulsations 
are the result of growing ultra- low- frequency waves which are generated in a volume inside the foreshock but relatively 
close to the ion-forcshock boundary. These waves grow to large amplitudes until evolving into pulsation which the 
flow carries toward the shock. Growth, slowing down, and accumulation then lead to the pile up of the pulsations 
at the shock location and formation of the turbulent shock structure. Of the magnetic fleld, in this figure on the 



37 




FIG. 27: The patchwork model of |126| of a quasi-parallel supercritical shock mentioned earlier. Left: Magnetic pulsations (SLAMS) 
grow in the ion foreshock and are convected toward the shock where they accumulate, thereby causing formation of an irregular shock 
structure. Note also the slight turning of the magnetic field into a direction to the shock normal that is more perpendicular, i.e. the 
magnetic field is more parallel to the shock surface with the shock surface itself becoming very irregular [after [1261 . Right: The same 
model with the pulsations being generated in the relatively broad ULF-wave-unstable region in greater proximity to the ion-foreshock 
boundary. When the ULF waves evolve to large amplitude and form localised structures these are convected toward the shock, grow, 
steepen, overlap, accumulate and lead to the build up of the irregular quasi-parallel shock structure which overlaps into the downstream 

direction. 



right, we plotted only the shock-tangential upstream field line. In the left part, several field lines are schematically 
shown exhibiting the fluctuations imposed by the background level of ultra- low- frequency fluctuations. Moreover, 
a certain bending of the field lines is included here in approaching the shock transition with the field lines turning 
more perpendicular the closer they come to the shock. This bending is what we claim to be a parallel shock turning 
quasi-perpendicular at a scale very close to the shock. In this schematic drawing, however, there would be no reason 
for the magnetic field to turn this way. What closer observations and simulations show is, however, that the turning 
of the field is the result of the presence of the large amplitude magnetic pulsations. This will become clearer below. 



Observations 



Of course, the model shown in Figure 27 is a schematic model only which, however, has some merits in explaining 
the observations. The signature of a quasi-parallel shock in the magnetic field is indeed quite different from that 
of a quasi-perpendicular shock. We have already seen in the electric recordings reproduced in Figure[24] that the 
quasi-parallel shock appears in the electric wave spectrum as a broadband emission of highest spectral intensities 
at the low frequency end. The magnetic signature of a quasi-parallel shock is quite similar in that it lacks a clear 
location of the shock front. Rather one detects a broad region of very large amplitude compressive oscillations in 
magnetic magnitude and in the direction of the magnetic field that subsequently is recognised as a passage across the 
quasi-parallel shock. 

An example is shown in Figure 28 as measured by the CLUSTER spacecraft. This figure shows eight hours of 
observation by CLUSTER in the immediate vicinity of the quasi-parallel shock. It is difficult to say where in the 
figure the shock transition is located as the large fluctuations in the magnetic field magnitude and directions mask 
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FIG. 28: Eight hours of CLUSTER magnetic field data during a long passage near and across the quasi-parallel supercritical (Alfvenic 
Mach number Ma ~ 12 — 13, ion inertial length ~ 140 km) bow shock. The time resolution was 4 s. The top panel shows the 
variation in the magnitude of the magnetic field. The two lower panels are the respective elevation and azimuthal angles 9,(j) in a GSE 
coordinate frame [data taken from I71| . Large variation in the magnetic compression and direction can be seen to be associated with this 
quasi-parallel shock crossing. Buried in these large variations on this highly time-compressed scale are many magnetic pulsations 
(SLAMS). The compressive large amplitude fluctuations in the upper panel are typical for a quasi-parallel shock transition. 



the various back and forth passages across the shock that are contained in the data. Clearly, at the beginning near 
1400 UT the spacecraft was in the downstream region. The fluctuations show that during almost the entire sequence 
the magnetic field is exhibits compressive fluctuations. These belong to the shock transition. At the same time large 
fluctuation in the direction of the magnetic field are also observed. In the compressions of the magnetic field buried 
are also upstream pulsations (SLAMS), and many of the changes in direction belong to the ultra-low-frequency waves 
present at and near the shock. The changes in direction indicate that the shock does not behave like a stationary flat 
surface. Instead, it shows structure with highly fluctuating local shock normal directions. 

|71j have checked this expectation by determining the local shock-normal angle &Bn and comparing it with the 
prediction for Q sn estimated from magnetic field measurements by the ACE spacecraft which was located farther out 
in the upstream fiow. The interesting result is that during the checked time-interval of passage of the quasi-parallel 
shock the prediction for the shock normal was around 20 — 30°, as expected for quasi-parallel shocks. However, this 
value just set a lower bound on the actually measured shock normal angle. The measured Qbu was highly fluctuating 
around much larger values and, in addition, showed a tendency to be close to 90°. This is a very important observation. 
It strengthens the claim that quasi-parallel shocks are locally, on the small scale, very close to perpendicular shocks, a 
property that they borrow from the large magnetic waves by which they are surrounded. In fact, we may even claim 
that locally, on the small scale, quasi-parallel shocks are quasi-perpendicular since the majority of the local shock 
normal angles was > 45°. By small-scale a length scale comparable to a few times the ion inertial length or less is 
meant here. 

The data suggest that, indeed, the quasi-parallel shock is the result of a build-up from upstream waves which 
continuously reorganise and reform the shock. Figure|29] shows three representative examples of such upstream waves 
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FIG. 29: CLUSTER magnetic field measurements of magnetic pulsation (SLAMS) near and remote from a quasi-parallel shock on 
February 2, 2001 [data from l7ll . Only the magnetic field magnitude is shown for all four CLUSTER (colour coded) spacecraft. The 
spacecraft separation was between a few 100 km and 1000 km. Top: Clustered pulsations in the shock transition. Three events of large 

amplitudes are shown. These structures are very irregular with steep fronts. Note that in spite of the small spacecraft separation the 
shapes of the structures differ strongly from spacecraft to spacecraft. Moreover, determination of the pulsation fronts and normals (not 
shown) indicates high variability over the spacecraft separation distance. Thus the structures are of relatively small scale and large 
amplitude. Middle: Isolated pulsation at greater distance from shock outside compression region. The structures are seen almost 
simultaneously at the spacecraft and thus must be of larger size. Copared to the embedded pulsations the amplitudes are lower, and the 
structures are more regular. Bottom: A shcklet observed outside the pulsation region in the domain of ultra-low-frequency waves. The 
steep shock-like front is well expressed with the qattached whistler waves it carries with it. Note the much lower amplitude than the 

pulsations. 



which are far from being continuous wave trains. The upper panel is taken from the large density fluctuation region 
in the shock transition. This region turns out to consist of many embedded magnetic pulsations (SLAMS) of very large 
amplitudes. In the present case amplitudes reach |b| ~ 25 nT. These pulsations have steep flanks and quite irregular 
shape, exhibit higher frequency oscillations probably propagating in the whistler mode while sitting on the feet or 
shoulders of the pulsations. It is most interesting that the different CLUSTER spacecraft - at spacecraft separation 
< 1000 km - do not observe one coherent picture of a particular pulsation. This implies that the pulsations in the 
shock are of shorter scale than spacecraft separation: the different spacecraft observe different structures respectively 
different pulsations. In addition, the magnetic field directions (not shown in the figure) are very different from 
spacecraft to spacecraft and thus from pulsation to pulsation, and even for one pulsation at its front edge and trailing 
edge different magnetic field directions are observed. The directions of the magnetic normal across a pulsation change 
on very short spatial scales. The quasi-parallel shock front has thus a rather irregular shape, which will be bent locally 
with changing direction of its normal. 

The second panel shows an isolated pulsation farther away from the shock transition. This pulsation exhibits a 
much more coherent way on the four different CLUSTER spacecraft. It seems as if the pulsation is still in evolution 
as three of the spacecraft see a nearly coherent structure while the fourth which is farther away sees it in a different 
state. Concluding from this event, isolated pulsations seem to have larger dimensions and lower amplitudes, which 
would be consistent with the assumed solitary properties of pulsations. 
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FIG. 30: Observation of electric field structures in large magnetic pulsations (SLAMS) in the quasi-parallel shock transition region 

[data taken from[7l[8]. Structures on three different time scales are shown, corresponding also to three different spatial scales. Top: 
CLUSTER passage across on (moderately large amplitude) magnetic pulsation in the shock transition. The (smoothed) magnetic field 
structure is a slightly steepened magnetic bump. The stationary parallel electric potential field across this structure shows a potential 
ramp with steep gradient at the leading edge of the pulsation. The potential drop of ~ 400 V corresponds to an electric field of 
~ 0.47, Vm~^. Note that the time scale in this panel is 90 s. Bottom left: Six seconds of a CLUSTER passage through the shock 
transition. The black dots show the spacecraft potential variation which maps the local density variation. Overlaid is the high frequency 

WHISPER trace of the plasma frequeny. In the magnetic pulsation regions (white) the plasma frequency exhibits huge excursions to 
both sides similar to those seen in the overview Figure [24| on wave observations. These excursions trace the BGK (nonsymmetric) modes 
and (symmetric) solitons. Bottom right: One example of one of the solitons on a 2 ms time scale. It is nicely seen how symmetric the 
parallel potential trough and the corresponding bipolar parallel electric field shape look like in the solitary wave structure. 



The third panel shows a shocklet, i.e. a structure which presumably has little in common with pulsations. It is 
embedded into long wavelength ultra-low-frequency wave trains, evolves into steep front and drives whistlers attached 
to this front across the flow. These waves were already observed by [110 . Their properties indeed resemble those 
of sub-critical little shocks which propagate against the flow, though with slower speed such that they effectively are 
slowly convected towards the shock. 

[7 report another interesting property of magnetic pulsations in the shock transition region where they overlap 
to form the quasi-parallel shock. Measurement of the electric cross-SLAMS potential identify a substantial unipolar 
drop in the electric potential of several 100 V, corresponding to a potential ramp, when passing from upstream to 
downstream across the pulsation. Such a drop signifies the presence of an electric field in one direction across the 
pulsation. Taking the mean size of a pulsation to be roughly 1000 km, the mean electric field is {E) ~ 400 mVm^^. 
However, this field drops mainly at the leading edge of the pulsation. Such a field presumably corresponds to a steep 
pressure gradient in the pulsation. It could also be generated by an anomalous collision frequency. This remains to 
be tested by further observations and data analysis. 

The measurements of [7] anticipated the more recent report of strong electric fields in the shock by [5J . In addition to 
this observation it was found [8 that the single pulsations were subject to a fairly large number of high frequency/Debye 
scale structures in the electric field seen in the WHISPER recordings (bottom panel), which belong to electron holes 
or solitons which form in the pulsation gradient regions as shown in Figure [30l The bipolar electric field and unipolar 
potential across one - indeed very symmetric ~ soliton is seen in this figure. These observations suggest that the 
pulsations are indeed the main constituents of a quasi-parallel shock with the dynamics on the micro-scale of a 
quasi-parallel shock going on mainly in the single pulsations of which the shock transition is built. 

The occurrence of these intense nonlinear electrostatic electron plasma waves at the quasi-parallel shock transition 
is intriguing. It forces one to draw another very important main conclusion from these observations (and other related 
observations like those of 0): that quasi-parallel shocks are sources of electron acceleration into electron beams, which 
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are capable to move upstream along the magnetic field over a certain distance and excite electron plasma waves at 
intensity high enough to enter into the nonlinear regime, forming solitons and electron holes (BGK modes). 

Presumably this is possible only when in the supercritical quasi-parallel shock transition region the magnetic field 
changes from quasi-parallel to quasi-perpendicular on the electron scale Xe- Indications of such a change on the 
ion scale have been noted above at a number of occasions, but the detection of solitary structures in the electron 
plasma waves in relation to quasi-parallel shock transitions provides a very strong argument for this to be true on a 
scale which is well below the ion scale. Only if this is the case, there will be ample reason for electrons to become 
reflected and accelerated into beams from the transition region in a quasi-parallel shock. As we already noted, we 
may, therefore, expect that quasi-parallel supercritical shocks on the electron scale are not anymore quasi-parallel 
but change to become locally quasi-perpendicular, while on the larger ion scale they still maintain properties of 
quasi-parallelity. 

If this conclusion will turn out to be true and will sustain future more sophisticated experimental tests, it will have 
important consequences for coUisionless shock physics. Supercritical coUisionless - nonrelativistic - shocks will, in 
fact on the electron scale, always behave quasi-perpendicularly - and it may be suspected that this conjecture will 
also hold for relativistic shocks though probably for other reasons (like the generation of transverse magnetic flelds 
by the Weibel instability, which becomes dominant in relativistic shocks [see, e.g., IHlllSSj)- This implies also that the 
true quasi-parallel shock physics cannot be properly elucidated when ignoring electron effects as is, for instance, done 
in hybrid simulations. 

Simulations of quasi-parallel shock reformation 

Nevertheless, before turning to the very few more realistic full particle PIC simulations with different mass ratios from 
small to about realistic, we are first going to discuss here hybrid simulations in one and two dimensions of the formation 
and behaviour of quasi-parallel shocks. Initial particle |100L IIOII I102j and hybrid simulations (TUl 11181 1119L 11421 1160] 
did already illuminate some of the particular interrelations between the dynamics of quasi-parallel shocks: energy 
dissipation by short wavelength whistlers at the shock transition, the presence of a diffuse hot ion component upstream 
of the shock, and the importance of upstream waves. Most of these have been the subject of reviews [the interested 
reader might be directed to the papers by TTJ [T31 EH which more or less systematically discuss selected aspects of 
these interactions]. The selection of results we will provide here is guided by the progress that has been achieved in 
the understanding of the quasi-parallel shock physics and in its relation to the observations. 

The first arising question is, whether nearly parallel shocks are at all capable of reflecting ions. This question 
is neither nonsensical nor academic. It makes sense, because the investigation of the critical Mach number by |54j 
becomes unreliable at small angles, say Osn ^ 30°, because it is based on linear whistler dispersion and does not 
take into account the completely modified plasma conditions near the shock ramp. Thus the simple conclusion that 
the critical Mach number approaches A^crit — *■ 1 for &Bn might be an unjustified extrapolation. Theoretically, 
the refiection of particles at small G_b„ must become entirely due to the electrostatic shock potential drop with the 
magnetic part of the Lorentz force being obsolete. The potential drop, however, should decrease with M, and thus the 
reflection of particles should cease, while being of vital importance for the generation of a coUisionless quasi-parallel 
shock. Hence, the question is also not purely academic. 

Of course, we do already know from the observations that quasi-parallel shocks exist at small Qbu with their 
foreshocks being populated by a diffuse ion component that excites upstream waves and mediates the beam-generated 
upstream ion-foreshock boundary waves. The impossibility for this diffuse component of being entirely due to shock 
reflection in the quasi-perpendicular part of the shock, immediately proves that the quasi-parallel (or even the nearly 
parallel shock) must be able to reflect particles upstream. Hence, either a quasi-parallel shock is capable of generating 
a large cross-shock potential, or it is capable of stochastically - or nearly stochastically - scattering ions in the shock 
transition region in pitch angle and energy in such a way that part of the incoming ion distribution can escape 
upstream, or - on a scale that affects the ion motion - a quasi-parallel shock close to the shock transition becomes 
sufficiently quasi-perpendicular that ions are reflected in the same way as if they encountered a quasi-perpendicular 
shock. 

Observations suggest that the latter is the case, while observations also suggest that large potential drops occur in 
the large-amplitude magnetic pulsations (SLAMS) where they accumulate in the shock ramp [8^. Hence, reflection of 
ions will be due to the combination of both effects, the electric potential drop and the magnetic deflection. In fact, 
this can be a quite complicated process for an ion passing across a number of magnetic pulsations, in each of which 
it is being retarded and at the same time deflected by a small angle until its normal velocity component is decreased 
sufficiently that a further deflection in pitch angle suffices to let it return into the upstream region. 

h. Hybrid simulations in ID. In agreement with what is known today, flrst one-dimensional hybrid simulations 
in an extended simulation box [ID] suggested that the reformation of quasi-parallel shocks is about cyclic and is 
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FIG. 31: Fast mode dispersion relation in simulations in the shock frame. Left: The Doppler shifted fast mode dispersion relation in a 
supercritical flow in the shock frame. The Dispersion relation assumes negative frequencies corresponding to the downstream convection 

by the flow of Mach number Ma = 2.8 and at @Bn = 30°. Waves at zero group velocity have energy at rest in the shock frame. 
Negative group velocities imply downstream transport of energy, positive group velocities imply upstream transport. Right: Simulated 
upstream wave dispersion spectra near the shock and upstream of the shock. Near the shock wave energy accumulated around standing 
and downstream transport. Away from the shock the wave energy still moves upstream [data from i61| . 



caused by the impact of large-amplitude upstream waves. |119j . using Qbu — 20° and A4a — 3.5 with an upstream 
ion thermal velocity Vi = va (corresponding to Pi = 1) in one-dimensional hybrid simulations (with small numerical 
resistivity) showed that the reflected ions are not coming from the core of the incident upstream ion distribution but 
originate in the shell of this distribution, having initial velocities v > 1.7vi. These ions escape from the shock quite 
far upstream and excite ultra-low frequency waves with upstream directed velocity of ~ 1.3Va at distances up to 
> SOOAi, which are convected downstream to reach the shock. In this one-dimensional hybrid simulation the only 
mode in which they can propagate is the compressive fast magnetosonic mode. 

These waves are in fact what in observations has been identified as pulsations (SLAMS) but is not yet recognised 
as such here. During downstream convection the waves grow and slow down in the interaction with the foreshock 
ion component. When approaching the shock they generate a large amount of new reflected ions. These slow the 
incident ion population down and steepen the wave crest, which becomes the new shock front. In the time between 
the arrival of the compressive waves the shock is about stationary and develops phase-locked upstream whistlers 
which the arriving next wave crest destroys. From these simulations it could not be concluded what process produced 
the reflected ions, however, as one-dimensional simulations among suffering from other deficiencies select only one 
particular direction of wave numbers and are thus not general enough for drawing final conclusions. 

The nature, generation and effects of the large-amplitude upstream waves have been further investigated in more 
detail in one-dimensional j6T | 111?) I118[ among others], and in two-dimensional hybrid simulations [2.0., .62, 117,. Since 
shocks are three-dimensional, it is clear that two-dimensional numerical simulations at same resolution come closer 
to reality. However, they suffer from restrictions in size of the simulation box and simulation time. Since reality 
does not confront us with an initial state, large boxes and long times are needed. However, for the investigation of 
particular questions, simulations have the great advantage of providing the possibility to prepare them for answering 
just those questions. Concerning the propagation of upstream waves in one direction with respect to the shock normal, 
one-dimensional simulations are just good enough. 

In order to identify a particular wave mode, the dispersion of the wave must be investigated. This dispersion relation 
depends on the frame in which it is taken, because the energy/frequency of a wave is not invariant with respect to 
coordinate transformations; in a medium moving with velocity V it is Doppler shifted according to uj' = u!(k) — k • V, 
where w(k) is the dispersion relation in the rest frame of the flow. While the Doppler shift at high frequency is 
negligible, it completely changes the dispersion of ultra-low frequency waves at large Mach numbers. Figure 31 on 



its left shows the deformation of the fast mode dispersion relation in the shock frame at large Mach number > 1 
and for waves propagating upstream in the plasma frame. The deformation causes negative frequencies of the waves 
which imply downstream directed phase velocities, which is nothing else but the intuitive downstream convection of 
the waves by the flow. However, the minimum in the dispersion relation implies that waves of a particular frequency 
and wave number have zero group velocities. In the shock frame the energy of these waves is stationary. Smaller wave 
numbers have energy moving downstream, larger wave number have energy moving upstream away from the shock. 
The right part of the figure shows simulations of upstream waves according to one-dimensional hybrid simulations by 
[61] for Qbu = 30° and a Mach number A4a = 2.8. The entire dispersion of the simulated waves is negative. The 
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waves are all convected toward the shock as their Mach number is less than the streaming Mach number. Near the 
shock most of the wave energy moves downstream and will cross the shock. Still some shorter wavelength waves (large 
k) move in energy upstream in the shock frame. Farther away from the shock most of the wave energy encountered 
is seen to move upstream. 

|117| investigated these waves further in onc-dimcnsional hybrid simulations performing several numerical experi- 
ments on them, taking away the shock and instead injecting a diffuse ion component from downstream. The main 
finding is that the large amplitude upstream magnetic pulsations (SLAMS) evolve out of the ultra-low frequency 
wave spectrum in the interaction with the diffuse ion component. In accord with observation the pulsations move 
upstream in the plasma frame. Thereby their upstream leading edge steepens and is right-hand circularly polarised 
like required for whistlers. However, dispersion is unimportant; the main cause of the evolution of large pulsations is 
nonlinearity when the wave interacts with the diffuse ion distribution. This distribution has a steeper shock directed 
density gradient than the pulsation wavelength. Moreover, the flow becomes decelerated at the leading edge of the 
pulsation (as is seen in Figure 25), and here the velocity difference between the flow and the diffuse ion component 
drastically decreases, which shifts the k vector of the resonant wave to larger values, and the wavelength decreases 
during convection of the pulsation toward the shock. (Note that no resonant ion beam-whistler interaction exists as 
the beam is hot and diffuse.) The standing whistlers are in the leading edge are simply generated by the current 
flowing in the steep edge. It is thus concluded that it is the gradient in the hot diffuse ion component over a length 
of the same order as the length of the wavelength which produces the pulsations. Ultimately these pulsation cause a 
quasi-periodic reformation of the shock, as we have described earlier. This is thus proved by one-dimensional hybrid 
simulations. The same result is obtained when the simulation starts right away without a shock but with an injected 
hot beam (which is no surprise as the generation of the shock, before it was removed in the former simulations in 
order to keep with the wave field, was due to the plasma flow-reflected ion beam interaction). 

i. Hybrid simulations in 2D. The two-dimensional evolution of the pulsation (SLAMS) was studied later by |120] 
and POJ (see Figure 25 ) . It basically confirmed the conclusions drawn from one-dimensional simulations with the three 
important modifications, first, the wave fronts of the pulsations (SLAMS) rotate into a direction that is more parallel 
to the shock thereby increasing the shock-normal angle locally to become quasi-perpendicular; second, the pulsations 
have short wavelength in shock normal direction, but are of substantially longer but finite lengths in the direction 
tangential to the shock, which provides structure to the shock in tangential direction; and third, shock reformation is 
a result due to the steepening and accumulation of the pulsations and is a quasi-periodic process but the downstream 
structure of the shock over some distance is caused by the downstream convection of the old shock front, i.e. the bulk 
of the pulsations that had accumulated at the location of the former shock transition. It is interesting to note that from 
reformation cycle to the next rather large fluctuations in the magnetic field and density exist in the transition from 
upstream to downstream which may be capable of trapping particles. The two-dimensional simulations do also confirm 
the conclusions that the diffuse upstream ion component is responsible for the growth of the pulsations (SLAMS). The 
simulations by |20j . in particular, followed the same scheme as the one-dimensional simulations, injecting hot diffuse 
ions into upstream in order to control the interaction between the diffuse ion component and large amplitude pulsations 

(slams). 

j. Full particle PIC simulations. So far we dealt just with hybrid simulations where the ions are macro-particles 
while the electrons represent a charge-neutralising background of zero mass. Clearly, such simulations are unrealistic if 
whistlers become involved. This is, however, the case, as we have discussed above, whenever large amplitude pulsations 
evolve at the leading edge of which phase-locked whistlers are attached. The question what role the electrons play 
in the evolution of the pulsations can only be answered by full particle simulations. These require large simulation 
boxes and at the same time high temporal and spatial resolutions. So far they could therefore only be performed 
in one dimension [SH 11221 1155] . [94] used a mass ratio of mi/rrie — 100. They had to small a box (just ~ 30Xi in 
the upstream direction) for following the evolution of upstream waves but stressed the importance of whistlers in 
shock reformation. One decade later it became possible to substantially enlarge the box and at the same time to 
switch to a larger mass ratio while staying with one dimension only. |122j , using the same mass ratio at Q Bn — 30° 
and Ma — 4.7, had an upstream entension of ^ 200Ai and could follow the shock evolution for a time tojd ~ 100. 
Similar simulations with mass ratio mi/rrie = 50 have been performed by |155j . Here we discuss in greater detail the 
calculations of |122j . An overview of their results is shown in Figure 32 in the fixed lab-frame for the main magnetic 
field component B^, electric shock potential <&, bulk stream velocity along the shock normal V , and density N , all as 
functions of distance x (measured in electron inertial lengths Ae). 

Far upstream from the quasi-parallel shock the magnetic field exhibits long wavelength ultra-low frequency waves 
(4) which, when approaching the shock, start steepening at their leading edges (note that in the plasma frame these 
waves are moving upstream, as also does the shock, i.e. to the left in the figure). The amplitude of the wave increases 
(3) , and the wave becomes a pulsation (SLAMS). Its amplitude is large enough to already substantially brake the 
upstream bulk flow, which causes a drop in V and in <I>, and an increase in density. The pulsation has slowed down 
the flow to a velocity ^ 2.8Fa already here. In fact the leading edge of the pulsation behaves like a 'baby shock', 
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FIG. 32: Full particle PIC simulations of the evolution of a quasi-parallel shock in one dimension only [after [T22l . From top: main 
magnetic field component B^, electric potential $, bulk plasma flow velocity V, density A'^, all in simulation units. The numbers indicate 
three pulsations (SLAMS). Pulsation 1 was the old shock. Pulsation 2 is the actual shock coinciding with the drop in V to zero and the 

steep increase in density and potential. Pulsation 3 is just evolving. It already has a steep leading upstream edge and decelerates the 
plasma flow. It will become the next shock ramp. Number 4 indicates a bump in the upstream waves that will become a pulsation. The 

actual shock ramp has some phase- locked whistlers attached to it which are not well resolved on the scale shown (see the next figure). 



which later will become the real 'adult shock'. The shock itself is formed further downstream at the local position 
of the leading steep edge of the previous pulsation (2). In front of this edge (i.e. the genuine shock at this instant 
and location) between it and the trailing edge of pulsation (3) a standing phase-locked large amplitude whistler has 
evolved. This whistler is spatially damped by the approaching pulsation (3). An indication of such a whistler has 
already been seen in front of the leading edge of pulsation (3) as well. Pulsation (2) (the instantaneous shock) has a 
substantial downstream extension. Further downstream of it the 'old shock' is seen, which was formed at an earlier 
time by pulsation (1); and even farther downstream a remainder of earlier shock ramps (pulsations) is recognised in 
the trace of B^. The instantaneous shock ramp (pulsation 2) is high enough to completely brake the upstream flow, 
the velocity of which drops to zero while the density steeply increases and forms a dense wall. 

Figure|33] gives and impression of the shock evolution in higher temporal and spatial resolution. It shows instanta- 
neous magnetic profiles in a box of lengths 400Ze including the shock ramp with time axis running upward. In this 
simulation frame the shock moves upstream to the left. The heavy step- like line indicates the approximate location 
of the nominal shock ramp. Also shown is an upstream ultra-low frequency wave that at time tujci — 82.5 has begun 
to evolve into a pulsation while convectively approaching the shock, growing in amplitude and developing a steep 
leading (upstream) ramp in front of which phase-locked magnetosonic whistlers start growing. When approaching the 
shock at tuJci = 92.5, the pulsation kills the phase-locked whistlers that were waiting in front of the ramp by damping 



45 




FIG. 33: Full particle PIC simulations of the evolution of a quasi-parallel shock in one dimension [data taken from 11251 . The profile of 
the main magnetic component is shown for subsequent simulation times shifted by AtWd = 2.5 upward. The representation is in the 
simulation frame, i.e. the shock moves to the left into the upstream direction. The simulation shows the reformation of the quasi-parallel 

shock resulting from the exchange of the shock with an incoming upstream wave which has steepened to become a pulsation (SLAMS). 

The magnetic field trace at time tujd = 100 has been overlaid on the ion phase space at this time. The heavy steps show the location of 
the nominal shock ramp (where the flow is stopped). It moves slowly upstream until a new pulsation arrives and when it suddenly jumps 
forward by roughly lOOAe. Also shown is the fast approach of an upstream pulsation starting at tujci = 82.5 and arriving at the shock at 

tcOci = 92.5 to take over the role of the shock. Note that in the minima of the 'whistler' fleld fluctuations (at tuici = 100) ions are 
trapped, oscillating back and forth and forming hole vortices in phase space centred around local minima of the electric potential <I> (not 
shown). The little boxes indicate where particle (ion or electron) phase space distributions have been determined. 



them out. At this time the pulsation takes over the role of the shock, and the nominal shock position jumps ahead 
to upstream by roughly a distance of the width of the pulsation ~ 100 Ag. This process repeats itself at tuJci = 102.5 
showing that the quasi-parallel on this time scale shock is not stationary but undergoes nearly periodic reformation 
which is mediated by the arrival of large amplitude pulsations. 

In order to illuminate what happens in the shock ramp, the ion phase space data have been overlaid at time 
tuJci = 100 just before a new reformation. The shock ramp coincides with the location where the upstream flow 
is completely braked and the hot downstream ion distribution begins. This location corresponds to the leading 
edge/ramp of the first broad downstream pulsation. On the other hand, the phase-locked upstream whistlers that 
are attached to its ramp are seen to trap ions, which form large phase space vortices (holes) with ions bouncing back 
and forth in the associated potential wells. These are not shown in the figure but resemble the magnetic oscillations, 
being minimum in the centres of the holes. However, what is really interesting to note is that there is little indication 
of any back streaming reflected ion beam in this simulation. Instead, a very dilute but broad completely diffuse ion 
component is seen all over the simulation box, which is denser downstream than upstream but of similar velocity 
dispersion. Somehow this component has been generated by the shock, and it is which is responsible for the growth 
of the pulsations. 



The evolution of the proton and electron distribution functions is shown in Figure 34 The two proton (ion) distri 



butions are taken in the phases when fully developed whistlers with trapped ions and phase-space vortex formation 
exist, and when these whistlers have been completely damped by the trapped component. In the former case at 
time tWci = 90 the (perpendicular) ion shows a substantial heating with respect to the upstream ion distribution, 
has been slowed down by a factor of roughly ~ 8 in velocity and, in addition, exhibits a nonthermal top-flat part in 
the upstream direction which is due to the presence of the trapped bouncing ions in the phase space vortices. At 
time tuJci = 95 when the whistlers are completely damped and the pulsation has taken over to become the shock, 
the trapped vortices have disappeared, as has been shown above, and the distribution evolved into a broad and hot 



46 



o 



I I I I r 

Protons tcD =90 

xa.= 1040-1090^ upstream. 



xyT. = 1050-1 100 - V \ \ 



f■'JLA.'^A'^^/^ bouncing ' i 

trapped ions ~ 
h, "^^^-On vortices -'rx^y^ ^ ^ 



n 

Electrons 



1000-1050 ' 



, ^ upstream^ 

^'^'vi^ __fi eld-aligned 
beam _ 



upstream 
nai 




y \ V evolv ing 

' \ \ flattop 
f(v ) * \ \ 



4 r .^f.'i^^ ^ 



1 010-1050 



o 



-6 



-80 







80 



FIG. 34: Velocity distribution functions fi{vx) for the full particle PIC simulations of Figurc |33| during different phases of shock 
evolution as indicated by the boxes in that figure [data taken from [122]. Left: Proton distribution functions at times tud = 90 and 
fojci = 95. The first period is in the arriving pulsation before it takes over to become the shock. Shown is the part of the pulsation where 

it encounters the whistlers. The upstream distribution has changed here into a heated distribution exhibiting a substantial upstream 
directed tail which is due to the trapped ion component. The whole distribution assumed a large velocity spread and is strongly heated. 
At time tujci = 95 the trapped ions have damped the whistlers. The distribution has now become very hot but has a two-temperature 
structure with high-energy tail in the leading part of the pulsation. At this stage the incident plasma is completely slowed down. Right: 
Two electron velocity distributions fe{vx) at time twd = 92.5 and fe{vz) at time tuici = 95, respectively. The former is right in the new 
shock built up of a fresh pulsation showing a hot but otherwise structureless perpendicular electron distribution, the latter is in the 
well-developed shock but along the main magnetic field. This parallel distribution is strongly nn-symmetric, heated, exhibits an 
upstream tail, indication of the formation of a top flat distribution and a remaining upstream-beam-likc part similar to those 

distributions as had been measured by | 34 |. 



distribution at about zero flow velocity. Thus the vortices are entirely due to the interaction between the trailing edge 
of the shock-arriving pulsation and the whistlers with the latter being compressed, amplified in amplitude until they 
trap ions, while the pulsation itself is free of those trapped ions. This was already clear from Figure[33] but its effect 
on the distribution is nicely seen here. 

The right part of Figure 34 shows two cases of electron distributions. The first in the perpendicular direction, the 
second in the parallel direction taken in the arriving pulsation when it becomes the shock, i.e. after the reformation. 



The boxes are given in Figure 33 The perpendicular distribution fe{v^) simply shows the heated, about symmetric 
electron distribution in the shock. However, the parallel distribution fe{vz) exhibits a number of interesting features. 
It is also hot, in fact substantially hotter than the perpendicular distribution, exhibiting nonthermal tails in both 
directions along the magnetic field with a stronger upstream tail. In addition, on the downstream side it shows the 
evolution of a broad flat top on the distribution with a remainder of an non-flattened upstream directed beam on the 
upstream side. This is ineresting as it is almost exact the type of distributions that had been reported by [33] from 
crossings of the super-critical quasi-parallel bow shock [see also [351 US] • 

k. Conclusions. The full particle simulations provide much deeper insight into the shock physics than did the 
hybrid simulations. Still, they have been only one-dimensional and thus are restrict to waves propagating solely 
into shock normal direction. Moreover, the mass ratio is still not realistic as it competes with the length of the 
box, and it is obvious that the length of the box is crucial in quasi-linear shock physics for properly investigating the 
evolution of the waves which build the shock, and for their interaction with the particle component. It is, in particular, 
of considerable interest that, even though the magnetic field should allow the streaming of particle upstream in a 
quasi-prarallel shock, no upstreaming beams have been seen in the simulations. Electrons exhibit a signature of an 
upstreaming beam, but ions do not. The diffuse component which is homogeneously distributed over the entire region 
depending only on the distance from the shock is obviously participating in a diffusive process and is not injected into 
the upstream region. This sheds light on the acceleration and confinement process acting in quasi-parallel shocks. 

These simulations have demonstrated that quasi-parallel supercritical shocks are nonstationary on the time scale of 
the simulations. They are subject to reformation. However, this reformation is quite different from the reformation 
of a quasi-perpendicular shock as it is not due to the presence of gyro-bunched reflected ions in the foot of the shock. 
Rather it is caused by the accumulation of large amplitude magnetic pulsations which are not separate entities but 
grow out of the upstream ultra-low frequency waves in the range where the wavelength of the wave is of the same 
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order as the density-gradient scale of the diffuse ion population in the foreshock. |120j . [20] and |122j have shown 
that extracting the diffuse ion component results in a lack of upstream pulsations, while injecting additional diffuse 
ions speeds the evolution of large amplitude pulsations up. Thus, pulsations are the result of the fine-tuning of the 
resonance of ultra-low frequency upstream waves with the diffuse ion component. These pulsations play an important 
role in quasi-parallel shock. Still, their dynamics has not been completely clarified and, in addition, is subject to some 
controversy. Closer investigation of pulsation dynamics is therefore of vital interest. 

To a certain extent this has been the subject of the one-dimensional full particle PIC simulations performed by 
|155j who use a modest mass ratio rrii/me = 50 [see also£8j which, however, allows them to extend the simulation 
time scale to times as long as tiOd ~ 38 and to enlarge the simulation box to ~ 700 A^. Their simulation is based on 
the magnetic piston method instead of the reflecting wall that has been used by |122j . In general these simulations 
confirm the previous results, adding some facts about the structure, size and life times of isolated pulsations whose 
sizes are determined to lie in the interval between (10 — 20)Ai with the tendency of a gradual shrinkage during the 
reformation process. This shrinkage can simply be attributed to the dynamical compression of the pulsations at the 
shock transition where many embedded pulsations compete for the available shock volume. 

Also confirmed is the rotation of the magnetic field from quasi-parallel to quasi-perpendicular in the reformation 
process. Large amplitude pulsations (SLAMS) behave like localised supercritical quasi-perpendicular shocks oj which the 
quasi-parallel shock is constructed. We stress the importance of this conclusion here again. Moreover, their evolution 
is accompanied by the generation of phase-locked whistler precursors at the leading edge of a pulsation, as has been 
shown above to be the case |122j . |155| report in addition the observation of large-amplitude spiky electric fields in 
the leading edge of a pulsation. In view of the observations by [5], their investigation, however, will require even much 
better time resolution on the electron plasma scale. Structures of this kind require the inclusion of the full electron 
dynamics and imply accounting for realistic mass ratios in the simulations. 

Inclusion of particle dynamics is one of the main problems in quasi-parallel shock physics. We have seen that the 
evolution of the building blocks of a quasi-parallel shock, i.e. the large amplitude pulsations (SLAMS), is completely 
determined by the interaction of the ultra-low upstream wave spectrum with the diffuse ion component. On the other 
hand, the pulsations build up the quasi-parallel shock transition region which is responsible for the generation of the 
diffuse ion component in a way which has not yet been satisfactorily clarified. The current knowledge will be discussed 
in the next paper in this series. 

One idea is that the quasi-perpendicular subshocks of which the quasi-parallel shock is constructed - and which are 
nothing else but the various pulsations of which the shock consists on the small scale of a few - reflect sufficiently 
many ions back upstream. However, because of the complicated magnetic structure of the shock transition region, 
these ions cannot simply escape along the magnetic field to upstream but remain trapped for long enough time in 
the shock transition, being scattered in pitch angle and energy until ultimately picking up enough energy in order to 
leave the shock either upstream or downstream in a diffusion process that works in the shock and partially also in the 
foreshock regions. The observation of a nearly isotropic diffuse ion distribution both upstream and downstream of the 
shock in the simulations and the lack of observation of ion beams, both in the measurements and in the simulations, 
provides a strong argument for a diffusive mechanism to act in the shock transition. The observations by |57j also 
support this conclusion as they show that the diffuse ion component exhibits an exponential density gradient towards 
the shock with the shock being their source. Such a gradient is typical for a diffusive process. 

Whether electrons undergo a similar process, remains an unresolved question that will be attacked with the increas- 
ing capacities of computers when the full electron dynamics can be included into the simulations. There are strong 
hints on the importance of electron dynamics in the presence of the observed spiky large amplitude electric fields in 
the pulsation and shock transition. These electric field structures point on the generation of electron beams. Since 
no such beams have been observed in the deep upstream foreshock, the structure of the shock on the small scale 
presumably inhibits the escape of electrons to far upstream in the form of beams. This can be due to either a chaotic 
magnetic field configuration on the electron scale which makes the escape of electrons along the magnetic field impos- 
sible. It can also be due to the trapping of electrons in the phase space holes which are related to the smallest-scale 
electric field structures. These, on the other hand, require strong electron currents to flow in the shock transition 
and, if present, will stochastically heat the electrons. Moreover, they might be responsible for the electromagnetic 
radiation generated in the shock ramp, a process that has been observed in solar and interplanetary type II radio 
bursts but never found a convincing explanation by any known mechanism. The finding in the PIC simulations of 
the formation of top-flat electron distributions and a rudimentary upstream directed field aligned beam in the shock 
transition |122j might indicate that the electrons are accelerated in the embedded pulsations by electric fields and 
can form a field aligned upstream beam which apparently cannot leave far upstream but possibly has sufficient free 
energy in the electrons for driving either a Buneman or modified two-stream instability which causes electron plasma 
waves, electron holes, and radiation - and at the same time scatters, heats, and confines the electrons. 
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IV. HOT FLOW ANOMALIES 



Would the coUisionless upstream flow be stationary at constant Mach number M > Merit, we could safely terminate 
the discussion of the supercritical shock transition regions at this place and focus our attention on the downstream 
part of the shock transition. Unfortunately, this is not the case. The upstream flow is usually in a highly disturbed 
fluctuating state with the smallest disturbances in it readily evolving into structures, which can appear in various 
form, either as turbulence on almost all scales, from the electron scales up to scales of a substantial fraction of the 
macroscopic dimension of the flow, as localised large amplitude waves, as shocklets, or as various kinds of boundaries 
in the flow, which appear as discontinuities and current layers. 

Then the question arises, what happens when any such forms impact on the shock and interact with it. As long as 
we are dealing with the interaction of a single wave with the shock one needs to investigate the conditions for wave 
reflection, refraction and transmission across the shock. This is most easily done for wavelength much longer than the 
shock transition scale and leads to modiflcations of Snell's law when taking into account the Rankine-Hugoniot jump 
conditions at the shock, for instance, when the MHD approximation holds. For shorter wavelengths the continuous 
change in the wave propagation conditions as given in the dispersion relation can be parameterically treated. The 
shock may then act as a spectral filter depending on which waves it allows to pass. 



A. Observations 



Of quite a different kind of problem is the interaction of a current sheet or discontinuity boundary with the shock. 
Current sheets and discontinuities are large amplitude distortions of the upstream flow and cannot be treated like 
waves interacting with a shock. Observations near the Earth's bow shock have shown that they lead to severe 
distortions of the shock, which have been termed 'Hot Flow Anomalies' (HFAs, the term now commonly used even 
though it does not appropriately refer to the physics involved), 'Active Current Sheets', 'Diamagnetic Cavities' or 
'3D-Plasma Structures with Anomalous Flows'. They have been discovered in spacecraft transitions of the Earth's 
bow shock from analysing AMPTE and ISEE measurements [3ll [97l I123[ 11241 1146} 1161(1162) . 

The first observation of the HFA phenomenon was made by |123j in the plasma and field observations of AMPTE UKS 
who detected an unusual very hot upstream plasma event resembling a broad layer of plasma in which the upstream 
magnetic field performed a large rotation as it is known from current sheets. The most surprising observation was, 
however, that the upstream fiow inside this current sheet was sharply defiected by a substantial angle close to ^ 90° 
from its original upstream direction, a defiection far in excess of the defiection the bow shock would provide and in 
addition not expected to occur in front of the shock in its upstream region. The upstream fiow was about stopped 
from its original fiow direction. This defiection was also different from the one predicted by the magnetic stresses 
across the boundary of the hot fiow anomaly, and it was clear that it required a violent momentum exchange between 
the upstream fiow, the shock, and some (unknown) agent. It was also found that the structure was not in its final 
state but was evolving. Hence, most probably it did not have come all the way along from the source of the upstream 
fiow, the Sun, and it was speculated that it was either caused by an upstream-fiow current sheet interacting with the 
supercritical bow shock or also could have had its origin in a violent change of the configuration of the obstacle, i.e. 
the magnetosphere in this case, an idea that was advocated later on [97 . It took about one decade of measurement 
and theoretical investigation until a consensus was reached that the HFA phenomenon was a shock phenomenon and 
was not caused by the obstacle. Clearly, violent changes in the obstacle configuration do also affect the shock [97 lll36j 
but most probably to a lesser extent than it happens to be observed in hot fiow anomalies. 

The signatures of a Hot Flow Anomaly in the upstream foreshock plasma are: (1) a strong defiection of the bulk 
fiow velocity, (2) reduced magnetic field strength, (3) reduced density, and (4) considerable heating of the plasma. 
The region where all these changes in the plasma parameters are observed is usually fianked by walls of enhanced 
plasma density and magnetic field strength when both, the upstream plasma and upstream magnetic fields pile up. 

All these properties are nicely seen in the most recent observation of a Hot Flow Anomaly by THEMIS plotted in 



Figure 35 Similar observations by the CLUSTER spacecraft have been reported by [72], but while all four CLUSTER 
spacecraft were upstream of the shock, the THEMIS spacecraft were both upstream and downstream such that they 
could follow the evolution of the HFA from upstream to downstream. 

The most spectacular effect in the observation of a Hot Flow Anomaly is seen in the fifth panel on the left of 
this figure. The velocity component Vx turns positive in the HFA, which corresponds to a complete stopping or even 
reversal of the fiow in the downstream direction. The other two fiow velocity components are enhanced simultaneously, 
indicating a defiection of the fiow in the (y, z)-direction. 

Inspecting the magnetic field, one finds that its magnitude drops to a small value. However the traces of its 
transverse components become dipolar, which is the sign of a current that is fiowing inside the HFA. To both sides 
of the HFA current both the magnetic field and density pile up. It is interesting that this pile up is not symmetric 
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FIG. 35: THEMIS spacecraft measurements of a Hot Flow Anomaly upstream and downstream of the quasi-parallel Earth's bow shock 
wave [data taken from [28]. Left: The THEMIS A observations in the foreshock close to the shock transition. Three regions of the HFA 
have been shaded: darker shading refers to the pile-up regions at the HFA edges, The magnetic depression region in the centre of the 
HFA is shaded light. The quasi-parallel shock transition is about 1 min before the HFA encounter. The deflection (retardation in Vx) of 
the upstream flow is clearly seen. It coincides with a drop in density and magnetic field and strong plasma heating. Bipolar shape of 

Bz, By indicates the presence of a current sheet. Note also the asymmetry of the HFA pile-up edges. The upstream edge which is 
exposed to the flow is much stronger piled up. These regions are themselves shocks. The lowest panel is the power spectral density of 
magnetic waves at different frequencies. Right: The simultaneous THEMIS E observations downstream of the shock behind the shock 
transition. The signature of the HFA is very irregular here. The later onset confirms that the HFA is an upstream phenomenon. 



but is stronger on the upstream side. This should be so since the HFA must stop the upstream flow, which causes a 
stronger compression on the upstream than on the downstream sides. In the figure the two pile up regions are shown 
in shading. They flank the current region where the plasma is dilute and hot, as is seen from the strong increase in 
temperature and from the ion energy flux in the uppermost panel which exhibits a broadening of the distribution 
that is comparable to the shock transition at the beginning of the panel. However, though the temperature is high, 
the mean energy of the ions is low, reflecting the strong deceleration of the flow. 

The right part of the flgure shows simultaneous measurements downstream of the shock. These indicate that the 
HFA downstream is very complex, causing a strong disturbance in the downstream flow. It, moreover, sets on later 
than in the upstream region which is a conflrmation for the HFA being an upstream phenomenon and is not caused 
by the obstacle or downstream of the shock. This rules out number of models that had been proposed in the past. 

We will not go into the details of this downstream observation other than point on the much more diffuse shape 
of the HFA after its "passage" across the shock to downstream. It has become dissolved into a number of different 
regions which still belong together while each differs from the others. The larger of these regions are numbered from 1 
to 4 in the flgure. Some of them showing very clear signatures (correlated magnetic plus density compressions, regions 
1 and 4) of fast modes, others of slow expansion modes (anticorrelated magnetic and plasma variations plus plasma 
heating, region 2). In any case, it becomes clear that the shock has split the HFA into different fllaments. The HFA 
has lost the compactness, consistence, and coherence it possessed in the foreshock and during the interaction with 
the shock. It would contradict the imagination (and the second law) if such a diffuse downstream structure would 
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FIG. 36: Schematic of the interaction of an upstream tangential discontinuity with the shock. The tangential discontinuity is 
convectcd along with the normal component Vn of the upstream flow which, in the shock frame, is thus the same on both sides. 
However, the tangential velocities and magnetic fields are different on both sides of the discontinuity. The latter is due to the inclined 
current flow J in the discontinuity. Here the rotation of the magnetic fleld is shown in such a way that the normal components of the 
electric fleld E„ point inward on both sides of the discontinuity thus confining any shock-reflected ions to the interior of the discontinuity 
which here is shown as an extended region. Due to the presence of the discontinuity and interaction with the shock the upstream flow is 
deflected from its original direction, because the tangential velocity has changed across the discontinuity. For other current direction the 
tangential magnetic field would rotate differently. E„ can then have other directions, pointing away from the discontinuity on both or 

also only on one side. 

combine to make up a coherent HFA of the kind we see in front of the shock. In particular, its later downstream 
than upstream onset, which causes a time delay between the upstream and downstream disturbances, provides a very 
strong causal argument against a downstream origin of the HFA, at least in this particular case and for the class of 
HFAs that exhibit similar properties. 

We note that [148] investigated a large number of Hot Flow Anomalies and determined the directions of their normal 
vectors, applying the minimum variance method. They found that the normal vectors are in good agreement with 
the assumption that upstream tangential discontinuities are involved into HFAs. This assumption since dominates 
the interpretation and theory of Hot Flow Anomalies. 

B. Models and simulations 

The mechanism of hot flow anomalies is believed to be related to the interaction of the shock with an upstream 
current sheet [TH I143j . most probably a tangential discontinuity, in the near-shock quasi-parallel foreshock. It is 
this model that is most strongly supported by the recent THEMIS observations shown in Figure[35] and also by the 
CLUSTER observations reported by |72]. The current sheet (tangential discontinuity) deflects the upstream reflected 
ion component, channels the ions back and focusses them along the current sheet. This produces a highly localised 
hot ion population which starts expanding, blowing off the upstream flow to both sides, deflecting the flow from its 
original direction, and piling the upstream plasma up to the sides of the current sheet where it causes new shock 
waves. 
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Tangential discontinuities have the property that they are convected with the upstream flow. The model is shown 



schematically in Figure 36 There is no flow across the discontinuity. Thus, in the upstream frame of reference the 
velocity component Vn perpendicular to the discontinuity is strictly zero; in the shock or laboratory frames it is the 
same on both sides of the discontinuity. The normal component of the upstream magnetic field i?„ = vanishes in all 
frames. The flow velocity Vt and the upstream magnetic field Bt are both tangential to the discontinuity to both 
of its sides but can be rotated by an arbitrary angle and can have different magnitudes on both sides, while the total 
pressure P = NT + B^/2^q is constant across the discontinuity in the co-moving upstream frame. Thus, the density, 
tangential magnetic field, and tangential velocity can change arbitrarily across the discontinuity. The temperature 
change is then fixed by the continuity of the pressure (we do not consider here an anisotropic pressure, however) . 

Any change in the magnetic field, which is a rotation of the field in the discontinuity plane plus a stretching 
of the tangential component, corresponds to an electric current flowing in the discontinuity. (Note that different 
magnitudes of the tangential magnetic fields to both sides of the discontinuity simply mean that the current flowing 
in the discontinuity is not distributed homogeneously over the width of the discontinuity, while a rotation of the field 
across the discontinuity implies an inclined current.) 

In addition, tangential discontinuities have one property that turns out to be of vital importance for the formation 
of a Hot Flow Anomaly. Tangential discontinuities possess a non-vanishing normal electric field component E„ = 
— Vt X Bt. This field vanishes only when Vt and Bt are either parallel or anti-parallel. However, its direction can 
be along the external normal to the tangential discontinuity (pointing away from the discontinuity) or anti-parallel to 
the external normal (pointing toward the discontinuity). In the first case, ions will be removed from the discontinuity 
by this field, while in the second case ions will be returned to the disontinuity. One expects that only in the second 
case a Hot Flow Anomaly is created in the interaction with the shock, while the first case will cause only a minor 
distortion of the shock. Hybrid simulations show that this is indeed the case. 

I. Two-dimensional hybrid simulation. An example is given in Figure [37] which shows the results of a two- 
dimensional hybrid simulation of the interaction of two discontinuities with a shock |143j . The only difference between 
the two discontinuities is that the normal electric fields E„ = — Vt x Bt to both sides of the discontinuity at a; = 25 
points away from the discontinuity while to both sides on the discontinuity at x = 75 it points into the discontinuity. 
The effect on the interaction between the discontinuities and the shock is dramatic. The x — 25 discontinuity causes 
only a minor distortion on the shock while the discontinuity at x = 75 digs a big hole into the shock and at the 
same time produces a violent effect in the upstream region which shows all the signs of an HFA. Note that because of 
technical reasons the simulation was performed assuming that the discontinuity impinges on a perpendicular shock, 
which is a reasonable assumption for the initial time of the interaction, because the discontinuity, being a tangential 
discontinuity and therefore containing in its plane the full upstream magnetic field, first touches the shock in the area 
of the perpendicular shock where the upstream field is tangential to the shock. Only later it is swept over the shock 
surface by the convective flow sitting for the longest time in the larger quasi-parallel region where HFAs are usually 
observed. 

The explanation of the simulation result is quite simple even though the complete physics of the interaction is still 
poorly understood. The quasi-perpendicular supercritical shock reflects ions. When _E„ = points away from the 
discontinuity all reflected ions in the interaction region between shock and discontinuity are channelled away from 
the shock. This has the effect, that shock reformation is inhibited by the discontinuity, and the shock becomes nearly 
stationary with the exception of minor distortion around the intersection point, as is seen in the lower part of the 
figure around x — 25. There is a slight increase in density N at the interaction point in front and a dilution behind the 
shock surface along the direction of the discontinuity where all the fluctuations are damped out by the discontinuity, 
and the shock front becomes narrow, weekend, and stable. A similar effect is seen in the temperature T. 

The more interesting violent distortion in the case En = E^ around x = 75 is caused by the channelling of 
the shock-reflected ions into the discontinuity. Inside of the tangential discontinuity the ions become approximately 
unmagnetized due to the presence of the current (which generates the discontinuity) and the correspondingly low 
internal magnetic fleld. The isotropic ions accumulate and cause a hot dilute ion cloud that compresses the shock in 
the downstream direction, excavate the upstream flow from the region, deflect it into tangential direction, and cause the 
density and magnetic fleld to pile up at the boundaries of the discontinuity. The shock is hereby violently deformed 
both in the upstream and downstream directions. In addition, upstream shocks are formed at the discontinuity 
boundaries (in the non-symmetric case only one shock boundary would form) ; downstream a HFA bubble structure is 
caused with sharp, dense and hot walls, as seen from the figure. These simulations are hybrid and thus not completely 
reliable as the effect of the electrons is not taken into account. Here the electron fluid is simply heated. Moreover, 
they use a perpendicular shock, i.e. the simulation is relevant only for the initial state of HFA formation. 

m. 'Global' hybrid simulations. The question arises whether the HFA survives or is modulated when moving 
into the quasi-parallel shock region. This question has been investigated with the help of a global two-dimensional 
hybrid simulation [93j which takes account of the entire obstacle, in this particular case the magnetosphere in its 
interaction with the supersonic solar wind. In such a simulation, because of the finite transverse size of the obstacle 
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FIG. 37: Two-dimensional hybrid simulation of the formation of a Hot Flow Anomaly (HFA) in the interaction of a tangential 
discontinuity with a quasi-perpendicular shock [after [T43l . The idea is that the tangential discontinuity first touches the shock at the 
quasi-perpendicular site because the upstream magnetic field is contained in the discontinuity plane, initially interacting with the 
shock-reflected ions. Two discontinuities are used. One with normal electric field E„ = pointing toward the discontinuity, the other 
with En = En pointing away from the discontinuity. Only the former generates an HFA by capturing the shock reflected ions. The other 
causes only a minor disturbance of the shock. The figure shows the density N and plasma temperature T in the simulation plane (a;, z). 
Top: Stacked profiles of density and temperature. The _E^-discontinuity produces a major disturbance in N and T at the shock surface, 
consisting of two dense and hot walls surrounding a diluted region. Bottom: Iso-contours of N and T showing the extension of the HFA 

disturbance on the shock profile. 



(the blunt magnetosphere) , the bow shock assumes its natural curvature around and distance from the blunt obstacle. 
A tangential discontinuity that is convected toward the obstacle thus starts interacting with the self-consistent shock 
first in the quasi-perpendicular region before reaching the quasi-parallel domain. On the other hand, global simulations 
suffer, in addition to being hybrid only, from the largeness of the scale that can be resolved. In these simulations 
the cell size is ~ lAi x lA^. Moreover, global simulations are not completely collisionless, which introduces some 
non-realistic element through the resistive scale = Xi^J v^jio, where is the numerical collision frequency, and 
T = Lo^^ is the time scale of the variation of the field quantities. 

Figure 38 summarises the results of the global simulation. On the left of this figure a large part of the simulation 
box is shown at time tiOd — 150 when the tangential discontinuity (TD) has arrived in the quasi-parallel shock region. 
Shown is the bulk flow velocity in colour coding. Upstream directed velocities are from white to red, downstream 
directed from green to blue. The two panels in the centre and on the right are enlargements of the white box on 
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FIG. 38: "Global" hybrid simulations of the interaction of a tangential discontinuity in the upstream medium (solar wind) with a 
supercritical shock (bow shock) [data taken from '93] • The obstacle is the magnetosphere. Left: The bulk velocity in colour coding. The 
simulation plane is the noon meridian plane (North-South). The magnetosphere is on the right, upstream is on the left. Originally the 
magnetic field had a shock normal angle 0Bn = 45° , i.e. the upper half of the shock is quasi-perpendicular, the lower half quasi-parallel. 
The tangential discontinuity (TD) is parallel to the field (white fat line in the lower left corner). Approaching the shock from left it moves 
down in the box from upper left corner to its position shown here at time tuid = 150. Behind the TD the magnetic field points into the 
simulation plane. Hence E„ points away from the TD behind and into the TD in front of the discontinuity. The case is non-symmetric. 
Centre: Enlargement of the white box in the panel on the left. Shown is the ion density. Right: Ion temperature in the enlarged box. 



the left as indicated. They shown in higher spatial resolution the ion density and ion temperature, respectively. The 
original upstream magnetic field (in front, i.e. downstream of the tangential discontinuity TD) was directed at a 
shock normal angle &Bn ~ 45° (as shown in the lower right corners of the central and right panels). Thus the original 
shock was quasi-perpendicular in the upper half of the left panel, while in the lower half it was quasi-parallel. Behind 
the TD, i.e. upstream of it the magnetic field points into the simulation plane, thus being rotated by 90° in the 
discontinuity plane. The TD being inclined like the original magnetic field at an angle of 45° arrived at the shock 
from above in the left panel touching it first at its quasi-perpendicular part, and then moved down along the shock 
into the quasi-parallel domain. The first and very interesting observation in this simulation was that during the entire 
pass along the quasi-perpendicular bow shock no HFA was created. This observation is somewhat disturbing in view 
of the simulations by |143j who used a perpendicular shock. Possibly this is because of the asymmetry of the normal 
electric field in the "global" simulation which inhibits accumulation of reflected ions on the backside of the TD. 

When the TD arrives at the quasi-parallel shock a Hot Flow Anomaly (HFA in the figure) is immediately formed. 
This is shown for the time tcuci = 150 in the figure. The left panel plots the local ion velocity. White-to-red colouring 
indicates upstream directed velocities, conversely green-to-blue colouring indicates downstream directed velocities. 
The magnitude of the velocity has been colour coded (as given by the bar on the right of the left panel in relative 
simulation units) . The Hot Flow Anomaly appears only on the front side of the TD and is seen as an extended region 
of fast upstream directed ions. These ions are concentrated and fiowing along the TD, which has caused a distortion 
of the quasi-parallel shock. 

The central and right panels are enlargements of the white box region in the left panel, but showing now the density 
and temperature variations, respectively. In both cases white is low, green and blue high. Again, in the density and 
temperature the HFA appears along the front side of the TD as a very dilute hot plasma region that is bounded by 
dense plasma walls which are approximately parallel to the TD. 

The bow shock is heavily distorted and deformed by the presence of the TD in a way similar to what had been 
found in the local hybrid simulations of [143 , bending into the downstream direction, away from its original shield 
shape. On the downstream side of the TD there is no sign of any distortion of the medium other than that the bow 
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FIG. 39: Variation of the plasma parameters and fields along the trajectory LI of a Active spacecraft passing the HFA in Figure |38| 
[data taken from 193 . The simulated HFA shows all the signs of an observed HFA. 



shock has become a broad transition region of upstream directed velocities with a sharp, dense upstream edge, which 
[93j identify as a "sohtary shock". This sharp edge is also seen in the density and temperature plots. 

Most interesting is the strong effect of the TD on the ion population. Upstream of the TD, i.e. behind it, the 
plasma is cold and dilute and of high downstream directed velocity. Downstream, in front of the TD, the plasma 
is hot and dense, being composed of the hot diffuse ion component that populates the upstream foreshock of the 
quasi-parallel shock. The TD sweeps this diffuse ion component along when moving with the flow, forming a very 
sharp boundary between the two regions. It thus seems that the diffuse upstream ion component at the quasi-parallel 
shock is heavily involved into the formation of the HFA. These ions are absent at the quasi-perpendicular shock. 
However, the mechanism is not clear in this case, and while it is very interesting to see these differences between the 
local and global simulations, the global simulations differ so much from the local simulations that it is very difficult 
to compare them on the same safe grounds. 

In order to compare the simulation results with observation j 93j let a fictive spacecraft fly along the line LI in their 
simulations. Figure[39] shows the changes in the plasma and field parameters along LI. On the left in this figure is 
the ion foreshock region (up to point A) with the usual fluctuations in the plasma and field quantities that are typical 
for the foreshock. The fictive spacecraft encounters the HFA about at point A, detecting a dense wall and strong 
magnetic field that separate the hot HFA interior from the environment. This interior has the property that the flow 
velocity has dropped to zero while the transverse flow velocity Vy has taken over the entire moment of the flow. 
Note that outside the HFA this component was zero. The magnetic fleld inside the HFA is weak, either corresponding 
to a magnetic hole or a current sheet. Leaving the HFA close to the tangential discontinuity TD, density and fleld 
are again enhanced, and the flow returns to its original direction. All this happens only on the front side of the TD. 
After passage across the TD, the flctive spacecraft flnds itself located in the upstream flow domain well outside the 
foreshock without any signature of an HFA. The plasma is dilute, cold and quiet and the magnetic fleld is weak, here. 

Most of these signatures agree quite well with the observations of Hot Flow Anomalies. As we have noted, the 
discrepancy is mainly the boundary property of the TD between the foreshock and the quiet upstream flow. Moreover, 
the HFA appears only on the front side of the HFA lying entirely inside the ion foreshock on the quasi-parallel shock. 
Whether this is an effect of the asymmetry of the normal electric fleld at the tangential discontinuity is an open 
question that cannot be decided at the present time but needs more and specialised simulations. Concerning the 
tangential discontinuity we note that the 'global' simulations seem to indicate that on the perpendicular side of the 
shock the tangential discontinuity simply passes the shock and survives in a modiflcd form behind the shock in the 
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FIG. 40: 'Global hybrid' simulation of the formation of a 'solitary shock' at the bow shock [data taken froml93|. Left: The solitary 
shock as a replacement of the bos whock in the quasi-perpendicular bow shock region appears as a narrow sharp boundary with a new 
downstream region behind it (labelled 'New sheath'). The 'Old sheath' is seen as a dense region. Right: Plasma data taken along the 
fictive spacecraft trajectory indicated in the left part of the figure by L. The solitary shock is shaded, while the new sheath region is light 

coloured. 



downstream flow, 
conclusion. 



This seems to be the case in the upper part of the left panel in Figure 38 which advocates such a 



"Solitary shock" 



[55] , in their 'global' hybrid simulations, discovered a new type of a thin shock, which they call a 'solitary shock' 
(because of its narrowness). The signature of this solitary shock is seen in Figure 38 as the sharp boundary between 
downstream and upstream directed flows behind the tangential discontinuity. 

[93^ investigated this new type of shock more closely, finding that it is completely unrelated to the presence (or 
absence) of the tangential discontinuity being a proper phenomenon of formation of a bow shock in front of a blunt 
obstacle in a magnetised flow (note that the 'global' simulations are not completely coUisionless, even though in 
investigating the 'solitary shock' |93j have tried to determine the dependence of its formation on the resistive scale!). 
Simulations without a tangential discontinuity also show the formation of a 'solitary shock'. 

The condition under that it is produced seems to be related to the direction of the motional upstream electric shock 
field only. If this has a component that is directed upstream (in the shock frame away from the shock), then 'solitary 
shocks' seem to form spontaneously. This motional electric field direction determines the drift and acceleration of 
the shock-reflected ions. Nonetheless, this observation does not yet lead to a complete understanding of the very 
mechanism that leads to the formation of this new kind of a shock that under the condition of the upstream directed 
motional electric field component might evolve and is attached to a bow shock. 

Figure 40 gives an example of the 'solitary shock' taken from the former simulations (still including the TD) plus 
a variation of the plasma and field quantities along the cross section L in the figure. The 'solitary shock' appears as 
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FIG. 41: Variation of the plasma parameters and fields along a direction normal to the 'solitary shock' and a model 'perpendicular 
shock' [from '93 . Indicated are the main differences: the motional electric field and the deflection of the flow which are both in opposite 

directions. 



a narrow line in the density which bounds the new sheath region sharply, separating it from the upstream flow. The 
plasma and field data on the right show the sharpness of the 'solitary shock' (grey shading), the spike in the electric 
field related to it, and the fiuctuating new sheath formed behind it all being separated from the old sheath. 

Simulations performed for a perpendicular shock without the presence of a tangential discontinuity but otherwise 
under exactly the same conditions with just different directions of the motional upstream electric field are compared 
in Figure|4l] The 'solitary shock' on the left behaves very similar to the perpendicular shock on the right. The 
two exceptions are the initial assumptions that the direction of the electric field in the 'solitary shock' is opposite 
to the direction of the electric field in the perpendicular shock. Related to this is the opposite direction of the flow 
deflection in Vy. In addition, the solitary shock transition is substantially narrower than the transition region of the 
perpendicular shock. 

We do not intend to discuss these observations/simulations in further detail as neither additional simulations or 
investigations in depth are currently available, nor is a comparison at hand to observations in real shocks like Earth's 
bow shock. Moreover, local simulations have demonstrated the importance of electron dynamics for shock formation 
and reformation. This dynamics is completely ignored in the global hybrid simulations. It is therefore not certain 
whether this finding will survive the experimental test in space or full particle PIC simulations with large mass ratios. 
Such 'solitary shocks' should, if they exist, be found at the bow shock, because the bow shock is a three-dimensional 
structure. Hence, at some location on its surface the motional electric field should have the correct direction to 
generate a 'solitary shock' . So far no indication of such a phenomenon has been seen, a lacking result which might 
be due to having escaped recognition. 



V. THE DOWNSTREAM REGION 



All shocks posses a downstream region which is located between the obstacle and the shock or the shock and the 
slow stream which has been overturned by a fast stream. This downstream region belongs to the shock transition 
with its properties being determined by and large by the shock. The downstream region has two boundaries, the 
shock and the (blunt) obstacle in the wide sense of the meaning, being a solid body, an unmagnetized atmospheric 
gas, a magnetosphere, a magnetic piston, the driver of a blast wave, the volume of a Coronal Mass Ejection (CME), 
or simply the slow stream fiow. In each case the properties of the downstream region will be different because of 
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the differences in the downstream boundary conditions at the obstacle. Therefore, the properties of the downstream 
region cannot be considered without including the obstacle. Nevertheless, a few general conclusions can be noted as 
far as they concern the effect of the shock boundary on the downstream region. 

It is clear that in the average downstream of the shock, whether quasi-perpendicular or quasi-parallel, the flow is 
decelerated to Mach number M <l, deflected to flow around the obstacle (independent on whether it can penetrate 
the obstacle in some cases up to a certain percentage Hke in the presence of a magnetosphere and diffusive processes 
or reconnection), the magnetic field direction changes, the plasma and magnetic fields are compressed, and the 
temperature and pressures are increased. Moreover, the pressure and temperature anisotropies which were very 
moderate upstream of the shock change across the shock because of some adiabatic heating in the shock and other 
effects. The precise processes and the degree to what this heating, change in anisotropy and changes in the parameters 
are predicted for the different kinds of shocks, shock normal angles, and Mach numbers are not precisely known, but 
the general trend is well described by the above list of changes. Thus, when asking for the effects of the presence of 
the shock on the downstream region we are less interested in these average quantities than in the dynamic behaviour 
of the downstream region as a function of the nature of the shock. 

Th dynamic behaviour meant is the turbulent behaviour of the plasma downstream of the shock. What re the 
turbulent properties? Can turbulence develop downstream of a shock? Is the turbulence, if it exists, a consequence 
of the presence of the shock, i.e. is it generated by the shock, or does it evolve locally in the downstream region, 
and the shock provides just the background properties for it such as, for instance, plasma pressure anisotropies? Are 
there differences in the properties of the turbulence between quasi-parallel and quasi-perpendicular shocks? Of what 
modes is this turbulence made of? What is the shape of the possible turbulent spectra? These axe some of the main 
questions that arise. 

Few of these questions can currently be answered less, because there would be insufficient observational material 

but, because turbulence theory is still incomplete and has not been developed for regions of such a limited transverse 
extent like the transition region between a shock and the obstacle in a high Mach number supercritical plasma flow. 
As usual, the bet investigated example is the Earth's magnetosheath plasma. But even here no consensus has been 
reached so far about the state of the turbulence and whether the state of the magnetosheath plasma can at all be 
called turbulent or not. The evolution of turbulence requires time. In a streaming plasma like the magnetosheath, 
which is one a few Earth radii in diameter, there is barely sufficient time available for the turbulence to develop up to 
a stationary state. In particular waves of large wavelength where the energy input is expected to occur at the shock, 
might not have enough time to cascade down to form a turbulent spectrum, hence any spectrum will be cut at wave 
numbers corresponding to the scale of the inverse diameter of the magnetosheath. 

Another reason is that the downstream flow is highly inhomogeneous on all scales larger than a few Ai. This 
implies that the theory of homogeneous turbulence does not apply, and it is not known very well what to make out 
of it, because at the large scales the effect of the presence of the downstream obstacle boundary cannot be neglected 
anymore, which inhibits to make general conclusions about the turbulent state of the downstream plasma. 

A. Sources of downstream turbulence 

Let us - for the moment - assume that the downstream plasma is indeed in a turbulent state, an assumption we may 
correct later. What are the sources of this turbulence? 



1. Upstream waves 

The first source that may contribute to the downstream turbulence are the upstream waves. If the upstream waves 
vanquish the shock barrier, they can enter the downstream region, interact with the downstream plasma population 
until either being damped or amplified; they can decay into other waves which the downstream plasma allows to 
propagate and, if of sufficiently large amplitude, they can cascade down into a broad spectrum of turbulence. As we 
have seen, there are several types of upstream waves: periodic ultra-low frequency waves, shocklets, huge pulsations, 
several types of whistlers from upstream propagating Alfven ion-cyclotron waves to phase-locked whistlers standing in 
front of the shock and its elements, the pulsations; there are waves which are generated in the ramp, high-frequency 
electron waves, BGK modes, and the waves in the feet of the quasi-perpendicular shock and, since the pulsations have 
turned out to behave very similar to quasi-perpendicular shocks, the waves in the feet of pulsations, viz. Buneman and 
modified-two stream modes, depending on the conditions and plasma composition. Of course, in the quasi-parallel 
shock upstream region the plasma conditions are modified by the presence of the diffuse ion component, which will 
change the conditions in the foot region of the pulsations. These questions have not yet been attacked properly, at 
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FIG. 42: Connection between upstream (foot) whistlers and downstream (transmitted/excited) long duration whistlers behind the 
(quasi-perpendicular) shock [data taken from 106 . The whistlers form a well distinct band in the downstream region. Top panel: An 
early dynamic spectrum observed by the IMP 6 spacecraft in 1972. The shock is quasi-perpendicular in this case. The vertical lines are 

time marks. The horizontal line upstream at 300 Hz is the electron gyro-frequency, downstream of the shock it fluctuates strongly 
around 1 kHz. Waves at the local whistler frequency appear in the shock foot. Some of these waves are whistlers. Behind the shock only 
the higher frequency part of this whistler band survives. Bottom left: The downstream magnetic spectrum showing the spectral power 
density above instrumental noise level. The two curves correspond to few second averages, the upper curve shows the 30 ms peak values 
within the average measuring time. Both curves show the presence of whistler waves at a fraction of the electron cyclotron frequency. 
Peak values are up to three orders of magnitude higher indicting the high variability of the downstream waves. Bottom right: The 
occurrence frequency of whistler versus frequency. The peak is at about 25% of the electron cyclotron frequency. 



least not to the extent as they have been investigated for quasi-perpendicular shocks [as was done in the papers by 
[751175]. 

a. ULF waves. Quasi-parallel shocks possess a large spectrum of upstream waves that propagates mostly in the 
upstream direction away from the shock, when seen from the plasma frame but is convected into the shock by the 
flow. These ultra- low frequency waves propagate in the fast, Alfven-ion cyclotron, and whistler modes. 

In the high-Mach number upstream flow they have no chance of propagating far ahead of the shock. However, 
arriving at the shock or passing it, the downstream convection becomes sub-magnetosonic and the wave speed has 
a chance to compensate for the convection. In this case the passing waves will accumulate in the shock transition 
just behind the ramp. However, the ramp is nonstationary. It reforms quasi-periodically and jumps ahead upstream. 
In a real system like that of the bow shock where the shock is found at about the same position this reformation 
implies that the shock ramp oscillates back and forth around its nominal position. It consists of an accumulation of 
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pulsations, and when a new pulsation takes over to become the shock, the old pulsation is expelled downstream with 
all the upstream waves that have accumulated in it and are attached to it. 

Reformation of a quasi-parallel shock thus implies that quasi-periodically a pulsation and a bunch of upstream 
waves is added from the shock to the downstream medium where they start their own shock-independent life. They 
then contribute to the downstream wave population that is in some kind of turbulent state and to wave dynamics, 
may participate in the resonant or nonresonant interaction with the particle populations or with other waves, or they 
contribute to the turbulence via the turbulent cascade. 

b. Whistlers. The higher frequency branch of the upstream waves propagates in the whistler mode. An example 



of an observation of (what is believed to be) whistler waves across a (quasi-perpendicular) shock is shown in Figure 42 



The upstream whistlers in the foot of the shock form a broad intense fluctuation band well below the electron cyclotron 
frequency (the about stationary horizontal line upstream of the shock ramp). There is a continuous connection 
between this whistler band across the shock to downstream, now with the low frequency part of it cut off. Only 
the higher frequencies, which are now found far below the electron cyclotron frequency, survive in the downstream 
region. It is, however, not known whether these waves have propagated from upstream across the shock (here the bow 
shock) to reach the downstream region or whether they are excited about locally in the downstream region (here the 
magnetosheath) . The continuity of the whistler band from upstream to downstream suggests the former: probably 
most of the downstream spectrum leaks in from upstream. If this is the case, then the lower cut-off of the spectrum 
is due to the propagation direction of the whistlers. 

The upstream whistler waves are attached to the leading edge of the shock, respectively the pulsation - in the case 
of a quasi-parallel shock. We already saw that the arrival of a new pulsation in the quasi-parallel case damps these 
whistlers out, in the simulations apparently killing them completely. The ejected 'old shock front' that has been the 
'old pulsation' shifts to downstream, but the whistlers which had been attached to it, have gone in the simulation. 
Hence, any whistlers which occur in the downstream region should have been produced in a different way, either by 
the downstream particle population, by wave-wave interaction, or in some way by the cascading of the large amplitude 
pulsation into smaller amplitude waves which build up a turbulent spectrum. Neither of these mechanisms has so far 
been explored. On the other hand, if the damping of the whistlers is not complete - which is highly probable in reality 
-, the pulsation which is ejected to downstream will carry the surviving whistlers along with it and add them to the 
downstream wave population. In any case the continuity of the whistler band across the shock in Figure 42 suggests 
that the downstream whistlers are not independent of the upstream whistlers, at least at the quasi-perpendicular 
shock. If this is the case, it should not be much different in the quasi-parallel case in particular that everything in 
the simulations points on the quasi-parallel shock behaving quasi-perpendicular the closer to the shock ramp. 

Local generation of whistlers downstream of the shock requires large electron temperature anisotropies Te±/Te\\ > 

1 and electron resonant energies £e,rcs > I'^^eV^e from linear theory, where Vac = B /y/JIoNrrie is the electron 
Alfven velocity, here taken in the downstream region. Figure|44] below indeed indicates an anisotropy in the electron 
temperature downstream of the shock. It also indicates the presence of counter streaming electron beams in the 
downstream region. Both, the anisotropy and the beams, might contribute to the excitation of whistlers locally. 
Whether this is the case has not yet been checked numerically. 

c. High frequency waves. Plasma waves like Langmuir waves, Buneman modes and ion-acoustic waves have a 
rather different behaviour at the shock. The sharp cut-off of Langmuir waves at the local plasma frequency prevents 
them from passing across the steep density increase at the shock ramp. Only such Langmuir waves can be convected 
downstream which are generated at the top of the density overshoot in the shock. Langmuir waves from the electron 
foreshock when being ultimately convected into the shock ramp, assuming that they have not been damped by the 
diffuse electron population in the foreshock, will stop long before arriving at the shock and will being cut off. 

Ion acoustic waves behave differently. As long as they approximately satisfy the linear dispersion relation the 
increase of the density when approaching the shock and in the shock ramp means that the upstream ion acoustic waves 
move down on the ion acoustic dispersion curve, i.e. their wave number decreases and they become long wavelength 
waves which can pass into the downstream region where they ion Landau damping in the warm downstream ion 
population will absorb their energy and contribute to heating the downstream ions stronger than electrons. It is, 
however, not clear what happens to the spiky ion acoustic waves and ion BGK-modes or holes which have been 
observed in the foreshock. Whether these wave survive the passage across the shock cannot be answered yet. 

Figure|43] shows an example of high frequency electric waves downstream of the bow shock as has been measured 
by the CLUSTER 4 spacecraft. The dynamic spectrum shows a weak indication of the plasma frequency in Langmuir 
waves which must be excited locally either as thermal noise or unstably by weak electron beams. The main signals 
are, however, the broadband short time emissions which do not show the slightest indication of a dispersive drift and 
in frequency go even beyond the plasma frequency in several cases. At a plasma frequency of 60 kHz taken from 
the dynamic spectrum, the ion plasma frequency, where the ion acoustic wave branch flattens out, is between 1 and 

2 kHz which corresponds to the intense red line at the bottom of the dynamic spectrum. 
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FIG. 43: Downstream magnetosheath high frequency waves observed by the CLUSTER spacecraft [data taken froml99|. The dynamic 
spectrum shows indication of the plasma frequency and broadband noise of very spiky nature, similar to that observed earher by |105) 

with the ISEE spacecraft. The high-time resolution wave electric field givn in the bottom panel shows that the latter noise is (as 
expected) produced from the occurrence of many bipolar electric field structures which indicate the presence of ion acoustic solitons in 

the downstream region. 



Such broadband signals can only be produced by highly localised wave packets in the streaming plasma. At a 
streaming velocity of ^ 100 km s^^ this implies wavelengths of the structures of the order of a few meters. In fact 
the extremely high time resolution of the wave electric field given in the lower panel confirms the presence of these 
local wave fields. They form (sometimes non-symmetric) bi-polar signals of amplitudes of some ~ 10~^ Vm^^. There 
are ~ 1 or 2 such structures per ms. At a streaming velocity of ^ 100 km s^^ this implies that the structures are of 
transverse size of < 50 m, i.e. they are Debye scale structures of the kind of solitons or BGK modes propagating on 
the ion-acoustic wave branch. These waves and structures may have passed the shock and entered the downstream 
magnetosheath region. The have similar or even slightly large amplitudes than the similar waves in the foreshock, so 
they should have been amplified in the interaction with the shock. It is, however, clear that they contribute to the 
high frequency turbulence in the downstream region. 



2. Shock generated waves 



The shock barrier can, on the other hand, permit a limited number of upstream waves to pass downstream, it can 
refract the waves and it can transform them into other wave modes. In the latter case the shock barrier will itself 
act as a wave generator becoming a source of downstream turbulence. One type of those waves that is known from 
the observation and simulation of quasi-parallel shocks are the large amplitude pulsations (SLAMS). Even though they 
are not generated in the shock, seen from the downstream region the shock is their source because they reach their 
largest amplitudes in the shock from where they spread downstream. 

Other proposals are that the shock itself generates low frequency waves due to interface instabilities. Streaming 
instabihties have been proposed which may arise due to the sheared velocity tangential to the shock. Other proposals 
concern the generation of interface waves due to the interaction between the downstream hot ion component and the 
upstream ion flow overlap in the shock ramp region. These interface modes would then be swept downstream by the 
flow. Typically those waves have long wavelength parallel to the shock surface and short wavelength of the order of 
the shock width normal to the shock. Along the shock they are included into ripple formation. Quasi-parallel shock 
reformation by pulsations is one kind of those oscillations of the shock surface. 

The shock transition is of sufficient width for short wavelength electron waves to be excited in the shock by the 
current flowing in the shock, and by the free energy contained in the shock heated electron distribution. Such waves 
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FIG. 44: Downstream magnetosheath electron distributions as measured by ISEE 2 [after|33]- The perpendicular distribution Fj_{vj_) 
is a hot symmetric Maxwellian, while the parallel distribution _F||{d||) is cooler, flat top and shows the indication of antiparallel electron 

beams. 



will be radiated downstream from the shock. However, as we have noted above, Langmuir waves can enter only if 
their frequency exceeds the high downstream plasma frequency, while ion acoustic waves excited in the shock can be 
convected downstream. 

3. Passing particles 

Part of the downstream turbulence has its origin in the upstream and shock generated waves. In addition, however, 
the free energy stored in the downstream particle distribution is another source for the excitation of downstream 
turbulence. In the first place it is the energetic ions which carry most of the free energy. They are shock-reflected 
either in the quasi-perpendicular shock or in parts of the quasi-parallel shock. After having been accelerated they 
can pass the shock ramp and move downstream where they occur as the energetic gyrating ion component which 
causes a high perpendicular excess in the ion pressure and thus a pressure anisotropy Pij_/Pi\\ > 1. These ions excite 
propagating waves on the transverse electromagnetic ion cyclotron branch with frequencies below the ion cyclotron 
frequency w < uJci and phase velocities uo/k < Va- This kind of anisotropy is largest close to the shock where the 
gyrating ion component has not yet merged into the plasma background. Further down, the plasma flow between 
shock and obstacle further increases the anisotropy adiabatically. Under certain conditions zero-frequency mirror 
modes can then be excited which require that P^±/Pii\ - 1 > l3~'^. Both types of waves have been observed, indeed. 

Electrons, on the other hand, do not appear as a transmitted gyrating component in the downstream plasma. 
They are heated and accelerated in the shock and are seen downstream as a single hot electron component of lesser 
temperature than the ions. Reasons for this difference in temperature have been given above. Accordingly, electrons do 
in the average provide a warm background component of little anisotropy. This might, of course, locally be incorrect, 
when current sheets are formed as the consequence of shock behaviour or by internal processes in the downstream 
plasma like acceleration in approaching current sheets, current dissipation processes, the internal processes acting in 
mirror modes, and circumstantially sometimes even reconnection in narrow current sheets. Figure[44] shows the typical 
electron distribution in the downstream region as had been measured by the ISEE spacecraft [34^. Two cuts through 
this distribution are given, the perpendicular cut is a symmetric Maxwellian Fe±(v±) of slightly larger temperature 
than the parallel component. The parallel distribution _F'e||(w||) has a flat top and shows an indication of the presence 
of two nearly symmetric but counter streaming electron beams. These measurements are taken close to the shock and 
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FIG. 45: CLUSTER data of magnetic fluctuations downstream of the bow shock [data taken from llir| . The top panel shows the 
normalised fluctuation modulus in the mean magnetic field frame for the four spacecraft which are colour coded. Note that on the 
separation of the spacecraft there is no difference in the fluctuations for any times longer than ~ 10 s. Hence, on these time scales the 
fluctuations are practically coherent, and there is no turbulence neither on this time scale nor on the corresponding spatial scale. The 
lower panels show the power spectrum of fluctuations (bottom left) exhibiting some incoherent power law shape that might indicate the 
presence of turbulence. Based on the spacecraft separation the spectral density in the wave number plane (kx, ky) has been determined 
(bottom right). The power peaks at long wavelengths corresponding to 'mirror modes' as indicated by the superimposed projection of 

the mirror and fast dispersion curves. 



are of interest for what concerns the electron wave dynamics of the near shock downstream region. Such weak beams 
might be involved into the generation of the plasma frequency close behind the shock which is seen in Figure|43] as 
well as for the formation of the BGK modes. They are one of the sources of the high-frequency turbulence in the 
downstream region. 

The formation of the particular form of the electron distribution shown in the figure has not yet been satisfactorily 
understood. Electrons are heated in the shock mostly parallel to the field while also being energised adiabatically 
in the perpendicular direction in the compressed magnetic fields when crossing the shock. Electron distributions of 
this kind can well be responsible for the local excitation of whistlers, as we have mentioned above. However, farther 
downstream of the shock transition the electron distribution looks different, much more isotropic with flat tops in the 
parallel and perpendicular directions and temperature anisotropies near zero [as seen in the CLUSTER measurements 
analysed bvFfi]. 

In the following we briefly discuss some observations of the downstream turbulence and compare them with simu- 
lations of its generation. 

B. Downstream turbulence 

Let us now check whether in the downstream region turbulence can evolve and which are the dominant modes that 
can be observed. This can be done in two ways, either by referring to observations or by referring to simulations. 
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FIG. 46: Left two panels: Dispersion relation of downstream low frequency waves in two-dimensional hybrid simulations at Mach 
number A4a = 3 [data taken from 62:. Shown is the normalised power in the Bz fluctuations plotted versus wave number and frequency. 
On the left is the dispersion relation in the shock frame indicating the two polarisation branches of the waves, on the right is the wave 
power in the {kx, fcy)-plane. The waves are mainly in the A+ mode (directed along the shock normal upstream on the plasma frame) 
splitting into propagation parallel to the magnetic field B and the shock normal n as seen in the right panel. The waves are thus 
Alfven-ion cyclotron waves into which the upstream waves had been converted when passing to downstream. Right two panels: A 
two-dimensional confirmation of these results for two wave propagation angles d = 0°,30° and two ion populations, protons (H+) and 

Helium nuclei (He++) [data taken from I79| . 



1. Observations 



Figure 45 shows an example of recent CLUSTER observations of magnetic field fluctuations downstream of the Earth's 
bow shock wave in the deep magnetosheath . The top frame is a plot of the normalised fluctuation amplitude for 
a time period of 160 s and for all four CLUSTER spacecraft which were at separation of a few 100 km. The fluctuations 
measured by the different spacecraft are colour coded, and the main field has been subtracted. For our purposes it 
is not necessary to know which colour represents which spacecraft. The important observation is, however, that on a 
time scale longer than 10s there is practically no difference in the signals between the different spacecraft. Thus, on 
this time scale and the corresponding separation distance which maps to the wave number, the downstream plasma 
is not in a turbulent state. Any ultra-low frequency waves having entered the magnetosheath or have evolved in the 
magnetosheath remain as such and have not cascaded down into turbulence while having evolved to nonlinearity, as 
is signalled by the large relative amplitudes. On the other hand, on the much shorter time scales the signals are 
incoherent among the spacecraft, and the high frequency fluctuations might have entered some turbulent state. 

This is reflected by the power spectrum in the lower left panel which evolves into a power law tail at frequencies 
higher than ^ 0.2 Hz. Whether or not this tail signals real turbulence remains to be answered. The time domain 
power read from this flgure is ^ —2.3 ± 0.2. If real, it might indicate fractal-dimensional intermittent turbulence, 
however. 

The panel on the right bottom in the flgure is a reconstruction of the fc-space for the ultra-low frequency waves, the 
only ones for which it can be done on the basis of the spacecraft distances. The spectral power peaks at very small 
values of the components kx, ky typical for long wavelength, but the structure does not lie on the fast or Alfven modes, 
they flt better to the mirror mode. Note also that the shape of the wave region in fc-space is elongated along the 
mirror dispersion curve. The method of determination of these curves depends strongly on the spacecraft separation 



and the related resolution and thus remains unaffected from the short wavelength turbulence. [157] tried to identify 
the contribution of modes to the higher frequency turbulence. Choosing a frequency of 0.61 Hz they found that the 
spectrum at this frequency was composed of the superposition of an Alfven wave, a weak contribution of a slow mode 
and some remainder of the mirror mode. It thus seems that the zero frequency mirror mode slowly decays into other 
wave modes. 

|85j analysed the plasma rest frame dispersion relation of downstream ultra-low frequency waves in dependence 
on the downstream distance from the shock (using CLUSTER observations across the bow shock when the spacecraft 
were on short separation of few 100 km). They conflrmed the flnding that close to the shock the waves seemed to 
be a mixture of Alfven ion cyclotron and mirror modes. Farther away from the shock closer to the obstacle (the 
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FIG. 47: Statistical behaviour of normalised plasma rest frame frequency Lo/u)ci = (i^sc ~ k ■ V)/a;ci (left panel) and wave number 
k\i (right panel) from upstream to downstream in quasi-parallel (solid lines) and quasi-perpendicular (light lines) shocks [data taken 
from (90]. Note the shift in frequency to lower frequencies which is most prominent in quasi-parallel shocks suggesting that a substantial 
part of upstream waves is transmitted with about constant frequency. Since u)ci increases, the ratio will decrease. The wave numbers 

increase instead. 



magnetosphere) the mirror mode dominated over the Alfven ion cyclotron mode. |90j extended this analysis to many 
observed events for the regions downstream of the quasi-parallel as well as quasi-perpendicular bow shock. They found 
a downshift in the frequency oj /ud normalised to the local cyclotron frequency and a decrease in kX^ from upstream to 
downstream in both cases (see Figure 47 1 . Generally, the downstream region of a quasi-parallel shock is more disturbed 
than the downstream region of a quasi-parallel shock, i.e. the number of fluctuation events is larger. Moreover, from 
upstream to downstream the direction of wave propagation changes from along the background (main) magnetic field 
to perpendicular to the magnetic field, upstream being in the fast-Alfven ion cyclotron mode, downstream showing 
properties of the slow or mirror modes. From this it might be concluded that the upstream waves (mainly in the 
case of a quasi-parallel shock) are not transmitted downstream [90]. The wave vector direction downstream seems 
to support this conclusion and contradicts the two-dimensional hybrid simulations of quasi-parallel shocks [55] which 
are shown in Figure [46] These simulations suggest that waves can be transmitted and that some of the downstream 
waves propagate like the upstream waves in the Alfvcn-ion cyclotron mode in upstream direction and some in the fast 
mode in downstream direction. 



2. Simulations 



[62j in their two-dimensional hybrid simulations consider a quasi-parallel shock at different Mach numbers from 
M.A < 2 to Ma > 3. As expected at the small Mach numbers upstream standing whistlers are initially attached to 
the shock. Later the backstreaming ion component adds long wavelength waves to them which are excited by the 
right-hand resonant ion-ion instability and propagate along the magnetic field, for the low Mach numbers they can 
escape upstream, but when the Mach number is large they are convected back to the shock and can be transmitted 
downstream. During the transmission they convert into upstream propagating (directed along the shock normal) 
Alfven-ion cyclotron modes. Because of the refraction of the magnetic field their wave vector assumes a substantial 
perpendicular component (Figure 46 1 . These results have essentially been confirmed by two-dimensional simulations 
at quasi-perpendicular shocks also shown in Figure[46] 

To study the transmission, |121| performed one-dimensional high-Mach number hybrid simulations at a nearly 
parallel &Bn — 5° shock. These authors wanted to infer about the role of the shock transition region in the transmission 
process. They found that the shock transition plays a role at medium super-critical Mach numbers Ma < 8 similar to 
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FIG. 48: Shock generated downstream waves as inferred from one-dimensional hybrid simulations at high Mach numbers [data taken 
from [121]. Left: Wave profiles of the negative helicity component of the magnetic field By for high Mach number Ma = 13.25 and two 

different regions 400 < x/Xi > 800 and 1000 < x/Xi < 1400 in the downstream frame as function of time (time is upward on the 
ordinate). The first domain contains the shock location. It is seen that the waves are generated in the shock, are moving in upstream 
direction but much slower than the shock such that they effectively are moving downstream toward the obstacle. They are large 
amplitude and practically undamped. The lower panel is taken in the downstream region. The waves have flattened and attain a longer 

wavelength. Some shorter waves have been produced such that the wave profile is less sinusoidal but still exhibits coherence. Right: 
Shock interface generated wave power spectra as function of wave number for three Mach numbers. The wave power increases with Mach 
number. At the same time the interface waves assume longer wavelengths with the spectral peak shifting to small k. Note also that for 
large Ma ^ 5 the spectrum has an exponentially decaying tail towards larger k. Hence there is no remarkable cascading. 



what had been described by [^2] (who did not see the formation of pulsations because of the small number of diffuse 
ions in their simulations) . For larger Mach numbers right hand polarised waves of both positive and negative helicity 
were found behind the shock ramp which could have been excited by the resonant ion-ion and the non-resonant right 
hand ion-ion instabilities in the very shock transition region. Linear theory predicts such a mechanism in the overlap 
region of the incoming fast ion beam and the dense downstream slow ion component in the very ramp. Here the 
relative velocities are large between the two ion components, and the ambient downstream ion density is high, as it 
should be for the instability to work. Hence the shock transition behaves actively in generating waves in addition to 
those which are already present. At high Mach numbers these waves have larger wavelengths and, because their phase 
velocities are large, they are weakly damped and can survive when propagating downstream. Downstream they may 
undergo an inverse cascade to generate longer wavelength waves in the course of a parametric decay instability. In the 
very high Mach number regime Ma > 10 no damping of these waves has been observed. This implies that for such 
high Mach numbers the waves will survive undamped over the entire downstream transition region between shock (for 
instance the bow shock) and obstacle (for instance the magnetosphere) . One would expect to observe monochromatic 
waves all over the transition region in this case. 

Figure 48 shows the negative helicity waves seen in this simulation at high Mach number. Clearly the waves are 
generated in the shock ramp as is visible from the upper left panel. The waves assume quite large amplitude and are 
moving upstream in the downstream frame, however much slower than the shock. They are thus flowing downstream 
toward the obstacle. With downstream propagation they broaden and attain longer wavelengths. The spectrum on the 
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right shows that the wave power and wavelength increases with Mach number. At large Mach number the spectrum 
exhibits an exponentially decaying tail which suggests that some shorter wavelengths are generated but that the waves 
do not participate substantially in a cascade. The effects observed in the latter simulations were not present in those 
of [55^ because of their much larger shock normal angles and further restrictions on the simulations. It thus seems 
that quasi-parallel shocks can themselves produce waves in their transition regions. It is, however, difficult to identify 
these waves in real observations if waves from upstream are mixed in into the transmitted waves. This will necessarily 
be the case for more oblique shocks that generate large amplitude pulsations upstream which themselves become the 
shock ramps and during reformation expel the old pulsation-old ramp downstream where it becomes downstream 
turbulence, respectively it becomes large amplitude coherent downstream waves with broad low frequency spectrum. 
In addition, whistlers might be attached to these pulsations and escape from them simply because the downstream 
Mach number is small enough for letting them go. Hence, the downstream turbulence is rather in a mixed state 
of some large amplitude waves than in a turbulent state, with the mixed state consisting of pulsations, detached 
whistlers, and interface waves from the shock transition, as well as waves which are excited by the instability of the 
diffuse downstream ion component in the downstream plasma background. These downstream diffuse ions should 
as well interact with the downstream wave component, either damping them, or exciting them, or stimulating their 
decay in a nonlinear wave-particle interaction process. Turbulence, on the other hand, will hardly have time to evolve 
in the narrow transition region between the shock and the obstacle for high Mach numbers, because the transition 
time for the waves will in general be shorter than the decay time of the waves or the typical time to become involved 
into a turbulent cascade. Therefore any conclusions about real turbulence must be taken with care. 

Another interesting conclusion that can be drawn from these considerations and simulations concerns the possibility 
of mirror modes near the shock. The simulations have all shown that the release of pulsations from the shock to 
downstream is accompanied by a sequence of large amplitude magnetic fluctuations and magnetic holes with the 
latter containing groups of hot nearly isotropic ions which compensate for the pressure balance. These holes are 
similar to mirrors but are in fact generated not by the mirror instability but by the shock process of accumulation of 
pulsations and trapping of ions between the pulsations. They are in fact enforced structures that have not evolved 
through an instability but in the violent impact process of an upstream pulsation onto the shock ramp. Their 
observation downstream of the shock implies that the pulsations can survive for quite a long distance downstream. 

We are not entering into any discussion of any processes in the downstream region farther away from the shock 
close to the obstacle. Such processes are strongly affected by the nature of the obstacle and have relatively little in 
common with the shock which for them provides an input source and a boundary condition. The latter consists for 
any waves generated near the obstacle in the inhibition of propagation across the shock into the region upstream of the 
shock. As we will see in the section on particle acceleration such a propagation is possible only for fast particles and, 
of course, for high frequency radiation. Fast particles are produces in almost all shocks to some degree. However, in 
non-relativistic shocks no efficient mechanism is known that would be capable of producing electromagnetic radiation 
downstream of the shock. 



VI. SUMMARY AND CONCLUSIONS 

Supercritical quasi-parallel shocks behave completely different from quasi-perpendicular shocks. They are subject to 
a very high fluctuation level with the shock ramp often very hard to distinguish from the environment in observations. 
They possess extended foreshocks which are divided into electron and ion foreshocks. Electron foreshocks contain a 
shock-reflected and accelerated upstream electron population that is flowing in upstream direction. These electrons 
are believed to have their origin at the quasi-perpendicular part of a bent shock. So far it remains unclear of how these 
electrons are generated and reflected. The shock potential is attracting and not repulsing for electrons. Therefore, 
the mechanism of generating reflected electrons is probably kind of a mirror mechanism in converging magnetic 
fields. Reflected electrons return into the upstream flow along the tangential fleld line as narrow fast electron beams. 
These can generate a spectrum of plasma waves via several instabilities but can survive due to the combined action 
of convection and escape from the self-excited wave population. Probably only part of the electron distribution is 
subject to quasi-linear depletion. 

The ion foreshock is less extended upstream. Similar to electrons an ion beam is flowing along the ion-tangential field 
line and is responsible for generation of ultra-low frequency waves. These waves fill the ion foreshock. In addition, the 
ion foreshock is filled with a diffuse energetic ion population which interacts with the wave spectrum. The upstream 
waves are instrumental in shaping the quasi-parallel shock. During their interaction with the upstream ion population 
they evolve into 'shocklets' and large amplitude 'pulsations', so-called 'SLAMS', which have a particular structure 
consisting of upstream leading fronts and downstream trails. In the plasma frame they move upstream at slow speed, 
in the shock frame they are convected toward the shock front, interact with the diffuse ion background gradient in a 
resonant way, steepen, become very large amplitude and, when arriving at the shock ramp, clump together with other 
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pulsations, accumulate, retard to zero speed and form the new shock ramp. They are responsible for quasi-linear 
shock reformation, which is completely different from quasi-perpendicular shock reformation. Quasi-parallel shocks 
reform due to accumulation of pulsations, while quasi-perpendicular shocks reform due to bunched gyrating foot ions. 

During reformation the old pulsation/ramp is expelled downstream. Whistler that had been standing in front of it 
become partially damped. These and the old pulsation form the initial source of downstream turbulence, however, 
other contributions to turbulence downstream are shock-interfac generated waves, leaking whistlers, and locally excited 
plasma waves due to anisotropy in the particle distribution and due to inhomogeneity. 

Tow important facts concerning quasi-parallel shocks should be noted here. The first is that the generation of 
pulsations in the foreshock implies that the upstream magnetic field is turned about into tangential direction to 
the nominal shock surface. This happens on the scale of a few ion inertial length or less. On the large scale the 
quasi-parallel shock remains to behave like a quasi-parallel shock, however, on the short scale the quasi-parallel 
shock becomes locally quasi-perpendicular, a fact that should be very important in considering particle dynamics, 
in particular the reflection and acceleration of electrons. It should be noted in this context that no ion beams have 
been seen in front of the quasi-parallel shock neither in observations nor in the available simulations. This suggests 
that quasi-parallel shocks remain to behave quasi-parallel for ions, which implies that the scale on that the shock 
becomes quasi-perpendicular leaves not enough space for ion reflection - or is strongly fluctuating such that ions 
become scattered but not really reflected. No consensus has yet been reached on this fact. 

The second fact is that the foreshock is populated not by beams but by energetic diffuse ions the density of which 
increases exponentially towards the shock. There is a close relation between the upstream wave growth and this 
density scale. It is most interesting that this diffuse component has maximum density at the shock ramp and is 
otherwise about homogeneously distributed both towards upstream and downstream of the shock. This implies that 
it is generated in the shock transition. 

All these observations leave a large number of question unanswered of which the production of diffuse energetic ions 
is only one. How are these ions generated? How are they injected? Why does the quasi-parallel shock (apparently) not 
generate ion beams, why not on the short scale? How does a quasi-parallel shock look like in two or three dimensions? 
It cannot behave stationary because it reforms continuously. It jumps back and forth around its nominal position that 
is governed by the arrival of new pulsations at the shock. How is the ejection of the old pulsation/ramp to downstream 
going on? We have seen that in between ion holes are formed and the magnetic field shows minima. What kind of 
structures are these? Are they related to mirror structures, as is naively believed, or aren't they rather violently 
produced ion holes instead of growing from an instability? What kind are the different particle distributions just 
upstream and just downstream of the shock ramp? Which waves can they drive unstable? What is the mechanism 
that generates the electron holes and spiky electric field structures that have been observed in quasi-parallel shocks? 
Is the generation of radiation at the shock ramp possible at all? 

The list of these question could be extended. It is clear that their answering requires considerably more observational 
work, mostly on the small scales. In addition, it requires performing full particle PIC simulations with realistic mass 
ratio in two or even three dimensions in order to include the variability of the shock along its surface and investigating 
the reaction of the various particle components to the shock ramp variations. 
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